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Kecorey  Action  Plan.  Lona-TctmRtpot.,'pubEtbtd  in  Match 1994 by  die  B.S.  Amy  Eognea  District  Portland  an!  satniitltd 
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IV 


PREFACE 


Under  authority  of  the  Economy  Act  (31  U.S.C.  1535)  and  Section  632  of  the  Foreign  Assam*  Act 
(22  US.C.  2357),  the  United  States  Agency  for  International  Development  (USAID)  requested  tne 
Department  of  the  Army,  acting  through  the  US.  Army  Corps  of  Engineers  (USACE),  to  prepare  a 
comprehensive  Recovery  Action  Han  (RAP)  for  Mount  Pinatubo  and  subsequent  hydrologic  events. 
The  RAP  is  being  prepared  in  accordance  with  a  Participating  Agency  Service  Agreement  (PAJ>A) 
signed  on  June  18, 1992,  between  USAID/Philippines  and  the  Department  of  the  Army. 

This  investigation,  “Post  Eruption  Hydrology  and  Hydraulics  of  Mount  Pinatubo,  The  Philippines,” 
was  begun  and  the  report  was  prepared  by  Mr.  Hilaire  W.  Peck  and  Mr.  Kari  W.  Eriksen,  US.  Army 
engine  District,  Portland;  MAJ  Monte  L.  Pearson,  US.  Army  Engineer  Waterways  Experiment 
Station  (WES);  and  Dr.  K.  Malcolm  Leytham,  Northwest  Hydraulic  Consultants,  Inc.,  Kent, 

Washington,  during  the  period  June  1992  to  March  1994.  Data  were  collected  and  analysis  tras 

onfoctefl  by  the  authors.  A  number  of  field  trips  were  made  to  the  study  ate.  Mount  Pinatubo,  The 
Philippines,  during  the  study  period. 


This  report  was  initially  published  as  Appendix  A:  Hydrology  and  Hydraulics  to  the  report  entitled 
“Mount  Pinatubo  Recovery  Action  Plan,  Long-Term  Report,”  published  in  Man*  1994  by  the  U.S. 
Army  Fnginrer  District,  Portland,  and  submitted  to  the  Department  of  State  in  March  1994. 


Tins  investigation  was  performed  under  the  direct  supervision  of  Mr.  Ron  Mason,  Chief,  River  and 
rwai  F^ginrrring  Brandt,  and  Mr.  Mike  Ron.  Program  Manager.  US.  Army  Eagmee:  Uistrmt, 
Vottlmd,  and  Mr.  Jerry  Cornell,  Project  Manager,  US.  Array  Engineer  Division,  Pacific  Dcean;  and 
Dec  W.  F.  Marcuson  HI  and  Paul  F.  Hadala,  Director  and  Assistant  Director,  Geotechnical  Laboratory. 
WES.  directly  supervised  MAJ  Pearson. 


At  the  time  of  publication  of  this  report.  Director  of  WE.1  was  Dr.  Robert  W.  Whalin.  Commander 
was  COL  Brace  K.  Howard,  EN.  Commander  of  the  US.  Army  Engineer  District,  Portland,  was 
COL  Charles  Hines,  EN. 


POST  ERUPTION  HYDROLOGY  AND  HYDRAULICS 
OF  MOUNT  HNATUBO,  THE  PHHJPPINS 


1.  INTRODUCTION 


1.1  Authorization 

Under  authority  of  the  Economy  Act  (31  U.S.C.  1535)  and  Section  632  of  the  Foreign 
Assistance  Act  (22  U.S.C.  2357),  the  United  States  Agency  for  International  Development 
(USAID)  requested  the  Department  of  the  Army,  acting  through  the  U.S.  Army  Corps  of 
Engineers  (USACE),  to  prepare  a  comprehensive  Recovery  Action  Plan  (RAP)  for 
controlling  sedimentation  and  flooding  resulting  from  the  June  1991  volcanic  eruption  of 
Mount  Pinatubo,  and  subsequent  hydrologic  events.  The  RAP  is  being  prepared  in 
accordance  with  a  Participating  Agency  Service  Agreement  (PASA)  signed  on  June  18,  1992 
between  USAID/Philippines  and  the  Department  of  the  Army. 

1.2  Purpose  and  Scope 

This  appendix  to  the  Long-term  Action  Report  is  to  present  hydrology  and  meteorology 
pertinent  to  the  design  of  measures  to  address  long-term  flooding  and  sediment  control 
measures  for  all  eight  major  river  basins  impacted  by  Mount  Pinatubo. 

U  System  of  Units 

The  hydrologic  output  from  this  study  is  reported  in  SI  units.  Volumes  are  most  often 
reported  in  cubic  decameters  (dam1).  One  decameter  is  equal  to  10  meters;  therefore,  1  dam’ 
is  equal  to  1,000  m’.  Table  1.3.1  provides  SI/English  conversions  for  all  SI  units  used  in 
this  appendix. 
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2.  REGIONAL  ANALYSES 


2.1  Study  Area 

Mount  Pinatubo  is  located  approximately  100  km  northwest  of  Manila  in  the  Zambales 
Mountains  on  the  west  coast  of  Central  Luzon.  The  eruption  of  Mount  Pinatubo  in 
June  1991  deposited  enormous  volumes  of  easily  erodible  fine-grained  pyroclastic  material  on 
the  flanks  of  the  mountain.  Debris  flows  and  shallow  flooding  worsened  by  blockage  of 
natural  drainages  have  caused  significant  economic  damage  and  loss  of  life. 

Mount  Pinatubo  is  drained  by  eight  principal  river  systems.  Clockwise  from  the  north,  they 
are: 


O’Donnell-Bangat  Rivers  (tributary  to  the  Tarlac  River) 

Sacobia-Bamban  Rivers 
Abacan  River 
Pasig-Potrero  Rivers 
Gumain-Porac  Rivers 
Santo  Tomas  River 
Maloma  River 
Bucao  River 

Drainage  basins  for  these  river  systems  are  shown  on  Plate  1. 

The  principal  drainages  can  be  conveniently  split  into  two  groups:  west-side  and  east-side 
drainages.  The  west-side  drainages  (the  Santo  Tomas,  Maloma,  and  Bucao)  drain  directly  to 
the  South  China  Sea.  Of  the  east-side  drainages,  the  O’Donnell-Bangat  Rivers  join  the  Bulsa 
River  to  form  the  Tarlac  River,  which  flows  north  to  the  Agno  River  and  thence  to  Lingayen 
Gulf.  The  remaining  east-side  drainages  (Sacobia-Bamban,  Abacan,  Pasig-Potrero,  and 
Gumain-Porac)  are  all  tributary  to  the  Pampanga  River  and  the  Pampanga  Delta  which  flow 
south  into  Manila  Bay. 

Data  used  in  the  analyses  were  largely  obtained  from  the  eight  principal  drainages  and  from 
their  immediate  surrounding  areas.  The  area  considered  in  the  regional  analyses  generally 
lies  between  14*45’  N  to  15*45’  N  and  119*45’  E  to  121*00’  E  of  each  drainage.  A  small 
amount  of  additional  data  was'obtained  from  mote  distant  locations. 

2.2  Physiography 

2.2.1  Pre-eruption  Conditions.  Prior  to  the  June  1991  eruption,  the  peak  of  Mount 
Pinatubo  was  at  1,745  meters  elevation.  The  upper  slopes  of  the  mountain  had  a  dense 
network  of  sleep  and  deeply  incised  drainages.  Above  1,000  meters  elevation,  slopes  ranged 
from  20*  to  65*  (Pierson  et  al.,  op.  cit.)  with  headwater  channel  gradients  in  excess  of 
400  meters  per  kilometer  (m/km).  Channel  gradients  on  the  lower  part  of  the  mountain 


flatten  out  to  about  10  to  20  m/km  at  elevations  of  200  to  300  meters.  Channel  gradients  in 
the  area  of  the  Fampanga  River  and  the  Pampanga  Delta  are  extremely  flat,  dropping  to  as 
low  as  0.1  m/km  to  0.2  m/km. 

The  upper  slopes  of  Mount  Pinatubo  were  generally  densely  covered  by  shrubs  and  tall  grass 
before  the  emption.  The  flatter  and  lower  areas  on  the  mountain  supported  a  variety  of 
crops  including  sugar  cane,  cassava,  and  maize. 

Much  of  the  low-lying  region  surrounding  the  Pampanga  River  east  of  Mount  Pinatubo  is 
intensively  cultivated.  This  is  one  of  the  Philippines’  principal  rice-growing  areas.  The 
more  deeply  flooded  parts  of  the  delta  are  used  for  fish-farming  and  other  types  of 
aquaculture. 

Little  detailed  information  is  available  on  soils  and  surficial  geology  of  Mount  Pinatubo. 
However,  there  is  widespread  evidence  of  past  eruptions,  including  large  expanses  of  old 
pyroclastic  deposits  in  the  Sacobia,  Abacan,  and  Pasig-Potrero  Rivers  on  the  east  side  and  in 
the  Marella  River  valley  on  the  southwest  side  of  the  mountain  (Pierson,  et  al.,  op.  cit.; 
JICA,  1978').  Analyses  of  available  hydrologic  data  indicate  that  soils  on  Mount  Pinatubo 
are  generally  highly  permeable,  though  exposures  of  hard  rock  have  been  reported  since  the 
eruption  in  the  Gumain  and  Bangat  basins. 

2-2.2  Impacts  of  the  June  1991  Eruption.  The  post-eruption  peak  of  Mount  Pinatubo 
is  approximately  1,600  meters,  a  reduction  of  more  than  130  meters  in  peak  elevation.  By 
necessity  the  hydrologic  analysis  is  based  almost  entirely  on  pre-emption  hydrometeorologic 
data.  The  eruption  of  Mount  Pinatubo  caused  substantial  changes  in  the  hydrologic  regime 
of  the  principal  drainages.  Some  changes  (such  as  reduction  in  infiltration  due  to  ashfall 
deposits)  are  believed  to  be  relatively  short-lived  and  are  not  reflected  in  the  hydrologic 
modeling.  Other  changes,  such  as  in  the  gross  configuration  of  the  drainage  systems,  are 
considered  to  be  permanent  relative  to  the  life  span  of  possible  engineering  measures,  and  are 
reflected  in  the  modeling  effort.  The  principal  hydrologic  impacts  of  the  eruption  are 
discussed  briefly  below. 

Changes  in  Headwater  Tributary  Areas.  The  June  1991  eruption  filled  much 
of  the  upper  portion  of  five  basins  (O’Donnell,  Sacobia,  Pasig-Potrero,  Santo  Tomas,  and 
Bucao)  with  pyroclastic  deposits.  The  drainage  patterns  that  developed  in  these  deposits 
resulted  in  numerous  changes  in  sub-basin  boundaries  within  the  Bucao  Basin  and  resulted  in 
much  of  the  headwaters  of  the  Sacobia  being  captured  by  the  Abacan  River.  However,  the 
Sacobia  beadcut  upstream  and  recaptured  its  pre-eruption  headwaters  plus  approximately  4- 
1/2  km  of  channel  length  at  the  upstream  end  of  the  Abacan  headwaters  (i.e.,  Abacan  above 


1  Japan  International  Cooperation  Agency,  1978,  Plowing  Report  on  the  Pasig-Potrero  River 
Flood  Control  and  Sabo  Project.  Main  Report  78-38-1/6. 


A-3 


the  Gates  of  Abacan).  It  appears  that  in  October  1993  the  Pasig-Potrero  captured  about 
21  km'-  of  the  Sacobia  River  headwaters. 

In  addition  to  the  changes  in  basin  and  sub-basin  boundaries  rest  ing  from  new  drainage 
patterns  through  the  pyroclastic  deposits,  the  eruption  itself  left  a  5.9  knr  caldera  that 
captured  portions  of  the  headwaters  of  the  O’Donnell,  Sacobia-Bamban,  Gumain-Porac, 

Santo  Tomas,  and  Bucao  Rivers. 

With  the  exception  of  the  headwater  drainage  areas  of  the  Pasig-Potrero  and  Sacobia  Rivers, 
hydrologic  modeling  for  post-eruption  conditions  is  based  on  drainage  patterns  and  drainage 
areas  inferred  from  aerial  photography  dated  November  1991  and  October  1992  and  assumes 
the  October  1992  drainage  patterns  are  relatively  stable.  The  portion  of  the  Sacobia  River 
headwaters  that  was  captured  by  the  Pasig-Potrero  in  October  1993  was  determined  by 
District  personnel  overflights. 

Blockages  and  Lake  Breakouts.  In  the  period  immediately  after  the  eruption, 
pyroclastic  deposits  caused  numerous  blockages  of  the  drainage  system  on  the  upper  slopes  of 
Mount  Pinatubo.  The  formation  of  lakes  behind  these  blockages  and  their  subsequent  failure 
contributed  to  debris  flows  along  the  principal  drainages.  Given  the  large  amounts  of 
unstable  material  remaining  in  the  headwater  drainages,  temporary  blockages  and  sudden 
breakouts  of  debris-dammed  lakes  may  be  a  continuing  hazard  in  some  basins.  No  attempt 
was  made  to  account  for  lake  breakouts  in  the  hydrologic  modeling. 

Following  the  June  1991  eruption.  Lake  Mapanuepe  formed  in  the  Santo  Tomas  River  basin 
near  the  mouth  of  the  Mapanuepe  River.  At  the  invert  elevation  of  its  outlet,  Mapanuepe 
Lake  has  a  surface  area  of  approximately  8  km5.  The  formation  of  Lake  Mapanuepe  resulted 
from  a  blockage  of  the  Mapanuepe  River  outlet  caused  by  recurrent  lahars  and  severe 
aggradation  on  the  Marella  River.  This  lake  was  the  source  of  repeated  lake  breakouts  in  the 
months  immediately  l  allowing  the  eruption,  but  a  man-made  outlet  channel  from  the  lake  has 
now  stabilized  the  situation.  Attenuation  of  the  flood  wave  moving  through  Mapanuepe  Lake 
due  to  temporary  storage  of  water  in  the  lake  was  accr  -ted  for  in  the  Santo  Tomas 
hydrologic  model. 

Channel  Degradation  and  Aggradation.  All  the  principal  rivers  draining 
Mount  Pinatubo  have  been  affected  by  extreme  channel  degradation  or  aggradation  at  some 
point  along  their  course.  Channel  degradation  and  aggradation  have  significantly  changed  the 
physical  configuration  of  some  of  these  rivers’  drainage  systems. 

The  lower  reaches  of  the  Pasig-Potrero,  Bamban,  and  Gumain  Rivers  have  aggraded  to  the 
extent  that  they  are  now  perched  between  levees  approximately  three  to  five  meters  above  the 
surrounding  terrain,  and  lateral  drainage  is  unable  tr-  enter  the  main  channel.  These  features 
are  reflected  in  hydrologic  modeling. 
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While  beyond  the  scope  of  this  study,  siltation  of  the  Guagua  River  ard  other  very  low- 
gradient  channels  leading  to  the  Pampanga  Delta  has  resulted  in  severe  flooding  of  low-lying 
land  around  Bacolor  and  San  Fernando,  and  is  reported  to  have  caused  a  general  increase  :n 
the  depth  and  duration  of  flooding  throughout  the  lower  reaches  of  the  Pampanga  River  and 
its  delta. 


Reduction  in  Infiltration.  The  June  1991  eruption  covered  an  area  of 
approximately  4,500  knr  with  airfall  deposits  of  fine  ash  greater  than  5  cm  in  depth,  with 
ash  depths  of  between  about  5  and  50  cm  on  the  principal  drainages  basins  covered  by  this 
study  (Pierson  et  al..  Fig.  2,  op.  cit.).  It  is  generally  believed  that  deposits  of  fine  ash  will 
reduce  pre-eruption  infdtration  rates,  hence  causing  an  increase  in  the  volumes  and  rates  of 
post-eruption  runoff.  In  the  two  years  since  the  eruption,  much  of  this  widespread  deposit  of 
ash  has  washed  off  during  the  monsoon  rains,  or  the  low-infiltration  crusted  surface  that 
forms  on  these  fine  grained  deposits  has  been  broken  up  by  new  plant  growth.  The  reduction 
in  pre-eruption  infdtration  rates  is  believed  to  be  a  relatively  short-lived  phenomenon  and 
was  not  considered  in  the  hydrologic  modeling  presented  in  this  report. 

Loss  of  Vegetation.  The  eruption  of  Mount  Pinatubo  buried  or  otherwise 
destroyed  all  vegetation  over  an  area  of  approximately  300  km'.  Loss  of  vegetation  greatly 
reduces  evapotranspiration  losses  and  eliminates  the  potential  for  both  interception  storage 
and  storage  in  leaf  litter  and  surface  soils  high  in  organic  material.  Loss  of  vegetation  will 
cause  an  increase  in  the  volumes  and  rates  of  post-eruption  runoff.  However,  the  climate  of 
the  area  around  Pinatubo  is  conducive  to  rapid  plant  growth,  and  it  is  generally  believed  that 
re-vegetation  of  all  but  llie  most  unstable  parts  of  the  mountain  will  be  relatively  rapid.  No 
attempt  was  made  to  account  for  loss  of  vegetal  cover  in  the  hydrologic  analyses. 

23  Climatology 

2.3.1  General.  The  Mount  Pinatubo  area,  on  the  west  coast  of  Central  Luzon,  at 

15 '  N  latitude,  has  a  tropical  climate  dominated  by  the  Northeast  Monsoon  during  the  winter 
months  (November  through  May)  and  by  the  Southwest  Monsoon  during  the  summer  months 
(June  through  October)  which  are  the  tainy-season  flood-producing  months.  The  seasonal 
reversal  of  airflow  results  in  a  pronounced  seasonality  in  prevailing  winds  and  rainfall. 

Severe  weather  conditions  (high  winds  and  heavy  rain)  are  associated  with  typhoons,  which 
most  commonly  occur  during  the  Southwest  Monsoon  season  (June  through  October). 

2.3.2  Climatic  Records.  Daily  rainfall  data  were  available  from  15  stations  in  the 
vicinity  of  Mount  Pinatubo.  Of  these,  data  for  13  stations  were  obtained  from  the  Philippine 
Atmospheric,  Geophysical,  and  Astronomical  Services  Administration  (PAGASA).  Data  for 
two  other  stations  were  obtained  from  the  U.S.  Navy  (Cubi  Point  Naval  Air  Station)  and  the 
U.S.  Air  Force  (Clark  Air  Force  Base).  The  stations,  periods  of  record  obtained,  and  other 
relevant  information  are  listed  in  Table  2.3.1.  Data  availability  is  summarized  in  the  time¬ 
line  on  Figure  2.3.1.  The  stations  are  located  on  Plate  1.  A  total  of  about  330  station  years 
of  daily  rainfall  data  was  obtained.  The  majority  of  the  PAGASA  stations  have  a  record 
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length  of  from  15  to  20  years.  Longer  term  records  (30  to  40  years  of  data)  are  available 
from  Iba,  Zam bales;  Cubi  Point  NAS;  and  Clark  AFB.  Data  on  station  elevations  were  not 
available.  However,  reference  to  1:50,000  scale  topographic  maps  published  by  the  National 
Mapping  and  Resource  Information  Authority  (NAMRIA)  indicate  that  all  daily  rainfall 
stations  have  elevations  in  the  range  of  0  meters  to  150  meters  (Ciaric  AFB).  No  daily  data 
were  available  from  stations  at  higher  elevations.  The  elevation  of  Mount  Pinatubo  prior  to 
the  eruption  was  1,745  meters. 

A  network  of  six  automatic  tipping  bucket  rain  gages  was  installed  at  relatively  high 
elevations  on  Mount  Pinatubo  by  the  Philippine  Institute  of  Volcanology  and  Seismology 
(PHTVOLCS)  after  the  June  1991  eruption  PHIVOLCS  high  altitude  rainfall  data  were 
obtained  for  the  period  August  1,  1991  to  January  10,  1993.  Data  from  these  gages  were 
transmitted  by  radio  telemetry  either  to  the  Pinatubo  Volcano  Observatory  at  the  former 
Clark  AFB  or  to  PHTVOLCS  main  facility  in  Manila  in  the  format  of  cumulative  bucket  tips 
recorded  at  pre-selected  time  intervals.  One  bucket  tip  corresponds  to  0.635  mm 
(0.025  inches)  of  rainfall.  The  station  names,  locations,  elevations,  and  other  relevant 
information  are  listed  in  Table  2.3.2.  The  stations  are  located  on  Plate  1.  There  have  been 
considerable  difficulties  in  operating  the  PHIVOLCS  gage  network  and  significant  periods  of 
data  are  missing  due  to,  for  example,  transmission  problems  and  accumulations  of  ash  in  the 
gages.  When  these  gages  are  operational,  data  are  reported  at  least  once  an  hour,  and  more 
generally  at  20-  or  30-minute  intervals. 

Hourly  data  were  sparse.  Gages  at  Iba,  Hacienda  Luisita,  and  Porac  provided  hourly  data, 
but  the  record  at  these  stations  is  described  ^  "fragmentary."  Hourly  rainfall  data  were 
obtained  for  three  storm  events  at  Iba,  three  events  at  Hacienda  Luisita,  and  five  events  at 
Porac.  Information  on  the  data  obtained  is  provided  in  Table  2.3.3.  The  stations  are  located 
on  Plate  1.  Hourly  rainfall  data  outside  the  immediate  study  area  were  available  for  one 
extreme  24-hour  period  from  Weather  Bureau  Technical  Paper  No.  4T,  which  provided 
hourly  data  from  Baguio  City,  located  approximately  140  km  north-northeast  of  Mount 
Pinatubo  at  1,370  meters  elevation  (4,500  feet),  for  a  typhoon-related  event  in 
September  1911  which  generated  the  then  world’s  record  24-hour  rainfall  of  1,168  mm 
(45.99  inches).  The  general  lack  of  good  hourly  rainfall  data  placed  severe  limitations  or. 
hydrologic  modeling  for  the  present  study. 

Daily  pan  evaporation  data  were  obtained  from  three  stations  in  the  vicinity  of  Mount 
Pinambo,  all  operated  by  PAGASA.  The  stations,  the  periods  of  record  obtained,  and  other 
relevant  information  are  listed  in  Table  2.3.1.  The  stations  are  located  on  Plate  1. 


3  Weather  Bureau,  U.S.  Dept,  of  Commerce,  1961,  Generalized  Estimates  of  Probable  Maximum 
Precipitation  and  Rainfall-Frequency  Data  for  Puerto  Rico  and  Virgin  Islands  for  Areas  to  400 
Square  Miles,  Durations  to  24  Hours,  and  Return  Periods  from  1  to  100  Years.  Technical  Paper 
No.  42. 
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Screening  of  Daily  Data.  Daily  rainfall  data  were  screened  to  identify  clearly 
anomalous  (high  or  low)  values  and  to  summarize  periods  of  missing  data.  The  completeness 
of  rainfall  records  is  highly  variable,  with  no  missing  data  reported  at  Clark  AFB  (station 
USA02)  and  over  30  percent  of  the  record  missing  at  Camiling  (Station  R0315).  The  records 
from  Cubi  Point  NAS  (USA01),  Clark  AFB  (USA02),  and  Iba,  Zambales  (D324),  the  three 
stations  with  the  longest  records,  are  all  relatively  complete. 

The  following  points  were  noted  during  screening: 

■  The  maximum  reported  daily  rainfall  of  748.5  mm  for  the  period  of  record 
at  Iba  (D324)  on  July  10,  1980  is  suspect.  No  other  neighboring  stations 
reported  daily  rainfall  greater  than  150  mm  on  that  date.  Examination  of 
hourly  data  for  this  event  showed  that  the  daily  amount  for  this  day  was  in  fact 
78.5  mm,  and  the  record  was  corrected  accordingly. 

■  The  rainfall  record  for  Magalang,  Pampanga  (A020)  is  too  short  to  be  of 
value.  A  longer  record  is  available  from  Bai  Magalang  (R0312), 
approximately  5  km  east  of  the  Magalang  gage. 

■  The  record  for  Bai  Magalang  (R0312)  is  essentially  complete  between 
February  1977  and  January  1992.  However,  the  maximum  recorded  daily 
rainfall  of  103.1  mm  an  November  14,  1977  is  significantly  lower  than  that  at 
other  neighboring  stations. 

■  The  record  from  Camiling  (R0315)  is  too  short  and  too  fragmented  to  be  of 
value. 

■  The  record  for  Mayantoc,  Tariac  (R0316)  is  quite  fragmented  with 
approximately  20  percent  of  the  daily  records  reported  missing.  The  record  is 
particularly  poor  for  the  period  1978  through  1984,  when  27  months  are 
missing,  including  most  records  for  July,  August,  and  September. 

■  The  record  for  Palawig,  Zambales  (R0318)  is  highly  suspect,  particularly  for 
1982  when  no  rainfall  was  recorded  even  though  the  station  was  reported  to  be 
in  operation.  Also,  the  Palawig  station  does  not  report  nearly  as  many  rainfall 
events  during  the  relatively  dry  months  of  January  through  April  as  do  other 
coastal  stations.  Over  the  15  years  of  record,  only  nine  measurements, 
including  trace  rainfalls,  were  reported  in  these  four  months.  It  is  possible  the 
station  is  only  operated  on  a  seasonal  basis. 

«  The  record  for  San  Felipe,  Zambales  (R0319)  is  suspect.  The  maximum 
reported  daily  rainfall  of  162  mm  for  the  period  of  record  is  far  below  that  for 
other  coastal  stations. 
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As  a  result  of  this  basic  screening,  data  from  Camiling,  Palawig,  and  Magalang  were 
dropped  from  further  consideration. 


Double  Mass  Analysis.  Double  mass  analyses  were  conducted  on  daily 
rainfall  data  using  the  following  groups  of  stations: 

Long-term  stations 

Group  1:  Clark  AFB  (USA02);  Cubi  Point  NAS  (USA01);  Iba  (D324). 

West  coast  stations 

Group  2:  Cubi  Point  NAS  (USA01);  Z-NAS,  San  Marcelino  (R0322);  San  Felipe 

(R0319);  Santa  Rita  Elementary  School  (R0320);  Iba  (D324). 

Northeast  stations 


Group  3:  Mayantoc  (R0316);  CLSU,  Munoz  (A0l7);  Hacienda  Luisita  (A016). 

Southeast  stations 

Group  4:  Clark  AFB  (USA02);  Bai  Magalang  (R0312);  Julian  Subdivision  (R0313); 

Masantol  (R0314). 

The  entire  concurrent  rainfall  record  for  each  group  of  stations  was  used  for  the  double  mass 
analysis.  In  computing  cumulative  rainfall  amounts,  days  were  ignored  for  which  data  were 
missing  for  any  station  in  the  group.  The  double  mass  analyses  of  the  long-term  stations 
indicated  that,  with  the  exception  of  an  unusually  high  rainfall  amount  recorded  at  Iba, 
Zambales  (D324)  in  1964,  the  records  for  these  stations  are  fairly  consistent  for  the  period 
1961-1990. 

The  double  mass  analyses  of  the  west  coast  stations  indicated  that,  with  the  exception  of  San 
Felipe,  the  records  are  relatively  consistent  over  the  period  1975-1990.  Removing  San 
Felipe  from  consideration  should  improve  the  results  for  other  stations. 

The  double  mass  analyses  of  the  northeast  stations  indicated  no  obvious  long-term  changes  in 
slope,  although  the  plot  for  Mayantoc  showed  an  unusual  shape  with  a  number  of  small  slope 
changes  and  an  overall  tendency  towards  diminished  rainfall  in  later  years.  Considering 
these  results  and  the  preliminary  screening  noted  earlier,  the  data  from  Mayantoc  (R0316) 
are  highly  suspect. 

The  double  mass  analyses  of  the  southeast  stations  indicated  that  the  data  for  Bai  Magalang 
underwent  several  major  slope  changes  between  1977  and  1991.  Again,  considering  the 
rerilts  of  the  screening  and  double  mass  analysis,  the  data  for  Bai  Magalang  are  suspect. 


For  Masantol  and  Clark  AFB  the  double  mass  analyses  do  not  indicate  any  major  shifts  in 
recorded  rainfall.  At  Julian  Subdivision  there  is  an  apparent  long-term  shift  at  about  1987. 
Comparatively  less  rainfall  is  recorded  at  Julian  after  this  date.  It  is  not  known  if  station 
relocation  or  other  changes  took  place  which  might  explain  this  shift. 

As  a  result  of  the  double  mass  analysis,  data  from  San  Felipe,  Mayantoc,  and  Bai  Magalang 
were  dropped  from  further  consideration. 

Snatial  Variation  of  Daily  Data.  Spatial  variations  in  rainfall  data  were 
examined  by  computing  cross  correlations  between  daily  data  at  all  stations  and  by  examining 
rainfall  amounts  in  individual  major  storm  events.  A  matrix  of  cross-correlation  coefficients 
for  all  stations  is  shown  in  Table  2.3.4.  Cross-correlations  were  computed  for  individual 
pairs  of  stations  using  all  available  data  for  which  daily  rainfall  exceeded  50  mm  at  either  or 
both  stations.  The  maximum  cross-correlation  is  0.53  between  Iba  (Station  D324)  and  Santa 
Rita  Elementary  School  (Station  R320),  which  are  approximately  20  km  apart.  If  stations 
with  suspect  records  are  eliminated  from  the  matrix  (e.g.  Palawig,  San  Felipe,  Bai 
Magalang),  the  cross-correlations  generally  decrease  with  distance  between  stations,  as  might 
be  expected.  The  low  cross-correlations  indicate  substantial  spatial  variations  in  daily  rainfall 
amounts  over  relatively  short  distances. 

Relationships  between  daily  rainfall  amounts  at  various  stations  during  major  storm  events 
were  examined  qualitatively.  Major  storms  in  western  Luzon  are  associated  with  large 
weather  systems  (the  Southwest  Monsoon  or  typhoons)  which  affect  the  entire  area  around 
Mount  Pinatubo.  However,  daily  rainfall  amounts  within  individual  events  show  great  spatial 
variability,  especially  when  rainfall  depths  from  stations  on  the  east  and  west  rides  of  Mount 
Pinatubo  are  compared.  Spatial  variations  of  rainfall  during  Typhoon  Didang  in  May  1976 
appear  to  be  typical.  This  event  produced  the  maximum  two-day  rainfall  depth  of  988  mm 
for  the  40-year  period  of  record  at  Iba.  Significant  rainfall  depths  were  recorded  at  all  other 
west  coast  stations,  although  storm  depths  were  relatively  less  severe.  The  maximum  two- 
day  rainfall  depth  at  Cubi  Point  during  this  event,  for  example,  was  485  mm,  a  depth 
exceeded  six  times  in  Cubi  Point’s  33-year  record.  Comparison  of  rainfall  depths  at  Iba  for 
this  event  with  rainfall  depths  from  the  east  ride  of  Pinatubo  (for  example,  at  Clark  AFB) 
show  even  greater  variation. 

Analysis  of  PHTVOLCS  Data.  Review  of  the  daily  rainfall  records  and 
monthly  totals  recorded  at  PHTVOLCS  stations  indicates  significant  problems  with  many  of 
the  stations.  PHTVOLCS  data  do  not  include  a  flag  to  indicate  missing  data  or  equipment 
malfunction.  Therefore  any  gaps  in  the  record  were  filled  with  a  value  of  0.0.  This  may 
lead  to  inconsistencies  when  comparing  PHTVOLCS  data  to  PAG  ASA  or  other  neighboring 
daily  rainfall  gages.  Following  is  a  review  of  the  available  PHTVOLCS  data. 
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Station  RG-1  —  Mt.  Caudrado:  Data  for  August  and  September  1991  appear  to  be  fairly 
complete.  Monthly  totals  are  not  significantly  greater  than  those  recorded  at  lower  elevations 
but  individual  storm  events  show  increased  daily  precipitation.  Data  beyond  October  1991 
ate  fragmented  and  are  not  useful  in  the  hydrologic  analyses. 

Station  RG-2  —  BUGZ:  Data  for  August  and  September  1991  are  relatively  complete. 
Similar  to  station  RG-1,  the  monthly  totals  are  not  substantially  larger  than  the  lower 
elevation  PAGASA  stations.  Data  from  mid- August  1992  to  mid-  September  1992  may  also 
be  of  some  value.  As  was  the  case  for  RG-1 ,  the  data  for  much  of  1992  were  recorded  as 
no  rainfall.  As  noted  earlier,  there  is  no  way  of  distinguishing  between  missing  data  due  to 
station  malfunction  and  data  that  actually  register  0.0. 

Station  RG-3  —  PI2:  Data  for  all  months  are  unrealistically  low.  The  maximum  monthly 
precipitation  recorded  was  535  mm  in  August  1992.  This  is  significantly  lower  than  the 
804  mm  recorded  at  Bai  Magalang,  for  instance. 

Station  RG-4  —  Mt.  Culianan:  Data  for  all  months  are  unrealistically  low.  The  maximum 
single  day  rainfall  in  1991  was  37.4  mm.  In  comparison,  the  PAGASA  station  at  Iba, 
Zambales  reported  six  events  greater  than  100  mm  during  August  and  September  1991. 

Station  RG-5  —  Gumain:  Rainfall  records  for  the  August  20,  1991  event  appear  to  be 
complete  through  midday  on  August  21,  1991.  With  the  exception  of  this  event,  the  Gumain 
data  are  fragmented,  unrealistically  low,  and  probably  not  of  value. 

Station _RG-6  —  Sacobia:  Data  collection  did  not  begin  until  mid-September  1991.  Rainfall 
totals  for  July  and  August  1992  are  substantially  less  than  those  reported  at  low  elevation 
stations.  The  record  is  fragmented  and  is  not  useful  in  the  hydrologic  analyses. 

The  PHTVOLCS  records  are  too  short  and  fragmented  to  be  used  for  determining  variations 
of  rainfall  with  elevation  on  Mount  Pinatubo. 

Evaporation  Data.  Daily  pan  evaporation  data  were  obtained  from  three 
stations  in  the  vicinity  of  Mount  Pinatubo,  all  operated  by  PAGASA.  The  stations,  the 
periods  of  record  obtained,  and  other  relevant  information  are  listed  in  Table  2.3.1.  The 
stations  are  located  on  Plate  1. 

Most  of  the  missing  data  in  the  evaporation  records  can  be  attributed  to  pan  overflows  during 
wet  weather.  Examination  of  the  available  data  showed  the  record  at 
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Hacienda  Luisita  to  be  fragmented  and  occasionally  reporting  unusually  large  (in  excess  of 
30  mm/day)  daily  evaporation  depths.  Data  from  this  station  were  consequently  dropped 
from  further  consideration.  The  data  from  Magalang  were  also  dropped  because  the  record 
was  too  fragmented  and  too  short  to  be  of  value. 

The  data  from  CLSU,  Munoz  are  believed  to  be  the  most  reliable  pan  evaporation  data 
available.  Missing  periods  in  this  record  were  filled  with  mean  daily  values  for  each  month 
to  obtain  estimates  of  annual  pan  evaporation.  The  estimated  annual  pan  evaporation  is  given 
in  Table  2.3.5. 

2.3.3  Precipitation.  The  rainfall  regime  of  the  Mount  Pinatubo  area  is  highly 
seasonal,  with  a  pronounced  wet  season  from  approximately  June  through  October, 
coincident  with  the  Southwest  Monsoon,  and  a  dry  season  from  approximately  November 
through  May.  Rainfall  amounts  along  the  west  coast  of  Central  Luzon  are  enhanced  by 
orographic  uplift  during  the  Southwest  Monsoon.  Mountains  to  the  west  and  southwest  of 
Mount  Pinatubo,  with  peaks  at  about  1 ,000  meters,  act  as  partial  orographic  barriers  during 
the  Southwest  Monsoon,  with  the  result  that  the  west  and  southwest  flanks  of  Mount 
Pinatubo  probably  receive  somewhat  less  rainfall  than  would  otherwise  be  expected.  Mount 
Pinatubo  itself,  and  the  Cabusilan  Mountains  to  the  south,  present  a  more  significant 
orographic  barrier  during  the  Southwest  Monsoon,  such  that  the  east  and  northeast  sides  of 
Mount  Pinatubo  lie  in  a  rain  shadow,  with  low-lying  interior  areas  east  of  Mount  Pinatubo 
receiving  only  about  half  of  the  total  annual  rainfall  experienced  on  the  coast.  Plots  of  mean 
monthly  rainfall  at  Cubi  Point  NAS  and  Clark  AFB  (Figures  2.3.2  and  2.3.3)  illustrate  the 
seasonal  distribution  of  rainfall  totals  on  the  west  side  and  east  side  of  Mount  Pinatubo, 
respectively.  The  mean  annual  rainfalls  at  Cubi  Point  NAS  and  Clark  AFB  are  about 
3,600  mm  and  2,000  mm,  respectively. 

Maximum  daily  rainfall  amounts  at  stations  in  the  Mount  Pinatubo  area  are  generally  caused 
directly  or  indirectly  by  tropical  cyclones,  which  are  most  prevalent  between  May  and 
November.  Data  available  from  PAGASA  indicate  that  between  1948  and  1991,  an  average 
of  16  tropical  cyclones  per  year  (tropical  depressions,  tropical  storms,  or  typhoons)  affected 
weather  conditions  at  various  regions  in  the  Philippines.  Of  these,  typically  three  or  four  per 
year  affected  weather  conditions  around  Mount  Pinatubo. 

While  large  one-day  rainfall  amounts  can  be  caused  by  the  direct  passage  of  a  typhoon  over 
the  area  (the  maximum  recorded  one-day  rainfall  at  Cubi  Point  NAS  of  442  mm  in  May  1966 
is  believed  to  be  one  such  example),  rainfall  events  of  longer  duration  with  greater  total 
event  volumes  and  comparable  one-day  depths  may  result  from  intensification  (or  surges)  in 
the  Southwest  Monsoon  flow  during  passage  of  typhoons  to  the  northeast  of  Luzon. 

The  Southwest  Monsoon  flow  can  bring  long  periods  of  near-continuous  heavy  rain.  Intense 
localized  rainfall  is  associated  with  intense  convective  activity  in  groups  or  pockets  of  storm 
cells  embedded  in  the  general  southwesterly  flow.  Surges  in  the  southwesterly  flow  due  to 
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the  influence  of  typhoons  may  increase  wind  speeds  by  10  to  20  knots  with  a  concomitant 
increase  in  rainfall. 


Normal  Annual  Precipitation  (NAP)  Man.  Figure  2.3.4  shows  isohyets  of 
mean  annual  rainfall  as  determined  from  assessments  of  available  long-term  rainfall  and 
streamflow  data.  Long-term  rainfall  records  are  available  only  for  elevations  generally  less 
than  100  meters  near  the  coast  and  in  the  interior  lowland  regions.  However,  the  study  area 
includes  much  higher  elevations:  the  pre-eruption  summit  elevation  of  Mount  Pinatubo,  for 
example,  was  1,745  meters.  Streamflow  and  evaporation  records  were  used  to  estimate 
mean  annual  rainfall  at  the  higher  elevations  which  correspond  to  the  watershed  areas  of 
interest,  as  described  in  the  following  paragraphs: 

1)  The  annua]  rainfall  data,  summarized  by  Table  2.3.6,  were  sufficient  to 
determine  an  isohyetal  line  corresponding  to  1,900  mm/yr  of  rainfall  in  the 
lowland  areas  east  of  the  Zambales  and  Cabusilan  mountain  areas,  and  to  show 
that  rainfall  decreases  with  eastward  distance  from  the  mountains.  Coastal 
rainfall  gages  to  the  west  of  the  mountains  indicate  rainfall  of  around 

3,600  mm/yr  to  4,200  mm/yr. 

2)  Annual  yields  from  watersheds  with  streamflow  gages  were  determined  in 
terms  of  annual  depth  of  ninoff  over  the  basin  area  as  summarized  by 
Table  2.3.7.  Data  for  the  Bagsit  and  (upper)  Camiling  stations,  key  stations 
for  rainfall  assessment  because  they  have  relatively  high-elevation  basins,  show 
average  annual  runoff  amounts  of  approximately  3,400  mm  and  3,700  mm, 
respectively. 

3)  Average  basin  rainfalls  were  estimated  by  adjusting  the  runoff  data  by 
1,550  mm  to  account  for  the  estimated  average  annual  evapotranspiration  in 
the  basins.  This  evapotranspiration  amount  was  computed  as  being  equal  to 
about  80  percent  of  pan  evaporation  reported  for  lower  elevations,  assuming 
there  would  be  no  significant  moisture  deficit  in  the  mountains.  In  transposing 
the  estimated  basin  rainfall  amounts  to  the  maps,  it  was  assumed  that  the 
average  basin  rainfall  would  occur  at  the  centroid  of  the  basin,  and  that  rainfall 
would  increase  with  elevation. 

4)  Knowing  the  rainfall  at  the  lower  basin  elevations  from  the  rain  gages,  and 
having  assumed  the  rainfall  for  the  centroid  of  the  basin,  rainfall  at  the  upper 
basin  was  determined  as  the  amount  necessary  to  generate  the  observed  runoff. 

The  analysis  resulted  in  findings  that,  for  mountains  next  to  the  coast,  without  any 
intervening  topographic  barriers,  rainfall  isohyetal  lines  of  5,000  mm/yr  and  6,000  mm/yr 
correspond  approximately  to  elevations  of  750  meters  and  1,200  meters,  respectively.  At 
Mount  Pinatubo,  which  is  separated  from  the  coast  by  other  mountains  and  may  experience  a 
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rain  shadow  effect,  isohyetal  lines  of  4,000  mm/yr  and  5,000  mm/yr  correspond 
approximately  to  elevations  of  500  meters  and  1,000  meters,  respectively. 


The  NAP  map  in  Figure  2.3.4  differs  significantly  from  those  published  with  other  studies  in 
the  area.  Other  available  NAP  maps  show  no  significant  increase  of  rainfall  with  elevation 
over  Mount  Pinatubo.  Most  of  the  annual  rainfall  in  this  region  is  associated  with  the 

f°nr  ?hic[\bringS  moist  air  from  'he  southwest  across  the  South  China  Sea 
nng  the  months  of  June  through  October,  striking  the  west  coast  of  Luzon  An 

!wn,^CatI°n  mo"s°°n  rain  on  to  west  side  of  the  Zambales  and  Cabusilan  Mountains 
due  to  orographic  uplift  is  expected,  and  hence  an  increase  of  rainfall  with  elevation  The 
role  of  orography  in  increasing  monsoon  rainfall  has  been  confirmed  by  meteorologists 
formerly  stationed  at  Cubi  Point  NAS.  S 


i  SrT  T  '  ro  eD  T*  Cau£ing  incrcasin8  rain  with  elevation  is  found 
“f  7Zri  J?ZZ  Cal  PTr?0--  42  ■  ™S  paPer  verified  the  Orphic  component 

Seotemh^qiT  R  "'T  ^  at  City,  Philippines  in 

elevatiorf  1  370  approximately  140  km  NNE  of  Mount  Pinatubo,  at 


.  2  \  A  Air  Temperature.  The  mean  annual  air  temperature  at  Cubi  Point  NAS  fat  sea 

eve  on  the  southwest  of  Mount  Pinatubo)  is  approximately  26-28*  C*  The  annual  average 

“  jese  Cubi  v *  nas  *  3i-4' c  211(1  L  “  4S: r 

Z.  A  plot  of  mean  monthly  temperatures  at  Cubi  Point  NAS  is  given  in  Figure  ”>  3  5 

thc  summit  of  Mount  Pinatubo  (pre-emption  elevation  of  ??«&,) 

ZTT  ,°fbeHa^Ut  V0  10‘  C  C0°Ier  1,1056  a‘  C“bi  poin.  and  Munoz  by  ^ 
consideration  of  adiabatic  lapse  rates.  ^ 

mnmhlv  Wbld  ^  ^  are  available  from  Cubi  Point  NAS  and  annual 

i  /UbUSh£d  by  ,hC  Naval  Weather  S^ice  Environmental 

S  rose  (e  s  lo  6To  T,  15  Sh0Wn  “  FigUre  Z3’6-  1116  numbeis  ^own  the 
-ZLTZlZ-  9  6'  l9'  ’ 9 ’  etC')  reprcs6nt 11)5  Percent  of  time  that  wind  comes  from  the 

.  ?  V2ri0US  122265  °f  Wi"d  SpCed  "*  th=  symbols  that  represent  these 
rnd  inTc  ^^H  r  right  h2nd  C°nJer  0f  Figure  2'3-6-  ^e  percent  time  that 
svmhni  n  %h  Sf>Td  12026  601,165  ^  a  given  dtrearan  can  be  scaled  from  the  appropriate 
rosc  (6f  19  6  P™  of  'he  time  wind  comes  from  the  north^t  and 

n^dl6%t  /  ^  ?”?  W1°d  6011165  from  the  northeast  at  a  speed  between 
1 1  and  16  KTS).  As  mdicated  m  the  center  of  the  wind  rose,  25.9  percent  of  the  time  wind 


1  See  Footnote  2. 

4  Wemstedt,  F.L,  1972,  World  Climatic  Data,  Climatic  Data  Press. 
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speed  is  less  than  4  knots  (KTS).  The  direction  of  winds  less  than  4  KTS  is  not  indicated  by 
the  wind  rose.  As  previously  indicated,  the  winds  come  from  two  predominant  directions: 
from  the  southwest  during  June  through  October  and  from  the  northeast  during  November 
through  May 

2.4  Hydrology 

2.4.1  Discharge  Records.  Daily  streamflow  data  were  available  from  18  stations  in 
the  vicinity  of  Mount  Pinatubo.  Data  were  obtained  from  the  Department  of  Public  Works 
and  Highways  (DPWH)  and  the  Bureau  of  Research  and  Standards  (BRS),  which  currently 
have  responsibility  for  collection  and  dissemination  of  streamflow  data.  Responsibility  for 
streamflow  data  collection  has  changed  hands  several  times  in  past  years  and  most  of  the  data 
obtained  were  originally  collected  by  either  the  Bureau  of  Public  Works  or  the  National 
Water  Resources  Council  (NWRC).  The  stations,  drainage  areas,  periods  of  record 
obtained,  and  other  relevant  information  are  listed  in  Table  2.4.1.  No  daily  streamflow  data 
were  obtained  prior  to  1957,  although  published  records  indicate  that  earlier  data  are 
available  for  some  stations,  as  indicated  in  Table  2.4.1.  The  last  year  for  which  data  were 
officially  published  was  1972.  A  small  amount  of  post-1972  data  were  obtained. 

Streamflow  stations  are  located  on  Plate  I. 

Peak  annua!  flow  data  were  obtained  at  the  18  stations  providing  daily  streamflow  data. 
Station  names,  periods  of  record,  and  other  relevant  information  are  given  in  Table  2.4.1. 

For  tile  most  part,  peak  annual  flows  were  extracted  from  two  publications: 

*  Philippine  Water  Resources  Summary  Data,  Volume  1  —  Streamjlow  and 
Lake  or  River  Stage  Ending  December  31,  *970.  Partial  copy.  Exact  title, 
publisher,  and  date  of  publication  unknown. 

*  Philippine  Water  Resources  Summary  Data,  Volume  2  —  Streamflow  and 
Lake  or  River  Stage  Ending  December  31,  1980.  Republic  of  the  Philippines, 
Department  of  Public  Works  and  Highways,  Bureau  of  Research  and 
Standards,  Quezon  City,  June  1991. 

The  first  of  these  publications  contains  peak  annual  flow  data  for  the  period  of  record 
through  'December  1970.  The  second  publication  covers  the  period  1971  through  1980. 

Some  additional  peak  annual  flow  data  subsequent  to  1980  were  obtained  from  unpublished 
sources. 

Most  of  the  streamflow  data  in  the  study  area  were  obtained  from  staff  gage  readings  that 
were  reportedly  read  two  or  three  tunes  per  day.  Data  at  some  stations  were  collected  using 
an  automatic  water  level  recorder.  Attempts  were  made  to  obtain  these  data  for  selected 
flood  events.  However,  no  data  could  be  obtained  other  than  mean  daily  flows  and  published 
annual  peak  flows. 
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Review  of  Rating  Curves.  Discharge  measurement  data  and/or  rating  curves 
were  available  for  review  for  about  one  half  of  the  stations  considered  in  the  analysis.  The 
available  information  yielded  the  following  general  observations  which  should  apply  for  all 
stations. 


■  The  rating  curves  tend  to  be  unstable,  and  most  stations  required  a  number 
of  different  rating  curves  over  time  as  channel  conditions  changed.  Substantial 
shifts  in  rating  curves  were  most  apparent  in  rivers  subject  to  aggradation, 
where  vertical  shifts  of  1  meter  or  more  were  observed  over  the  period  of 
record. 

■  Measured  (gaged)  discnarge  points  are  available  only  at  flows  significantly 
less  than  peak  annual  discharges.  Measured  data  were  apparently  used  as  the 
basis  for  shifting  previously  defined  rating  curves  to  match  the  measured 
points  as  well  as  to  define  new  rating  curves  when  the  extent  of  shifting 
exceeded  certain  (unspecified)  limits.  In  general,  the  low-flow  records  are 
considered  to  be  quite  good  for  those  periods  when  discharge  measurements 
were  made  on  a  regular  basis. 

■  Prior  to  1970,  dischaige  measurements  were  made  on  a  regular  basis  at  most 
stations,  typically  at  least  three  measurements  per  year.  After  1970,  the 
frequency  of  discharge  measurements  decreased  significantly,  with  most 
stations  showing  one  or  more  years  with  no  measurements  at  all. 

■  Plotted  stage-discharge  curves  appear  to  be  biased  towards  matching  the 
lower  limit  of  measured  discharges  for  a  given  stage.  If  this  observation  is 
accurate,  the  implication  is  that  available  water  supplies  are  conservatively 
reported  but  that  peak  discharges  tend  to  be  underestimated. 

■  The  basis  for  extrapolating  the  rating  curves  beyond  the  measured  discharges 
is  generally  not  known.  The  following  three  methodologies  were  indicated  in 
the  literature  for  some  of  the  stations: 

1)  "The  upper  portion  of  the  curve  was  extended  by  area-velocity  method." 

2)  "Due  to  the  unavailability  of  the  cross-section  of  the  river,  the  upper 
portion  of  the  curve  was  extended  by  logarithmic  plotting." 

3)  There  are  infrequent  references  to  some  peak  discharges  as  having  been 
determined  by  "slope-area"  method. 

In  summary,  pre-1970  reported  low  discharges  are  considered  to  be  quite  accurate  but  the 
discharges  may  not  reflect  natural  hydrologic  conditions  due  to  irrigation  withdrawals  on 
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many  of  the  rivers.  Peak  flows  are  subject  to  considerable  error  due  to  a  lack  of  data  for 
estimation  plus  shifting  of  the  stage-discharge  relationships  during  flood  events. 

Screening  of  Daily  Data.  Daily  stneamflow  data  were  screened  to  identify 
clearly  anomalous  (high  and  low)  values  and  to  summarize  periods  of  missing  data.  Data 
availability  is  summarized  in  the  time-line  in  Figure  2.4. 1 . 

The  principal  observations  from  screening  the  data  are  as  follows: 

■  Useful  streamflow  data  are  sparse  after  1972.  Data  available  after  1972 
suggest  staff  gages  were  read  not  on  a  regular  daily  basis,  but  at  irregular 
intervals  several  times  a  month.  For  example,  the  record  for  Maloma  for 
1987  contains  a  period  of  21  days  in  the  wet  season  with  a  constant  flow  of 
7.2  m’/s  with  an  abrupt  change  on  August  17  from  7.2  to  372  m’/s.  The 
records  for  1987  and  other  yean  contain  similar  periods  of  "constant"  flow. 

■  As  can  be  seen  from  Figures  2.3.1  and  2.4.1,  there  is  very  limited  overlap 
of  daily  rainfall  and  daily  streamflow  data.  This  makes  it  impossible  to 
reconstruct  storm  isohyetal  maps  for  which  reliable  concurrent  streamflow  data 
are  available. 

■  Data  from  the  two  stations  on  the  Porac  River  near  Del  Carmen  and  Valdez 
show  inconsistencies.  It  is  possible  that  the  data  near  Valdez  reflect  irrigation 
divenions  or  divenion  of  flows  into  the  Gumain  Floodway.  However,  no 
information  is  available  on  the  nature  of  diversions  or  the  physical 
configuration  of  the  floodway.  As  a  result,  data  from  the  Porac  River  near 
Valdez  and  the  Gumain  Floodway  were  dropped  from  further  consideration. 

■  Records  from  many  stations  show  unreasonable  abrupt  changes  in  flow  rate. 
These  occur  at  many  places  throughout  the  record  and  are  too  numerous  to 
document  individually. 

Double  Mass  Analysis.  Double  mass  analysis  was  conducted  on  daily 
streamflow  data  using,  in  the  first  instance,  the  following  group  of  four  stations: 

Porac  River  near  Del  Carmen  (084A) 

Gumain  River  (086A) 

Bucao  River  (093A) 

Santo  Tomas  River  (094A) 

These  stations  were  selected  because  they  are  of  direct  interest  to  the  hydrologic  analyses  (all 
are  affected  by  the  eruption  of  Mount  Pinatubo)  and  have  a  relatively  long  common  record. 
Hie  entire  concurrent  record  for  the  group  of  stations  was  used  for  double  mass  analysis. 

The  analysis  was  done  in  terms  of  mean  monthly  cumulative  'unoff  in  millimeters  using 
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published  drainage  areas  to  convert  from  discharge  rate  in  m’/s  to  runoff  depth  in 
millimeters.  The  plot  for  Porac  River  near  Del  Carmen  indicates  a  progressive  increase  in 
flows  relative  to  other  stations  throughout  its  period  of  record.  The  plot  for  the  Santo  Tomas 
indicates  a  dramatic  reduction  in  flows  starting  in  1964.  The  reason  for  the  reduction  in 
flows  is  not  known  but  could  result  from  gage  errors  or  irrigation  diversions.  However,  no 
information  was  available  as  to  whether  such  diversions  occurred.  Plots  for  the  Gumain  and 
Bucao  Rivers  showed  a  reasonably  consistent  record. 

For  completeness,  double  mass  analyses  were  done  for  all  other  available  records  by  plotting 
cumulative  monthly  runoff  at  the  station  of  interest  against  the  cumulative  mean  monthly 
runoff  from  a  control  group  made  up  of  stations  on  the  Porac,  Gumain,  Bucao,  and  Santo 
Tomas  Rivers.  All  plots  exhibited  inconsistencies  in  the  records.  Records  for  the  two  gages 
on  the  Camiling  River  appeared  to  be  the  most  reliable  outside  the  control  group. 

Review  of  Peak  Annual  Flow  Data.  Table  2.4.2  presents  a  summary  of  all 
available  peak  flow  data  expressed  in  m’/s.  Table  2.4.3  presents  the  same  data  on  a 
normalized  yield  basis  computed  by  dividing  each  peak  flow  (m’/s)  by  the  published  basin 
drainage  area  (km1)  at  the  gage.  These  tables  include  data  found  to  be  unrepresentative  of 
actual  annual  peak  flow  amounts. 

Some  of  the  published  peak  annual  flow  data  were  determined  to  be  unrepresentative  of 
actual  annual  peak  flow  amounts  when  the  peak  was  published  for: 

•  a  year  in  which  there  was  only  a  parual  record  of  flow,  and  for  which  there 
were  no  discharge  records  for  significant  portions  of  the  high-flow  months  of 
July  through  December, 

■  a  year  or  longer  period  for  which  there  were  serious  questions  as  to  the 
adequacy  of  the  rating  curve  used  at  the  station. 

Tables  2.4.4  and  2.4.5  summarize  peak  flow  and  normalized  peak  flow  data  which  have  been 
screened  to  exclude  doubtful  records. 

Examination  of  data  in  Tables  2.4.2  through  2.4.5  show  very  low  peak  flow  yields  for  the 
Pasig-Potrero  and  O’Donnell  Rivers  and  very  high  yields  for  the  Caulaman  River.  The  low 
yields  on  the  Pasig-Potrero  and  O’Donnell  could  result  from  local  geologic  conditions 
although  no  information  is  available  to  confirm  this.  The  reason  for  the  unusually  high  yield 
on  the  Caulaman  is  unknown. 

Most  streamflow  gaging  stations  used  in  the  hydrologic  analyses  rely  on  a  staff  gage  (which 
is  read  between  once  and  three  times  a  day,  depending  on  the  gage  site)  and  its  associated 
rating  curve.  The  reported  "peak  annual  flow"  at  these  stations  appears  to  be  the  largest  of 
the  discrete  number  of  available  observations.  There  are  no  known  crest-stage  gages  at  the 
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gage  sites,  and  few,  if  any,  of  the  repotted  peak  flow  measurements  appear  to  be  based  on 
observed  high  water  marks.  The  following  points  are  noted: 

■  Because  continuous  stage  records  are  lacking  at  many  sites,  many  reported 
“peak  annual  flows"  may  understate  true  instantaneous  peak  annual  flows. 
However,  due  to  the  considerable  uncertainty  in  the  magnitude  of  high  flows 
resulting  from  extrapolation  of  rating  curves,  it  cannot  be  determined  if  the 
reported  "peak  annual  flows"  are  in  fact  understated. 

■  Measurement  procedures  appear  to  be  inconsistent  from  year  to  year.  In 
some  years,  some  stations  report  the  peak  annual  flow  as  having  the  same 
value  as  the  maximum  daily  annual  flow,  implying  that  only  one  stage 
measurement  was  made  on  that  day,  even  though  the  gage  was  reported  to  be 
read  two  or  three  times  a  day. 

Review  of  Annual  Runoff  Data.  The  annual  runoff  (mm)  for  each  complete 
year  of  streamflow  record  was  computed  for  all  stations  and  is  shown  in  Table  2.3.7.  Some 
serious  inconsistencies  are  evident: 

*  The  Porac  River  near  Valdez  (drainage  area  118  km1)  shows  a  significantly 
lower  yield  than  the  Porac  River  near  Del  Carmen  (drainage  area  111  km5). 
The  Valdez  record  is  presumably  affected  by  irrigation  diversions  or  diversions 
to  the  Gumain  Floodway. 

■  The  Santo  Tomas  River  shows  an  abrupt  drop  in  yield  in  1967  (drop  was 
also  indicated  on  the  double  mass  plot). 

■  The  Pasig-Potrero  River  at  Hacienda  Dolores  (drainage  area  28  km5)  shows 
an  abrupt  drop  in  yield  in  1969. 

Because  of  significant  periods  of  missing  data,  the  annual  runoff  data  from  several  stations 
could  only  be  estimated  for  two  or  three  years  and  were  not  considered  to  be  of  value  to  the 
hydrologic  analyses. 

Due  to  lack  of  overlapping  rainfall  and  streamflow  data,  only  limited  comparison  of  annual 
runoff  and  annual  rainfall  data  was  possible.  Plots  of  annual  rainfall  at  Clark  AFB  against 
annual  runoff  on  the  Gumain,  Porac,  and  O’Donnell  Rivers  are  given  in  Figures  2.4.2 
through  2.4.4,  and  a  plot  of  annual  rainfall  at  Cubi  Point  NAS  against  runoff  on  the  Santo 
Tomas  is  given  in  Figure  2.4.5.  With  the  exception  of  the  Gumain  River  (Figure  2.4.2), 
these  plots  do  not  show  a  good  relationship  between  annual  rainfall  and  annual  runoff.  The 
plot  for  the  Santo  Tomas  River  (Figure  2.4.5)  would  improve  somewhat  if  suspect  data  after 
1967  were  dropped  from  the  analysis. 
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Review  of  Gaping  Stations.  The  records  available  for  analyses  are  listed  in 
Table  2.4.6.  The  principal  findings  of  the  review  of  stations  from  which  these  records  were 
obtained  are  provided  below  for  each  gage  site: 

Bulsa  River  (Station  W010A):  The  headwaters  for  the  Bulsa  River  are  on  the 
east  slopes  of  the  Zambales  Mountains,  15  to  45  km  north  of  Mount  Pinatubo. 
The  record  shows  unusually  high  yield  (i.e.  normalized  peak  flow  in  m'/s/km') 
relative  to  other  stations  whose  headwaters  originate  on  Mount  Pinatubo,  with 
some  periods  reporting  extreme  (and  likely  erroneous)  monthly  runoff  (e.g. 
runoff  for  July  1972  was  reported  as  3.8  meters,  with  monthly  rainfall  at  Iba 
of  1.7  meters  and  at  Hacienda  Luisita  1.6  meters).  The  high  yield  relative  to 
basins  on  Mount  Pinatubo  may  result  in  part  from  different  geologic 
conditions.  However,  the  Bulsa’ s  double  mass  curve  is  concave  upward, 
indicating  a  progressive  increase  in  flows  with  time  relative  to  other  stations, 
and  its  rating  curves  show  a  progressive  upward  shift  due  to  aggradation. 

O'Donnell  and  Bangat  Rivers  (Stations  W011A,  W011B,  and  W012A):  The 
records  from  the  O’Donnell  and  Bangat  Rivers  show  very  significant  and 
irreconcilable  inconsistencies.  Peak  flows  on  the  Bangat  (Station  W012A, 
drainage  area  90  km1)  are  invariably  much  higher  (by  as  much  as  a  factor  of  6) 
than  those  recorded  at  the  downstream  gaee  on  the  O’Donnell  (dramage  area 
240  km1).  A  possible  high  flow  diversion  exists  from  the  O’Donnell  above 
gage  W011B  into  the  Bangat  above  gage  W012A.  However,  no  diversion  has 
been  identified  out  of  the  system  between  the  Bangat  River  gage  W012A  and 
the  O’Donnell  River  gage  W011A.  Records  for  these  stations  are  less  than 
10  years  in  length. 

Camiline  River  (Stations  W023A  and  W023B):  The  headwaters  of  the 
Camiling  River  are  on  the  east  slopes  of  the  Zambales  Mountains  from  40  to 
60  km  north  of  Mount  Pinatubo.  The  records  show  an  unusually  high  yield 
(i.e.  normalized  peak  flow  in  m’/s/km1)  relative  to  other  stations  originating  on 
Mount  Pinatubo.  The  record  for  both  these  stations  is  less  than  10  years  in 
length. 

Pasig-Potrero  River  (Stations  W081A  and  W082A):  The  record  from  station 
W081A  (drainage  area  242  km1)  may  have  been  tidaliy  affected  and  stage 
records  only  are  available  for  much  of  the  record.  Discharge  records  are 
available  for  only  five  years  from  Station  W081A  and  six  years  from  Station 
W082A. 
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Porac  River  (Stations  W083A  and  W084A):  The  records  from  the  two  Porac 
River  gages  (Station  W083A  with  drainage  area  118  km'  and  W084A  with 
drainage  area  111  km1)  show  significant  and  irreconcilable  inconsistencies, 
with  the  downstream  gage  frequently  showing  significantly  higher  peak  flows 
than  the  upstream  gage  despite  the  very  small  difference  in  drainage  area.  It  is 
likely  that  the  records  are  affected  by  the  operation  of  both  irrigation  and  flood 
control  projects.  However,  no  information  could  be  obtained  on  the 
configuration  or  operation  of  these  schemes. 

Gumain  Floodwav  (Station  W085A):  This  station  is  at  the  downstream  end  of 
the  Gumain  Floodway  and  the  record  is  affected  by  upstream  flood  control  and 
irrigation  projects.  No  information  could  be  obtained  on  the  configuration  or 
operation  of  these  projects  and  how  they  affect  the  flow  record  at  gage 
W085A. 

Gumain  River  (Station  W086A):  The  headwaters  for  the  Gumain  River  are  on 
the  southeast  slopes  of  Mount  Pinatubo.  The  available  daily  record  is  15  years 
in  length.  There  are  no  known  upstream  diversions  into  or  out  of  the  system. 
The  record  appears  to  be  of  reasonable  quality.  A  water  level  recorder  was  in 
operation  for  most  of  the  period  of  record. 

Caulaman  River  (Station  W087A):  The  Caulaman  River  originates  on  the  east 
slopes  of  Mount  Bitnung  i5  km  south  of  Mount  Pinatubo.  It  has  not  been 
possible  to  determine  the  exact  location  of  this  gage.  Assuming  the  reported 
drainage  area  is  correct,  the  records  show  an  extremely  high  yield  relative  to 
other  stations.  Reported  monthly  and  event  runoff  is  occasionally  (and  likely 
erroneously)  extremely  high.  For  example,  the  reported  runoff  for  September 
1963  was  1.85  meters  with  reported  rainfall  at  Cubi  Point  NAS  of  0.93  meters 
and  at  Clark  AFB  of  0.51  meters. 

Colo  River  (Station  W088A):  The  Colo  River  basin  lies  approximately  30  km 
south  of  Mount  Pinatubo.  The  station  is  downstream  from  an  irrigation  dam. 
The  effect  of  this  dam  on  the  record  is  not  known. 

Ba°s:t  River  (Station  W092A):  The  Bagsit  river  originates  on  the  west  slopes 
of  the  Zambales  Mountains  approximately  40  km  north  of  Mount  Pinatubo. 
The  basin  was  judged  to  be  too  far  from  Pinatubo  to  be  representative  of 
hydrologic  conditions  of  interest. 
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Bucao  River  (Station  W093A):  The  Bucao  River  originates  on  the  northwest 
slopes  of  Mount  Pinatubo.  The  available  record  of  daily  flows  is  15  years  in 
length.  Although  the  record  quality  is  uncertain  because  of  extreme 
extrapolation  of  the  available  rating  curves,  the  record  prior  to  1972  appears  to 
be  relatively  consistent.  There  are  occasional  periods  of  suspiciously  high 
reported  daily  ninoff  volumes.  A  water  level  recorder  was  apparently  in 
operation  for  part  of  the  record. 

Santo  Tomas  (Station  W094A):  The  Santo  Tomas  River  originates  on  the 
southwest  slopes  of  Mount  Pinatubo.  The  available  record  of  daily  flows  is 
15  years  in  length.  The  double  mass  curve  for  the  Santo  Tomas  showed 
discharge  volumes  after  1967  to  be  significantly  reduced  by  irrigation 
diversions  upstream  of  the  gage.  However,  the  11  years  of  data  prior  to  1967 
appears  to  be  relatively  consistent. 

Maloma  River  (Station  W099B):  The  Maloma  River  originates  on  the  west 
slope  of  Mount  Pinatubo.  Seven  years  of  flow  data  are  available;  however, 
the  record  is  judged  to  be  exceedingly  poor  and  too  unreliable  to  be  of  use. 

2.4.2  Stream  flow  Characteristics.  Table  2.4.7  lists  the  15  streamflow  gages 
considered  in  the  hydrologic  analyses,  together  with  rainfall  ard  evaporation  gages.  The  list 
excludes  those  gages  determined  to  be  unreliable  or  not  useful  to  the  analyses.  Gage 
locations  are  shown  on  Plate  1 .  The  data  available  for  each  streamflow  gage  consist  of 
reported  mean  daily  discharges  and  annual  peak  instantaneous  discharges.  No  continuous  or 
short-duration  flow  hydrograph  data  are  available. 

The  convention  adopted  to  identify  streamflow  gage  data  consists  of  a  five-character  code: 

•  the  first  character  is  *W"  to  signify  a  streamflow  gage; 

•  the  second  through  fourth  characters  identify  the  stream  gage  number,  based 
where  possible  on  gage  numbers  published  by  Philippine  agencies; 

•  the  fifth  (and  final)  "A”  or  "B”  character  signifies  whether  the  gage  actually 
had  a  published  gage  number  "A"  signifies  that  the  number  was  published, 
and  "B"  signifies  that  no  gage  number  had  been  published  and  that  one  was 
arbitrarily  assigned  for  purposes  of  this  analysis. 

Periods  of  record  of  available  daily  streamflow  data  are  summarized  by  Table  2.4.7  and 
Figure  2.4.6  together  with  rainfall  and  evaporation  data.  Record  lengths  vary  from  five  to 
16  years  and  are  mostly  within  the  period  1957  to  1972.  Peak  instantaneous  streamflow  data 
considered  in  this  analysis  are  summarized  by  Table  2.4.8.  Record  lengths  vary  from  five  to 
33  years,  with  very  little  data  available  after  1972. 
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Summary  hydrographs  at  each  of  the  gages  are  presented  on  Figures  2.4.7  through  2.4.21. 
The  high  points  shown  on  the  figures  reflect  discrete  major  flood  events;  differences  in  the 
presence  or  absence  of  specific  peaks  from  gage  to  gage  are  due  in  part  to  differences  in  the 
periods  of  record. 

Plots  of  mean  monthly  minimum,  average,  and  maximum  daily  discharges  are  presented  on 
Figures  2.4.22  through  2.4.36.  The  mean  monthly  maximum  daily  discharge  for  the  month 
of  June,  for  example,  was  determined  by  averaging  the  maximum  daily  values  for  June  from 
each  year  of  record. 

The  plots  on  Figures  2.4.7  through  2.4.36  indicate  a  strong  seasonal  pattern  of  low  flows 
during  the  months  of  January  through  April  and  high  flows  during  the  months  of  June 
through  October.  High  flows  indicated  for  May  mostly  reflect  a  major  storm  in  May  1966. 
The  Bucao  River  (W093A)  has  no  data  for  this  month;  hence,  the  Bucao  River  data  as  shown 
on  Figures  2.4.19  and  2.4.34  are  misleading  by  wrongly  suggesting  that  high  flows  have  not 
occurred  on  this  river  in  May. 

Plots  on  Figures  2.4.37  through  2.4.51  sliow  the  daily  flow  duration  curves  for  each  of  the 
streamflow  gages.  Plots  on  Figures  2.4.52  and  2.4.53  show  normalized  flow  duration  curves 
for  all  gages,  in  which  the  curves  were  normalized  by  dividing  the  curve  ordinates  by  the 
average  daily  flow  for  each  gage. 

The  normalized  flow  duration  curves  show  pronounced  differences  between  the  gages  at 
higher  discharges.  While  the  reasons  for  these  differences  are  not  known,  inaccuracies  in  the 
reported  peak  discharges  are  probably  one  factor. 

2.4.3  Runoff.  Figures  2.4.54  through  2.4.68  summarize  the  annua!  runoff  for  each 
of  the  streamflow  gages. 

2.4.4  Frequency  Analyses.  The  HEC-FFA  Flood  Frequency  Analysis  program 
Version  3.0,  which  computes  flood  frequencies  using  a  Log  Pierson  Type  3  fit  in  accordance 
with  guidelines  described  in  Bulletin  17B  of  the  U.S.  Water  Resources  Council1,  was  used  to 
assess  frequency  characteristics  of  both  flood  flow  and  rainfall  data. 

Rainfall  Data.  HEC-FFA  frequency  analyses  of  maximum  annual  1,  2,  5,  10, 
and  15-day  duration  rainfall  amounts  were  conducted  for  all  daily  data  rain  gage  stations. 

Tables  2.4.9  through  2.4.13  summarize  the  rainfall  frequency  data  as  computed  from  the  raw 
data  and  also  after  multiplication  by  factors  to  convert  from  observational  day  data  to  n-hour 


1  U.S.  Water  Resources  Council,  March  1982,  Guidelines  for  Determining  Flood  Flow 
Frequencies.  Bulletin  17B. 
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amounts.  The  following  empirical  factors  suggested  hy  the  U.S.  Weather  Bureau*-7  were 
used: 


n-davs 

n-hours 

conversion  factor 

1 

24 

1.13 

o 

48 

1.04 

5 

120 

1.02 

10 

240 

1.01 

15 

360 

N/A 

The  plotted  analytical  frequency  curves  fit  the  data  well  for  all  n-day  durations  at  all  gages  as 
typified  hy  Figures  2.4.69  and  2.4.70. 

The  frequency  data  support  the  earlier  observation  that  daily  and  multi-Jay  rainfall  is 
significantly  greater  along  the  coast  than  in  the  interior. 

Streamflow  Data.  The  criteria  for  selecting  stations  for  frequency  analysis 
were  generally  as  follows: 

•  minimum  of  10  years  record  (excluding  periods  or_  clearly  erroneous  data); 

•  headwaters  originating  on  Mount  rhsatubo  (records  from  streams  originating 
in  other  areas  were  not  selected  for  frequency-based  calibration  because  of 
lack  of  information  or,  surficial  geology  and  Concern  about  whether  such 
records  would  be  representative  of  hydrologic  conditions  on  Mount  Pinatubo); 

•  clearly  defined  drainage  »jex  with  consistent  records  and  no  significant 
upstream  diversions  into  or  cut  of  the  system  (e.g.,  for  irrigation  or  flood 
control). 

It  was  determined  in  the  review  of  gaging  stations  in  Section  2.4.1  that  ttuee  stations  met  the 
above  criteria.  These  were  the  Gumain  River  (WU8bA),  the  Bacao  2jvet  (W0934,\  and  the 
Santo  Tomas  River  (W094A). 


*  Weather  Bureau,  U.S.  Dept,  of  Commerce,  no  date,  Rai.ful  Ftiquency  Alias  oj  ‘.he  United 
Stales  for  Durations  from  30  Minutes  to  24  Honrs  and  Return  Periods  from  !  to  100  Years. 
Technical  Paper  No.  40. 


7  Weather  Bureau,  U.S.  Dept,  of  Commerce,  1964,  7 wo-  to  Ter.-Phty  Pmnptteils it  for  Return 
Periods  of  2  to  100  Years  in  the  Contiguous  United  States.  Technical  Paper  No.  49. 
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For  each  of  these  gages,  frequency  analyses  were  conducted  on  the  following  maximum 
annual  data: 


A)  Peak  instantaneous  discharge 

B)  Mean  daily  discharge  (and  corresponding  1-day  volume) 

C)  Mean  3-day  discharge  (and  corresponding  3-day  volume) 

Also  tor  each  of  the  gages,  different  data  sets  were  assessed  according  to  data  availability 
and  reliability,  as  follows: 

A)  Screened  data  set  witnin  the  period  1957  through  1972.  As  previously 
discussed,  these  data  were  screened  by  double-mass  analysis  and  other 
methods  to  eliminate  obviously  unreliable  data. 

8)  Extended  data  set  including  the  screened  data  plus  additional  pre-1957  and 
post-19r/’  annual  maximum  peak  instantaneous  and  daily  dischaiges.  The 
extended  portion  of  this  data  set  was  obtained  from  Philippine  summary  sheets 
listing  only  annual  maximum  peak  instantaneous  and  daily  dischaiges.  Data 
outside  tfte  1957  through  1972  period  ate  considered  "not  screened"  and 
therefore  of  unknown  quality.  Table  2.4.7  shows  periods  of  record  available 
at  each  gage  using  tie  extended  data  set. 

Relative  to  the  ,creened  data  set,  the  extended  data  set  added  16  data  points 
for  frequency  analysts  of  peak  and  daily  flows  on  the  Gumain  River  and  nine 
dau  prints  for  frequency  analysis  of  peak  and  daiiy  flows  on  the  Santo  Tomas 
River.  No  additional  data  points  were  gained  for  the  Bucao  River  or  for 
analysis  of  multi-day  (i.e.,  3-day)  flows  and  volumes  at  any  of  the  gages. 

C)  Subsets  of  'he  above  data  sets  (A  and/or  3),  created  by  excluding  years  in 
which  data  anomalies  w«c  noticed  on  or  near  tits  day  of  maximum  annual 
flow.  In  practice,  data  subsets  were  analyzed  for  the  Santo  Tomas  and  Bucao 
Rivers  onlv,  and  in  each  case  excluder;  data  for  year  1962. 


Gmtrain  River  Frequency  Analysis. 

Two  date  sets  were  analyzed  for  the  Gumain  River  (W086A):  1)  screened  data  for  1957 
through  i971  and  2)  extended  data  for  1947  threugh  1979.  The  screened  data  set  is  believed 
to  )x?  of  high  quality  because  the  data  were  col'ected  using  an  automatic  water  level  recorder. 
Post-1971  data,  and  possibly  much  of  the  pre-1957  data,  were  based  on  staff  gage  readings 
only.  Furthermore,  all  post-1972  data  are  subject  to  greater  uncertainty  due  to  a  significant 
deterioration  in  the  gaging  program  after  that  date.  The  most  significant  difference  between 
the  two  data  sets  is  that  the  maximum  peak  instantaneous  discharge  reported  under  the 


extended  data  set,  740  m’/s  in  1976,  is  nearly  double  the  highest  value  reported  under  the 
screened  data  set,  375  mVs  in  1964. 

Figures  2.4.71  through  2.4.76  show  the  results  of  the  frequency  analyses  on  data  for  the 
Gumain  River  (W086A),  plotting  the  maximum  event  data  together  with  expected  probability 
frequency  curves  as  computed  by  the  HEC-FFA  program.  These  figures  are  discussed 
briefly  in  the  paragraphs  which  follow. 

1)  Figure  2.4.71  shows  the  frequency  analyses  results  arranged  to  show'  the 
peak  instantaneous,  1-day,  and  3-day  flows  together  for  the  screened  data  set. 
Extended  data  set  results  (for  which  3-day  data  were  not  available)  are  shown 
by  Figure  2.4.72.  Figure  2.4.71  shows  that  the  computed  frequency  curves  fit 
the  peak,  1-day,  and  3-day  values  from  the  screened  data  set  reasonably  well, 
although  the  lines  are  not  parallel  as  might  be  expected.  For  the  extended  data 
set.  Figure  2.4.72  shows  a  good  fit  for  the  1-day  data  but  only  a  fair  fit  for  the 
peak  instantaneous  data,  due  to  the  extraordinarily  high  peak  instantaneous 
data  point  for  year  1976,  which  has  a  considerable  influence  on  the  shape  of 
the  frequency  curve  and  on  estimates  of  flood  flows. 

2)  Figures  2.4.73  and  2.4.74  show  the  frequency  analyses  results  arranged  to 
compare  the  effects  of  the  alternative  data  sets  for  the  peak  instantaneous  and 
1-day  flows  respectively.  The  figures  show  that  flows  estimated  for  a  given 
return  period  are  sensitive  to  the  data  set  used.  The  frequency  curves 
computed  for  the  extended  data  set  are  steeper  than  those  for  the  screened  data 
set,  and  differences  between  the  curves  are  greatest  at  higher  return  periods. 
For  example,  expected  probability  estimates  of  the  100-year  peak  instantaneous 
flow  vary  from  about  500  m’/s  based  on  the  screened  data  set  to  660  m’/s 
based  on  the  extended  data  set.  Estimates  of  the  100-year  1-day  flow  vary 
from  about  340  m’/s  based  on  the  screened  data  set  to  420  m’/s  based  on  the 
extended  data  set. 

3)  Figures  2.4.75  and  2.4.76  show  the  1-  and  3-day  volume  frequency 
analyses.  These  data  and  curves  are  the  same  as  for  1-  and  3-day  flows,  after 
conversion  to  volume  units. 

Bucao  River  Frequency  Analysis. 

Two  data  sets  were  analyzed  for  the  Bucao  River  (W093A):  1)  screened  data  for  1957-1965 
and  1967-1971  and  2)  a  data  subset  consisting  of  the  same  screened  data  in  (1)  but  excluding 
data  from  year  1962.  Data  from  year  1966  were  excluded  from  both  of  the  data  sets  because 
there  were  no  data  for  that  year’s  peak  flow  month  (based  on  other  gages)  of  May. 

The  year  1962  was  excluded  from  the  second  Bucao  River  data  set  because  a  single  value, 
2,220  m’/s,  was  reported  to  be  both  the  peak  instantaneous  discharge  and  the  maximum 
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average  daily  discharge  for  the  year.  As  a  peak  instantaneous  discharge,  2,220  m’/s  ranks 
second  in  the  station’s  history.  However,  as  a  maximum  average  daily  discharge,  this  value 
would  be  more  than  double  the  second  highest  daily  value  of  992  mVs  reported  in  1961. 
Furthermore,  total  runoff  from  the  basin  for  the  five-day  period  in  1962  of  July  20-24  was 
942  mm,  compared  to  a  total  rainfall  at  Iba  for  the  same  period  of  only  530  mm.  Even 
allowing  for  spatial  variations  in  rainfall  and  increases  in  rainfall  with  elevation,  the  total 
runoff  volume  appears  unreasonably  large  in  comparison  to  Iba  rainfall. 

The  Bucao  River  data  are  of  mixed  quality.  An  automatic  water  level  recorder  was  in  place 
for  1957  through  1964;  post-1964  data  are  based  on  staff  gage  readings  made  twice  daily. 
The  reported  maximum  annual  discharge  for  year  1962  was  measured/computed  by  the 
"slope-area”  method. 

Figures  2.4.77  through  2.4.83  shew  the  results  of  the  frequency  analyses  on  data  for  the 
Bucao  River  (WC93A),  piotting  the  maximum  event  data  together  with  expected  probability 
frequency  curves  as  computed  by  the  HEC-FFA  program.  These  figures  are  discussed 
briefly  in  the  paragraphs  which  follow. 

1)  Figures  2.4.77  and  2.4.78  show  the  frequency  analyses  results  arranged  to 
show  the  peak  instantaneous,  1-day.  and  3-day  flows  together  for  each  of  the 
data  sets,  i.e.  the  screened  data  set  with  and  without  year  1962  data.  The 
curves  generally  fit  the  data  beat  for  the  data  set  widen  excludes  year  1962 
data. 

2)  Figures  2.4.79  through  2.4.81  show  the  frequency  analyses  results  arranged 
to  compare  the  effects  of  the  alternative  data  sets  for  the  peak  instantaneous,  1- 
day,  and  3-day  flows  respectively.  The  figures  show  that  all  estimated  flows, 
?nd  the  1-  and  3-day  flows  in  particular,  are  sensitive  to  tlte  data  set  used,  and 
that  curves  computed  on  the  data  sets  including  year  1962  data  are  steeper  than 
these  for  which  that  year  was  excluded.  Differences  between  the  curves  are 
greatest  at  higher  return  periods.  For  example,  expected  probability  estimates 
of  the  105- year  peak  ir.sntntaneous.  flow  vaty  from  about  3,740  tn’/s.  excluding 
year  1962  data,  to  4,330  m*/s.  including  1962  data.  Estimates  of  the  ICO-yeat 
1-day  flow  vary  from  about  1,260  mVs  excluding  1962  data  to  2.730  m’/s 
when  including  1962  data. 

3)  Figures  2.4.82  and  2.4.83  show  the  1-  and  3-day  volume  frequency 
analyses.  These  data  and  curves  are  the  same  as  for  1  -  and  3-day  flows,  after 
conversion  to  volume  unit;. 

Santo  Tomas  River  Frequency  Analysis . 
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Four  data  sets  were  analyzed  for  the  Santo  Tomas  Rive  (W094A):  1)  screened  data  for 
1957  through  1967;  2)  extended  data  for  1948  through  1967;  3)  screened  data  set 
(1)  excluding  year  1962  data;  and  4)  extended  data  set  (2)  excluding  year  1962  data. 


All  Santo  Tomas  River  data  are  based  on  staff  gage  readings.  Year  1962,  which  includes  the 
historical  maximum  peak  instantaneous  and  daily  values,  was  excluded  from  the  final  sets  of 
data  because  of  anomalies  in  the  daily  data:  the  minimum  daily  discharge  for  1962  is 
reported  to  occur  just  two  days  after  the  maximum  flow.  As  the  Santo  Tomas  River  gage 
was  located  400  meters  downstream  of  a  structure  that  is  believed  to  be  a  major  irrigation 
project,  it  is  possible  that  high  flows  during  the  flood  may  have  been  supplemented  by 
releases  of  reservoir  storage,  and  that  the  very  low  flows  after  the  flood  reflect  the  recovery 
of  reservoir  storage. 

Figures  2.4.84  through  2.4.92  show  the  results  of  frequency  analyses  on  data  for  the  Santo 
Tomas  River  (W094A),  plotting  the  maximum  event  data  together  with  expected  probability 
frequency  curves  as  computed  by  the  HEC-FFA  program.  These  figures  are  discussed 
briefly  in  the  paragraphs  which  follow. 

1)  Figures  2.4.84  and  2.4.85  show  the  frequency  analyses  rt  ults  arranged  to 
show  the  peak  instantaneous,  1-day  and  3-day  flows  together  for  the  screened 
data  set  with  and  without  year  1962  data.  Extended  data  set  results  (for  which 
3-day  data  were  not  available)  are  shown  by  Figures  2.4.86  and  2.4.87.  The 
curves  are  generally  more  parallel  for  the  extended  data  sets  than  for  the 
screened  data  sets.  The  best  fit  of  expected  probability  curve  to  the  data  is  for 
1-day  discharges  with  the  extended  data  set  excluding  1962  (Figure  2.4.87). 

2)  Figures  2.4.88  through  2.4.90  show  the  frequency  analyses  results  arranged 
to  compare  the  effects  of  the  alternative  data  sets  for  the  peak  instantaneous,  1- 
day,  and  3-day  flows  respectively.  The  figures  show  that  all  estimated  flows 
are  sensitive  to  the  data  set  used,  and  that  curves  computed  on  the  data  sets 
including  year  1962  data  are  steeper  than  those  in  which  that  year  was 
excluded.  Differences  between  the  curves  are  greatest  at  higher  return 
periods.  For  example,  expected  probability  estimates  of  the  100-year  peak 
instantaneous  flow  vary  from  about  800  m’/s  based  on  the  extended  data  set 
excluding  year  1962  data  to  1,560  m’/s  based  on  the  screened  data  set 
including  year  1962  data.  Estimates  of  the  100-year  1-day  flow  vary  from 
about  520  mVs  based  on  the  extended  data  set  excluding  year  1962  data  to 
1160  m’/s  based  on  the  screened  data  set  including  year  1962  data. 

3)  Figures  2.4.91  and  2.4.92  show  the  1-  and  3-day  volume  frequency 
analyses.  These  data  and  curves  are  the  same  as  for  1-  and  3-day  flows,  after 
conversion  to  volume  units. 


Summary. 
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Frequency  analyses  of  the  maximum  strearr.flow  data  are  greatly  complicated  by  the 
uncertain  quality  of  the  available  data  and  related  uncertainty  over  which  of  the  possible  data 
sets  is  more  representative  of  "true"  conditions.  The  choice  of  data  set  significantly 
influences  estimated  high  return  period  flows  and  volumes. 

For  all  three  gages  assessed,  ihe  fit  of  the  computed  curves  to  the  data  was  generally 
improved  by  disregarding  high  peak  values  over  which  there  was  some  justifiable 
uncertainty.  However,  a  review  of  rainfall  data  showed  that  these  and  other  (uncertain) 
extreme  flows  were  associated  with  very  heavy  rainfall  expected  to  result  in  extreme  flows. 

It  is  inappropriate  to  either  fully  accept  or  disregard  the  extreme  value  data. 

Given  the  available  record  of  uncertain  data  and  the  sensitivity  of  the  analyses  to  single  data 
points,  single-value  estimates  of  100-year  or  other  return  period  flows  are  not  made. 

Instead,  the  families  of  expected  probability  frequency  curves  shown  by  Figures  2.4.71 
through  2.4.92  are  used  to  indicate  a  plausible  range  of  reasonable  values.  Ranges  for  2- 
and  100-year  return  periods  are  summarized  in  Table  2.4. 14.  The  results  of  this  analysis 
were  used  to  calibrate  the  hydrologic  model,  HEC-1,  which  was  used  to  generate  design 
flood  hydrographs. 

2.4.5  Design  Storms.  The  objective  of  the  design  siorm  analysis  was  to  develop 
rainfall  byetographs  for  2-,  10-,  50-,  100-,  and  500-year  hypothetical  storms  of  appropriate 
duration  to  result  in  similar  return  period  flooding  of  streams  affected  by  the  Mount  Pinatubo 
eruption.  The  assessment  was  complicated  by  the  fact  that  no  rainfall  data  are  available  to 
describe  conditions  in  the  mountain  watershed  areas  of  interest. 

The  approach  taken  to  develop  the  desigu  storms  is  described  in  the  following  report 
sections.  In  summary: 

1)  Isopluvial  maps  cf  2-,  10-,  50-,  100-,  and  500-year  rainfall  over  durations 
of  1-,  2-,  and  5-days  were  developed  from  the  frequency  analysis  of  available 
nun  data  and  with  the  assumption  that  the  ratios  of  coastal-to-mountain  rainfall 
amounts  would  be  the  same  as  for  annual  rainfall  amounts  shown  by  the  NAP 
map. 

2)  U.S.  Weather  Bureau  area-reduction  factors  were  reviewed  in  light  of  the 
available  rainfall  data,  and  a  methodology  for  applying  area-reduction  factors 
jointly  with  elevation-based  variations  in  rainfall  was  developed. 

3)  Short-duration  (less  than  24-hour)  storm  characteristics  were  assessed  from 
available  published  intensity-duration-frequency  data  and  from  the  limited 
available  hourly  data.  A  methodology  was  oeveloped  to  construct  and 
distribute  design  storm  byetographs  throughout  the  basins  of  interest. 
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Freouencv-Du  ration  Isopluvial  Maos.  Figures  2.4.93  through  2.4.107  show 
isopluvial  maps  of  2-,  10-,  50-,  100-,  and  500-year  rainfall  over  durations  of  1-,  2-,  and 
5-days. 

The  isopluvial  lines  shown  for  the  coast  and  interior  lowlands  are  based  on  the  rainfall 
frequency  analyses  results  summarized  by  Tables  2.4.9  through  2.4.13.  Isopluvial  lines  for 
the  mountain  watershed  areas  are  based  on  the  assumption  that  rainfall  near  the  1 ,000-meter 
contour  around  the  summit  of  Mount  Pinatubo  is,  for  all  return  periods,  equal  to  the  average 
rainfall  at  the  coastal  gages  multiplied  by  a  factor  of  1.35  for  all  durations  and  return 
periods.  The  1, 000-mcter  contour  follows  the  top  of  a  narrow  ridge  to  the  south  of  Mount 
Pinatubo  and  an  irregular  circular  patn  about  5  km  in  diameter  around  Pinatubo’s  2.5  km 
diameter  (post-eruption)  crater  rim.  Tne  1.35  factor  is  based  on  average  annual  data  for 
coastal  rain  gages  and  on  the.  NAP  map. 

Tables  2.4.15  through  2.4.17  provide  data  to  examine  the  assumption  that  the  ratio  of  same- 
retum-period  rainfall  at  the  coast,  at  high  elevations,  and  in  the  interior  is  relatively  constant 
for  all  durations  and  return  periods.  Table  2.4.15  presents  data  and  ratios  of  interior  gages 
to  coastal  gages  for  multi-day  duration  events.  Tables  2.4.16  and  2.4.17  present  short- 
duration  data  and  ratios  for  stations  representing  interior,  coastal,  and  high-elevation  sites. 
Figure  2.4.108  plots  the  ratios  for  the  short-duration  data;  Figure  2.4.109  plots  the  ratios 
between  Hacienda  Luisita  and  Iba  for  both  short-  and  long-duration  events. 

Several  observations  are  made  from  the  data: 

1)  For  durations  of  six  hours  and  longer,  the  ratio  of  same-frequency-duration 
rainfall  between  any  two  stations  is  essentially  constant,  independent  of  the 
duration.  This  ratio  is  approximately  equal  to  the  ratio  of  average  annual 
rainfall  between  any  two  stations. 

2)  For  durations  of  one  hour  or  less,  rainfall  depth  is  essentially  independent 
of  station  location  or  elevation.  For  example,  the  100-year  1-hour  rainfall 
depth  at  coastal  sites  is  essentially  the  same  as  the  depth  at  sites  in  the 
mountains  and  the  interior. 

3)  There  is  no  clear  indication  of  the  effect  of  return  period  on  the  ratio  of 
same-frequency-duration  rainfall  between  two  stations.  The  short-term  data 
suggest  that  differences  between  two  stations  become  slightly  greater  (ratios 
further  from  1.0)  with  increasing  return  period.  However,  the  long-term  data 
suggest  the  opposite,  that  differences  become  slightly  less  (ratios  closer  to  1.0) 
with  increasing  return  period.  The  discrepancy  may  result  from  the  fact  that 
the  data  analyses  were  conducted  using  different  theoretical  distributions.  The 
short-duration  data  were  assessed  by  PAGASA  using  a  Gumbell  distribution, 
which  has  a  fixed  skew.  The  long-duration  data  were  assessed  as  part  of  this 
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study  using  a  Log  Pierson  HI  distribution  (HEC-FFA)  in  which  skew  is 

calculated  from  the  data. 

In  summary,  for  rainfall  durations  of  six  hours  and  greater,  the  data  show  that  mountain 
rainfall  of  a  given  frequency  can  reasonably  be  estimated  as  a  constant  multiple  of  same- 
frequency  average  coastal  rainfall.  A  constant  multiplier  equal  to  the  ratio  of  average  annual 
rainfall  appears  reasonable  for  durations  of  six  hours  or  longer.  No  adjustment,  i.e.  a 
multiplier  of  1.0  applied  to  coastal  rainfall,  is  appropriate  for  durations  of  one  hour  or  iess. 

Area  Reduction  Factors.  Area-reduction  factors  are  used  to  estimate  average 
basin-wide  or  area  rainfall  amount  from  point  values.  In  relatively  flat  topography,  the  100- 
year  1-hour  rainfall  over  a  basin  area  is  expected  to  be  less  than  the  100-year  1-hour  rainfall 
at  any  point  within  the  basin. 

Area-reduction  factors  for  the  continental  U.S.  have  been  developed  by  the  U.S.  Weather 
Bureau’  based  on  data  from  dense  rain-gage  networks.  These  factors  have  previously  been 
considered  "reasonable"  to  be  applied  to  conditions  in  Hawaii*,  despite  a  lack  of  data  to  test 
the  relationships  for  Hawaiian  conditions.  After  an  inconclusive  check  for  reasonableness 
relative  to  the  available  study  area  data,  the  Weather  Bureau  factors  were  adopted  for  use  in 
the  hydrologic  analysis. 

Figure  2.4.110  shows  distance-reduction  curves  de  :'  ed  from  the  Weather  Bureau  area- 
reduction  factors.  The  latter,  which  were  based  on  circular  areas,  were  converted  to 
distance-reduction  factors  by  assuming  that  peak  rainfall  occurs  at  the  center  of  the  circle  and 
that  rainfall  decreases  linearly  with  radial  distance  from  the  center.  This  conversion  was 
done  to  facilitate  comparison  with  the  available  data  and  to  provide  a  methodology  which 
could  be  applied  jointly  with  factors  to  account  for  elevation-based  variations  in  rainfall. 

Construction  of  Storm  Hvetographs.  Table  2.4.18  summarizes  frequency- 
depth-duration  data  for  a  hypothetical  rain  gage  located  near  the  summit  of  Mount  Pinatubo 
at  about  1000  meters  elevation.  These  data  provide  the  basis  for  construction  of  design 
storms  hyetographs. 

The  24-hour,  2-  and  5-day  data  in  Table  2.4.18  are  approximately  equal  to  1.35  times 
average  coastal  gage  values  shown  in  Table  2.4.15;  the  data  vary  slightly  from  the  1.35 
multiplier  due  to  rounding  effects  when  tabulating  hourly  data  for  HEC-1  model  input.  The 
1-,  2-,  3-,  6-,  and  12-hour  data  in  Table  2.4.18  are  based  on  the  short-duration  data  of 
Table  2.4.16  with  emphasis  given  to  the  high-elevation  station  at  Baguio  City. 


*  See  Footnote  7. 

’Weather  Bureau,  U.S.  Department  of  Commerce,  no  date.  Two-  to  Ten-Day  Rainfall for  Return 
Periods  of  2  to  100  Years  in  the  Hawaiian  Islands.  Technical  Paper  No.  51. 
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While  peak  discharges  on  rivers  draining  Mount  Pinatubo  generally  result  from  heavy  rains 
lasting  less  than  12  hours,  a  design  storm  duration  of  five  days  was  selected  to  provide 
additional  data  on  total  runoff  volumes  as  might  be  required  for  later  design  work.  Design 
storms  of  a  specified  return  period  were  constructed  by  embedding  same-frequency  events; 
i.e.,  the  100-year  24-hour  rain  was  embedded  within  the  100-year  2-day  rain,  which  was 
embedded  within  the  100-year  5-day  rain.  This  general  approach  was  consistent  with  the 
available  daily  rainfall  data. 

During  early  efforts  to  calibrate  the  HEC-1  model  to  the  region,  it  became  apparent  that  the 
conservative  approach  of  fully  embedding  all  same-frequency  durations  within  a  single  design 
storm  yielded  excessively  high  peak  discharges.  Unfortunately,  peak  discharges  on  the 
basins  of  interest  result  from  short-duration  (typically  3-hour  to  12-hour)  events,  while  hourly 
data  were  available  for  only  two  significant  storm  events. 

Figure  2.4. 1 1 1  plots  hourly  data  and  associated  frequency-duration  characteristics  for  a 
severe  storm  recorded  at  Iba  (RD324)  in  May  1976.  While  the  maximum  48-hour  rain 
corresponds  to  a  200-year  storm  and  the  maximum  24-hour  rain  corresponds  to  a  25-year 
storm,  the  maximum  hour  corresponds  to  only  a  2-year  storm. 

Figure  2.4.112  plots  hourly  data  and  associated  frequency-duration  characteristics  for  a 
severe  storm  recorded  at  Baguio  City  in  September  1911.  The  rainfall  data  were  obtained 
from  Weather  Bureau  Technical  Report  No.  42“  and  reflect  a  (former)  world  record  24- 
hour  rainfall  amount.  While  the  maximum  24-hour  rain  corresponds  to  about  a  200-year 
storm,  the  maximum  hour  corresponds  to  only  about  a  20-year  storm. 

The  available  historic  storm  data  support  the  proposition  that  a  single  design  storm  should 
not  fully  embed  all  same-frequency  depth -durations,  i.e.,  that  a  1  in  100  year  24-hour  rainfall 
should  not  contain  a  1  in  100  year  1-hour  rainfall.  Possible  physical  explanations  for  this  are 
that  the  monsoon  or  typhoon  conditions  responsible  for  large  daily  rainfalls  are  not  conducive 
to  the  formation  of  intense  cloudbursts,  or  that  the  monsoon  or  typhoon  winds  associated 
with  large  daily  rainfalls  do  not  allow  intense  cloudbursts  to  stay  in  one  place  long  enough 
for  significant  rainfall  amounts  to  accumulate  at  a  point. 

Figure  2.4.113  plots  the  2-,  10-,  50-,  100-,  and  500-year  hypothetical  storms  adopted  as 
"base  storms"  before  adjustments  for  elevation-based  total  storm  depth  and  fcr  depth-area 
reduction.  These  hyetographs  are  for  a  hypothetical  station  located  at  approximately  1,000 
meters  elevation  near  the  summit  of  Mount  Pinatubo,  and  embed  same-frequency  events  from 
six  hours  duration  through  five  days  duration.  A  constant  hourly  intensity  is  assumed 
throughout  the  maximum  6-hour  period. 


“  See  Footnote  2. 
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A  constant  hourly  intensity  over  the  maximum  6-hour  period  produces  1-hour  to  24-hour 
characteristics  which  are  similar  to  those  observed  in  the  two  historic  storms  for  which 
hourly  data  were  available.  The  maximum  hourly  rain  in  a  25-year  24-hour  storm  has  a 
return  period  of  about  two  years,  and  the  maximum  hourly  rain  in  a  200-year  24-hour  storm 
has  a  return  period  of  about  20  years. 

The  procedure  described  below  was  developed  to  construct  multi-basin  rainfall  hyetographs 
for  a  single  storm  event.  In  this  description,  "summit  rain"  refers  to  rainfall  characteristics 
at  the  hypothetical  station  at  approximately  1 ,000  meters  elevation  near  the  summit  of  Mount 
Pinatubo.  "Base  storms”  refer  to  design  rainfall  hyetographs  for  this  hypothetical  station  as 
shown  by  Figure  2.4.113. 

1)  The  basin  of  interest  was  subdivided  into  smaller  sub-basins,  considering 
isohyetal  gradients  through  the  basin  as  one  factor. 

2)  The  average  annual  rainfall  for  each  sub-basin  was  determined  from  the 
NAP  map  (Figure  2.3.4),  and  the  ratio  of  sub-basin  rain  to  summit  rain  was 
determined.  For  any  sub-basin,  this  ratio  is  a  constant  for  all  return  periods 
and  for  all  durations  greater  than  six  hours. 

3)  For  each  sub-basin,  the  base  storm  was  adjusted  to  correct  for  the  ratio  of 
sub-basin  rain  to  summit  rain.  The  ratio  was  applied  directly  for  all  durations 
of  six  hours  and  longer,  and  factored  within  the  peak  6-hour  period  so  that  a 
single  peak  hour  was  not  adjusted.  The  justification  for  maintaining  the  peak 
hourly  intensity  at  the  average  6-hour  intensity  for  the  base  storm  is  the 
previously  stated  observation  that  peak  hourly  rainfall  intensities  appear  to  be 
independent  of  location  and  elevation.  Figure  2.4.114  illustrates  this  step  in 
sub-basin  hyetograph  development  by  showing  a  50-year  base  storm  before  and 
after  the  adjustments  for  a  subbasin  normally  receiving  80  percent  of  the 
summit  rain. 

4)  A  storm  center  was  assumed  to  be  located  in  the  middle  of  the  sub-basin 
with  the  highest  rainfall.  The  distances  from  the  assumed  storm  center  to  the 
mid-points  of  all  other  sub-basins  were  determined.  For  the  sub-basin 
containing  the  assumed  storm  center,  the  average  distance  to  the  storm  center 
was  estimated  as  2/3  of  the  radius  of  a  circle  having  an  area  equal  to  the  sub¬ 
basin  area. 
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5)  Finally,  depth-area  corrections  were  applied  to  each  sub-basin  storm  (i.e., 
the  base  storm  after  correction  for  the  ratio  of  sub-basin  rain  to  summit  rain). 
These  corrections  followed  the  depth-duration-distance  curves  shown  by 
Figure  2.4.110  for  durations  of  six  hours  and  greater.  The  6-hour  adjustment 
curve  of  Figure  2.4.110  was  applied  to  all  durations  within  the  peak  6-hour 
period  because  of  the  earlier  assumption  of  a  constant  hourly  intensity  within 
the  peak  6-hour  period  for  the  base  storm.  Figure  2.4.114  illustrates  this  final 
step  in  sub-basin  hyetograph  development  by  showing  a  50-year  sub-basin 
storm  before  and  after  depth-area  corrections  for  a  sub-basin  located  10  km 
from  an  assumed  storm  center. 

Sub-basin  storm  hyetographs  derived  in  this  manner  were  used  as  input  to  HEC-1  models 

used  in  the  basin  analyses. 


3.  BASIN  ANALYSES 

3.1  Introduction 

Hydrologic  basin  analyses  were  performed  on  all  eight  major  river  basins  impacted  by  the 
eruption  of  Mount  Pinatubo.  Hydrologic  analyses  products  were  required  at  specific 
locations  within  the  basins  for  use  in  the  design  of  mitigation  measures.  The  primary 
products  required  from  the  hydrologic  analyses  are,  for  the  purposes  of  the  foliowing 
discussion,  grouped  as  follows:  (1)  2-.  10-,  50-,  100-,  and  500-year  instantaneous  flood 
peaks,  1-day  volumes,  and  3-day  volumes  and  the  hydrographs  corresponding  to  these  peaks 
and  volumes,  (2)  flow  duration  curves,  and  (3)  flow  velocities  and  depths.  The  modeling 
methodologies  used  to  obtain  products  in  groups  (1),  (2),  and  (3)  are  discussed  in  Sections 
3.2,  3.3,  and  3.4  respectively.  The  results  of  the  modeling  are  presented  in  Sections  3.5  and 
3.6. 

3.2  Rainfall/Runoff  Modeling  Methodology 

3.2.1  Introduction.  Rainfall/runoff  processes  were  simulated  with  the  use  of  an 
HEC-1  computer  model.  Model  runs  were  made  with  the  2-,  10-,  50-,  100-,  and  500-year 
hypothetical  storm  events  obtained  from  the  regional  analyses  in  order  to  obtain  estimates  of 
the  2-,  10-,  50-,  100-,  and  500-year  instantaneous  flood  peaks,  1-day  volumes,  and  3-day 
volumes  and  the  corresponding  hydrographs. 

Methods  used  to  obtain  HEC-1  model  parameters  are  discussed  in  Section  3.2.2.  The 
methodology  used  in  the  calibration  of  the  HEC-1  model  is  discussed  in  Section  3.2.3.  The 
methodologies  used  in  the  construction  of  the  basin  models  are  described  in  Section  3.2.4. 

3.2.2  HBC-1  Rainfall/Runoff  Model  Parameters.  Hydrologic  modeling  for  this 
study  was  done  using  an  extended  memory  version  of  HEC-1.  The  extended  memory 
version  of  HEC-1  allows  event  simulation  for  up  to  2,000  computational  time  steps,  as 
compared  to  only  300  time  steps  in  earlier  releases  of  the  model.  The  extended  memory 
version,  however,  assumes  the  time  base  of  the  unit  hydrogtaph  is  less  than  300  time  steps. 
All  modeling  work  was  done  using  metric  units.  Minor  modifications  to  the  extended 
memory  program  were  made  by  the  Hydrologic  Engineering  Center  (HEC)  to  allow  output  of 
unit  hydrogtaph  ordinates  to  three  significant  digits  when  working  in  metric  units. 

Computational  Time  Step.  Hydrologic  modeling  for  rivers  draining  Mount 
Pinatubo  is  complicated  both  by  the  lack  of  good  hydrometeorologic  data  and  by  the  flat 
response  of  event  hydrographs.  Storm  event  hydrog  pits  are  often  quite  flat,  with  recorded 
ratios  of  instantaneous  peak  to  maximum  daily  aver  e  flows  of  as  little  as  1 . 1  or  1 .2. 
Furthermore,  major  storms  affecting  this  area  are  f.  -quently  several  days  in  duration, 
resulting  in  a  need  to  model  hydrologic  response  from  a  storm  event  for  periods  of  up  to 
about  10  days. 
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The  fla<  sponse  of  the  drainage  basins  under  study  naturally  produces  a  flat  unit  hydrograph 
with  an  extended  time  base.  HEC-1  assumes  that  the  length  of  the  unit  graph  time  base  is 
less  than  300  time  steps.  Any  volume  in  a  unit  hydrograph  beyond  this  point  is  lost,  such 
that  the  program  fails  to  conserve  mass.  In  order  to  model  the  flat  response  of  these  rivers 
while  avoiding  unacceptable  loss  of  volume  in  the  tails  of  the  unit  hydrographs,  modeling  had 
to  be  done  at  a  1-hour  time  step. 

While  a  1-hour  time  step  is  appropriate  for  modeling  most  of  the  larger  study  basins, 
modeling  was  also  required  to  simulate  flows  from  a  number  of  small  headwater  catchments 
with  drainage  areas  as  small  as  4.4  km’,  and  with  times  of  concentration  as  low  as  0.4  hour. 
While  this  situation  would  normally  require  modeling  at  a  time  step  of  perhaps  as  small  as  10 
minutes,  tests  conducted  during  the  analysis  found  the  model  response  to  be  relatively 
insensitive  to  time  step  due  to  the  flat  response  of  hydrographs  and  the  relatively  large 
amounts  of  storage  implied.  Tests  conducted  on  small  sub-basins  with  modeling  time  steps 
from  five  minutes  to  one  hour  showed  a  difference  in  simulated  peak  flows  of  about  2 
percent  —  substantially  less  than  the  degree  of  uncertainty  in  the  basic  data  used  for  model 
development. 


Unit  Hydrogranh.  Three  synthetic  unit  hydrographs  are  available  in 

HEC-1: 


•  Clark  Unit  Hydrograph 

•  Snyder  Unit  Hydrograph 

•  SCS  Dimensionless  Unit  Hydrograph 

The  shape  of  the  SCS  unit  hydrograph  is  controlled  by  a  single  parameter,  namely,  basin 
time  lag.  Preliminary  simulations  showed  it  to  be  incapable  of  reproducing  the  peak-to- 
volurne  characteristics  of  observed  hydrographs  in  the  Pinatubo  region.  The  SCS  unit  graph 
produced  a  much  sharper  response  than  observed,  such  that  it  was  not  possible  to  maintain 
both  observed  peaks  and  volumes.  Most  of  the  HEC-1  modeling  was  therefore  done  using 
the  Clark  Unit  hydrograph.  This  method  uses  three  parameters  to  compute  a  unit  graph: 
time  of  concentration,  a  storage  coefficient,  and  a  time-area  curve.  The  greater  flexibility 
inherent  in  the  Clark  method  was  found  to  greatly  improve  modeling  results.  The  Clark 
storage  coefficient  provides  an  avenue  for  modeling  significant  attenuation  of  peak  flows  in 
the  highly  permeable  volcanic  deposits  on  Mount  Pinatubo.  The  SCS  unit  graph  was  retained 
and  used  to  represent  runoff  from  the  relatively  large  impervious  paved  areas  of  two  sub¬ 
basins  in  the  Abacan  River  catchment  that  incorporate  parts  of  Clark  AFB  and  nearby  urban 
areas. 

Parameters  for  the  Clark  unit  hydrograph  were  determined  as  follows: 
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•  Time  of  Concentration 


In  the  absence  of  observed  data,  time  of  concentration  was  determined  by  an  empirical 
equation  provided  by  the  U.S.  Bureau  of  Reclamation"  as: 


where 


T,  =  time  of  concentration  (hours) 

L  =  length  of  longest  watercourse  from  the  point  of 
interest  to  the  watershed  divide  (miles) 
tJH  =  deration  change  along  the  longest  water  course 
(feet) 


•  Storage  Coefficient 

The  Clark  storage  coefficient  was  determined  by  calibration  to  available  data,  as  will  be 
discussed  in  detail  in  Section  3.2.3,  and  expressed  as  a  function  of  7,  to  facilitate  application 
to  ungaged  sub-basins. 

•  Time-Area  Curve 

The  time-area  curve  describes  the  cumulative  area  of  a  sub-basin  contributing  runoff  to  the 
sub-basin  outlet  as  a  function  of  time,  expressed  as  a  fraction  of  the  time  of  concentration. 
The  default  time-area  curve  provided  by  HEC-1  was  used  for  all  sub-basins.  HEC-1  uses  a 
dimensionless  time-area  curve  given  by  the  following: 


AI  =  1.414PJ 

0^T<0.5 

1  -  AJ  =  1.414(1  -  T)1-’ 

0.5£T<1 

where 

AI  =  cumulative  area  as  a  fraction  of  the  total  sub-basin 

area 

T  =  fraction  of  the  time  of  concentration 

Parameters  for  the  SCS  unit  hydrograph  were  determined  as  follows: 


11  United  States  Department  of  the  Interior,  1977,  Design  of  Small  Dams,  Bureau  of 
Reclamation,  Washington  D.C. 
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•  Time  Lag 

The  time  lag  is  the  lag  (in  hours)  between  the  center  of  mass  of  rainfall  excess  and  the  peak 
of  the  unit  hvdrograph.  The  lag-time  v'as  determined  using  an  empirical  expression 
recommended  by  the  U.S.  Soil  Conservation  Service: 


T  = 

0.6 -T, 

where 

tl  = 

lag  time,  hours 

T.  = 

time  of  concentration  in  hours,  estima'ed  using  the 
same  expression  as  described  above  for  the  Clark 
unit  hydrograph. 

Base  Flows.  Base  flows  in  the  context  of  the  basin  analyses  are  representative 
flood  season  low  flows  expected  to  occur  at  simulation  sites  prior  to  the  onset  of  the  flood 
hydrograph.  They  are  significantly  higher  than  the  low  base  flows  which  occur  during  the 
dry  season. 

Flood-season  base  flows  at  stream  gages  were  estimated  by  inspection  of  time  series 
hydrograph  plots.  Average  annual  rainfall  for  each  basin  was  estimated  (as  before)  to  be  the 
average  annual  runoff  plus  1,550  mm  to  account  for  evapotranspiration.  Comparison  of 
these  data  yielded  the  following  approximate  relationship  for  estimation  of  flood-season  base 
flows  on  streams  originating  on  Mount  Pinatubo  as  a  function  of  average  annual  basin 
rainfall  and  basin  area: 


Qbase  =  A  (R/17,500  -  0.12) 

where  Qbase  =  Typical  flood  season  base  flow,  m’/s 

A  =  Basin  area,  knf 

R  =  Average  annual  rainfall  over  basin,  mm/yr 

Base  flows  determined  by  this  relationship  for  each  sub-basin  were  assumed  to  be  constant 
(without  any  recession)  throughout  each  event  simulated  in  the  HEC-1  modeling. 

Loss  Rates.  The  HEC-1  model  allows  the  computation  of  loss  rates  by  a 
number  of  different  methods.  These  methods  include: 

•  initial  and  uniform  loss  rate 

•  exponential  loss  rate 

•  SCS  Curve  Number  method 

•  Holton  loss  rate 
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All  modeling  for  this  study  assumed  a  uniform  loss  rate  with  no  initial  loss.  Major  floods  in 
this  area  occur  primarily  during  the  Southwest  Monsoon.  Significant  amounts  of  rainfall  can 
be  expected  prior  to  flood-producing  events  which  are  therefore  assumed  to  occur  under 
essentially  saturated  conditions. 

The  uniform  loss  rate  was  estimated  by  calibration  to  available  data  as  will  be  described  in 
Section  3.2.3. 


Routing  Parameters.  Routing  of  flood  hydrographs  through  the  stream 
channels  in  the  basins  of  interest  was  done  in  HEC-1  using  the  Muskingum  method.  There 
is  a  lack  of  basic  data  from  which  to  determine  channel  routing  parameters.  Only  limited 
information  exists  on  channel  cross-sections  or  floodplain  geometry,  no  information  exists  on 
flood  wave  travel  times,  and,  as  discussed  earlier,  there  is  considerable  uncertainty  in 
discharge  rates.  In  the  absence  of  basic  data,  all  flow  routing  was  done  assuming  a  flood 
wave  velocity  of  2.5  m/s  and  a  Muskingum  "X”  coefficient  of  0.2.  A  detailed  description  of 
the  routing  scheme  implemented  in  HEC-1  is  provided  in  the  HEC-1  User’s  Manual. 

3-2.3  Calibration  of  HEC-1  Model.  HEC-1  model  calibration  for  basins  draining 
Mount  Pinatubo  was  limited  by  the  lack  of  reliable  streamflow  data.  From  the  review  of 
gaging  stations  in  Section  2.4.1,  it  was  determined  that  only  three  stations  had  data  reliable 
enough  for  model  calibration.  These  stations  were  the  Gumain  River  (W086A),  the  Bucao 
River  (W093A),  and  the  Santo  Tomas  River  (W094A).  Model  calibration  was  further 
limited  by  the  lack  of  short-interval  rainfall  data  that  were  concurrent  with  data  from  the 
above  three  stations. 

Figures  3.2.1  through  3.2.3  show  the  calibration  basins  as  they  existed  prior  to  the  eruption 
of  Mount  Pinatubo.  Sub-basins  upstream  of  the  stream  gages  were  delineated  to  reflea 
major  tributaries  and  to  reflea  the  elevation-related  rainfall  gradients  up  the  mountain  slopes. 

Sub-basins  and  simulation  output  sites  defined  for  pre-eruption  mod^Lng  of  each  of  the  three 
calibration  basins  were  assigned  a  two-letter  prefix  beginning  with  “P"  (for  pre-emption)  for 
identification  purposes.  The  two-letter  prefix  for  each  of  the  three  calibration  basins  is: 

PG  —  Gumain  River  (W086A) 

PB  —  Bucao  River  (W093A) 

PT  —  Santo  Tomas  River  (W094A) 

Sub-Basin  Identifiers.  Identifiers  PG1  through  PG9  refer  to  sub-basin  areas 
within  the  pre-eruption  Gumain  basin  as  shown  by  Figure  3.2.1.  The  prefix  letters  “PG" 
denote  the  pre-eruption  Gumain  basin,  and  the  following  numbers  identify  sequentially- 
numbered  sub-basins.  These  identifiers  are  used  in  tables  and  figures.  The  convention 
generally  followed  for  sequential  numbering  of  subbasins  was  to  start  at  the  most  upstream 
sub-basin  and  to  end  at  the  most  downstream  subbasin. 
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Output  Site  Identifiers.  Identifiers  PG5US  through  PG9DS  refer  to  the 
simulation  output  sites  located  at  specific  points  within  the  pie-eroption  Gumain  basin  as 
shown  by  Figure  3.2.1.  The  first  two  characters  of  the  identifier  indicate  the  sub-basin 
within  which  the  site  is  located.  The  last  two  characters,  either  "US"  or  "DS”,  indicate  that 
the  site  is  located  at  either  the  upstream  or  downstream  end  of  the  sub-basin.  Calibration 
output  results  are  provided  in  tables  and  figures  only  for  the  site  corresponding  to  the.  stream 
gage  location,  this  being  the  most  downstream  site  in  each  calibration  basin. 

Tables  3.2.1  through  3.2.3  summarize  physical  and  computed  parameters  for  each  of  the 
calibration  basins.  Some  of  these  computed  parameters  (Unit  hydrograph  method,  storage 
coefficient,  loss  rate)  were  determined  through  the  calibration  process  as  described  in  the 
following  sections. 

Calibration  to  Historic  Event.  Calibration  of  the  HEC-1  runoff  model 
comprises  two  basic  steps: 

1)  Identifying  historic  events  for  which  concurrent  rainfall  and  runoff  data  are 
available  and  sufficient  to  describe  rainfall  over  the  basin(s)  as  well  as  the 
shape  and  volume  of  the  resultant  runoff  hydrograph  from  the  basin. 

2)  Running  the  HEC-1  model  with  known  (historic)  basin  rainfall  events  for 
input,  and  adjusting  model  parameters  until  the  model  suitably  reproduces  the 
known  (historic)  runoff  hydrograph  for  each  event. 

Because  most  of  the  basins  draining  Mount  Pinatubo  have  maximum  times  of  concentration 
of  the  order  of  three  to  12  hours,  hourly  data  were  desirable  for  calibration. 

HEC-1  model  calibration  for  basins  draining  Mount  Pinatubo  was  complicated  by  a  lack  of 
suitable  concurrent  rainfall  and  streamflow  data.  Available  stream  flow  data  are  often  of  poor 
quality  and  are  limited  to  average  daily  discharges  and  annual  peak  discharges.  Only  very 
limited  hourly  rainfall  data  are  available  from  any  of  the  rain  gage  stations. 

There  is  a  general  mismatch  between  available  streamflow  data  and  available  rainfall  data. 
Streamflow  data  are  very  sparse  after  1972,  and  only  limited  rainfall  data  are  available  prior 
to  1972.  Hourly  rainfall  data  could  be  obtained  for  only  a  small  number  of  storm  events 
from  three  stations  in  the  general  vicinity  of  Mount  Pinatubo:  Iba,  Hacienda  Luisita,  and 
Porac. 

Only  one  historic  event,  on  September  1,  1970,  was  found  with  (marginally)  sufficient 
concurrent  rainfall  and  streamflow  data  for  model  calibration.  The  data  consisted  of  hourly 
rainfall  at  Porac  and  streamflow  data  for  the  Gumain  River  at  Pabanlag  (W086A),  which 
included  daily  flows  and  a  peak  instantaneous  discharge  (NOTE:  In  the  publication 
Philippine  Water  Resources  Summary  Data,  Vol.  II  Streamflow  and  Lake  or  River  Stage, 
December  31,  1970,  the  peak  instantaneous  discharge  used  for  the  historic  event  calibration 
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is  reported  to  have  occurred  on  September  1 1  at  1:30  p.m.  However,  no  flows  o) 
significance  were  reported  on  September  1 1  a!  any  other  nearby  gages  as  were  reported  on 
September  1.  The  form  on  which  the  date  and  time  of  the  peak  were  recorded  did  not  have 
enough  spaces  to  record  a  2-digit  date  and  a  4-digit  time.  It  is  believed  that  the  peak 
reported  on  September  1 1  at  1:30  p.m.  actually  occurred  on  September  1  at  1 1:30  p.m.) 

The  s' ream  gage  oil  tlir  Gumain  River  at  Pabaniag  i»  reported  to  have  had  a  water  level 
recorder  hi  opetanon  during  this  event,  and  the  estimates  of  the  daily  average  flows  and  the 
peak  instantaneous  discharge  are  probably  more  reliable  than  ihose  from  any  other  nearby 
gage.  Although  a  water  level  recorder  was  apparently  in  operation  on  the  Gumain,  the 
recorder  strep  cnarts  cannot  be  located  and  hence  no  short-interval  (e.g.,  hourly)  flow  data 
are  available. 

The  calibration  event  was  the  Gumain  River  (W086A)  peak  annual  flow  of  267.5  m'/s  on 
September  1,  1970,  which  had  a  return  period  of  about  one  in  four  years.  Concurrent  hourly 
rainfall  data  were  available  from  a  gage  at  Porac,  located  about  20  km  east  of  headwater 
areas  of  the  Gumain  River  basin.  Figure  3.2.1  shows  the  relative  locations  of  the  basin  and 
the  Porac  rain  gage. 

Ideally,  hourfy  rainfall  data  for  model  calibration  would  be  available  from  rain  gages  within 
the  watershed,  or  at  least  from  several  gages  closely  surrounding  the  watershed.  This 
methodology  would  allow  for  an  accurate  estimation  of  the  spatial  and  temporal  distribution 
of  rainfall  over  the  basin.  Unfortunately,  hourly  data  were  available  for  only  the  Porac  rain 
gage,  and  there  were  no  daily  data  rain  gages  located  sufficiently  close  to  the  basin  to 
improve  estimation  of  actual  basin  rainfall. 

In  the  absence  of  additional  information,  hourly  rainfall  at  each  of  the  sub-basins  was 
estimated  as  being  hourly  rainfall  at  Porac  multiplied  by  an  adjustment  factor  to  account  for 
increasing  rainfall  with  elevation  as  shown  by  the  2-year  24-hour  isopluvtal  map  of  the  study 
area.  For  each  sub-basin,  the  adjustment  factor  was  computed  to  be  the  ratio  of  the  2-year 
24-hour  rainfall  at  the  middle  of  the  sub-basin  to  that  at  Porac.  By  this  approach,  rainfall 
amounts  on  the  sub-basin  varied  from  1.3  to  1.9  times  the  rainfall  at  Porac,  averaging  about 
1.6  times  the  Porac  rainfall. 

Figure  3.2.4  summarizes  key  results  from  three  simulations  for  the  historic  event.  Variables 
considered  in  these  simulations  were  loss  rate  and  the  Clark  storage  coefficient.  In  each  case 
these  parameters  were  adjusted  iteratively  until  the  historic  peak  discharge  was  reproduced. 

A  Clark  unit  hydrograph  storage  coefficient  of  15  to  25  times  the  time  of  concentration,  with 
appropriate  loss  rates,  provides  a  fairly  good  approximation  of  the  historic  hydrograph. 

When  considered  together  with  the  frequency-event  calibration  work,  a  storage  coefficient 
equal  to  25  times  the  time  of  concentration  and  a  constant  loss  rate  of  3  mm/hour  were 
adopted  for  subsequent  hydrologic  analyses  of  basins  throughout  the  study  area. 


Calibration  to  Frequency  Characteristics.  In  the  absence  of  other  suitable 
historical  data,  additional  calibration  efforts  were  limited  to  frequency  events.  The 
frequency-event-based  calibration  consisted,  for  example,  of  adjusting  HEC-1  model 
parameters  so  that  a  synthesized  100-year  design  storm  applied  over  a  gaged  basin  would 
generate  100-year  flow  volumes  ana  peak  instantaneous  discharge  at  the  gage  site. 

The  stations  chosen  for  the  frequency-event-based  calibration  were  the  same  three  chosen  for 
the  frequency  analyses:  the  Gumain  River  (W086A),  the  Bucao  River  (W093A),  and  the 
Santo  Tomas  River  (W094A).  The  pre-eruption  basins  above  these  gages  are  shown  by 
Figures  3.2.1  through  3.2.3.  Due  to  irrigation  diversions  after  1967,  post-1967  data  from 
the  Santo  Tomas  were  excluded  from  the  frequency-event-based  calibration. 

Frequency  analyses  of  streamflow  data  from  the  calibration  basins  are  described  in 
Section  2.4.4.  The  results  of  those  analyses,  in  light  of  the  available  record  of  uncertain  data 
and  the  sensitivity  of  the  analyses  to  single  data  points,  were  presented  as  families  of  flow 
frequency  curves  to  indicate  plausible  ranges  of  reasonable  discharges  and  volumes  at  various 
return  periods.  Ranges  of  reasonable  values  for  2-  and  100-  year  return  period  floods  were 
adopted  as  the  targets  to  be  attained  in  the  frequency  event  calibration. 

The  procedure  for  calibrating  the  two  available  parameters,  loss  rate  and  storage  coefficient, 
was  first  to  assign  a  storage  coefficient,  and  then  to  vary  the  loss  rate  until  the  target  value 
was  matched.  Tables  3.2.1  through  3.2.3  summarize  parameters  for  the  three  calibration 
basins.  Table  3.2.4  and  Figures  3.2.5  through  3.2.13  summarize  the  results  of  the  frequency 
event  calibrations,  showing  the  target  ranges  of  reasonable  values  together  with  the  results 
obtained  from  HEC-1  simulations  of  2-  and  100-year  events.  The  simulations  summarized 
are  based  on  a  one-hour  time  step  simulation  of  design  storms  of  five  days  in  duration. 

The  2-year  peak  instantaneous  discharge  for  the  Gumain  River  (W086A)  was  taken  as  the 
primary  target  value  for  calibration,  as  it  was  considered  to  be  the  most  reliable  of  possible 
targets.  Parameters  determined  for  the  Gumain  River  were  then  used  in  models  for  all  three 
calibration  basins  for  both  2-  and  100-year  hypothetical  storms. 

In  reviewing  the  final  calibration  run  outputs,  the  greatest  weight  was  given  to  matching  peak 
instantaneous  discharge,  then  to  three-day  volume.  The  frequency  data  and  curves  for  1-day 
volumes  were  given  the  least  weight  because  the  source  data  were  based  on  arbitrary 
calendar-day  periods  (and  often  on  staff  gage  readings)  which  would  tend  to  underestimate 
true  24-hour  maximum  values. 

Based  on  the  calibration  results,  it  was  decided  to  make  subsequent  tuns  of  the  HEC-1  model 
using  a  constant  loss  rate  of  3  mm/hr  and  the  Clark  unit  hydrograph  method  with  storage 
coefficient  computed  as  25  times  the  time  of  concentration  for  each  sub-basin.  These 
parameters  were  considered  to  yield  generally  good  results  for  both  peak  flows  and  total  flow 
volumes  for  all  three  calibration  basins.  The  simulated  100-year  peak  flow  for  the  Bucao 
River  appears  low  (Figure  3.2.8),  but  may  be  within  reason  when  considering  that  the  two 


A-41 


extreme  high  recorded  flows  may  have  actual  return  periods  in  excess  of  about  10  to 
20  years  as  suggested  by  the  plotting  positions  based  on  the  available  period  of  data  record. 

Sensitivity  Analysis.  Most  of  the  sensitivity  analysis  was  conducted 
concurrently  with  model  calibration  discussed  in  the  previous  sections  and  as  summarized  by 
Table  3.2.4  and  Figures  3.2.5  through  3.2.13.  That  analysis  consisted  of  varying  paired 
combinations  of  constant  loss  rate  and  Clark  unit  hydrograph  storage  coefficient. 

Increasing  (or  reducing)  the  loss  rate  within  reasonable  bounds  from  the  adopted  3  mm/hr 
value  has  a  relatively  large  impact  on  total  runoff  volume  and  a  relatively  small  impact  on 
peak  discharge.  A  typical  2-year  5-day  storm  near  the  summit  of  Mount  Pinatubo,  for 
example,  would  have  an  average  rainfall  intensity  of  about  6.5  mm/hr  considering  the  full  5- 
day  period,  and  a  peak  hour  rainfall  intensity  of  about  30  mm/hr.  Increasing  the  loss  rate 
from  3  to  6  mm/hr  for  determining  the  rainfall  excess  would  obviously  have  a  far  greater 
impact  on  the  total  runoff  volume  than  on  the  peak  hourly  flow. 

Increasing  (or  reducing)  the  Clark  unit  hydrograph  storage  coefficient  has  a  significant 
impact  on  the  peak  discharge,  but  no  impact  on  the  total  runoff  volume.  Increasing  the 
storage  coefficient  dampens  the  peak  flow  and  flattens  the  shape  of  the  hydrograph.  The 
impact  of  the  storage  coefficient  on  the  shape  of  the  hydrograph  is  illustrated  by 
Figure  3.2.4. 

With  the  adopted  unit  hydrograph  (Clark),  storage  coefficient  (25  times  time  of 
concentration),  and  loss  rate  (3  mm/hour),  model  results  were  found  to  be  generally 
insensitive  to  variations  in  model  time  step  and  routing  parameters.  Tests  using  the  Gumain 
Basin  of  varying  the  model  time  step  from  five  to  60  minutes,  the  Muskingum 
X  coefficient  from  0.2  to  0.4,  and  the  Muskingum  "K"  coefficient  within  a  range  reflecting 
flood  wave  velocities  of  1  to  5  m/s  all  yielded  peak  discharge  results  which  were  within 
about  2  percent  of  each  other.  Again,  this  lack  of  sensitivity  results  from  the  large  storage  in 
the  system  implied  by  the  laige  storage  coefficient  used  with  the  Clark  unit  hydrograph. 

3.2.4  HEC-1  Basin  Models 


Sub-basin  and  Output  Site  Numbering  Convention.  Sub-basins  are  designated 
by  a  2-  or  3-character  alpha-numeric  identifier.  The  first  character,  which  designates  one  of 
the  eight  major  basins,  is  one  of  the  following: 


Sacobia-Bamban  -  "S" 
Abacan  -  *A” 

O’Donnell  -  *0” 
Gumain/Porac  -  *G” 


Pasig-Potrero  -  "P” 
Santo  Tomas  -  "T" 
Bucao  -  "B” 

Maloma  -  "M” 


The  second  character,  which  is  numeric,  identifies  the  sub-basin  within  the  major  basin.  The 
convention  generally  followed  for  sequential  numbering  of  the  sub-basins  within  a  major 
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basin  was  to  start  at  the  most  upstream  sub-basin  and  to  end  at  the  most  downstream  sub¬ 
basin. 

Output  sites  are  designated  by  a  4-  or  5-character  alpha-numeric  identifier.  The  first  two  or 
three  characters  identify  the  sub-basin  in  which  the  site  is  located.  The  last  two  characters, 
either  “US"  or  “DS”,  indicate  that  the  site  is  located  at  either  the  upstream  or  downstream 
end  of  the  sub-basin. 

Structure  of  Basin  Models.  In  constructing  the  models,  sub-basins  were 
delineated  so  that  output  would  be  provided  for  at  required  locations.  Some  of  these  sub¬ 
basins  were  further  divided  to  define  major  tributary  streams  and  to  provide  better  definition 
of  rainfall  gradients  through  the  basins.  Additionally,  (1)  sub-basins  S7  in  the  Sacobia- 
Bamban  Basin,  P7  in  the  Pasig-Potrero  Basin,  and  G19  in  the  Gumain-Porac  Basin  were 
delineated  to  isolate  narrow  dike-confined  sections  of  the  channels,  and  (2)  sub-basin  T7  in 
the  Santo  Tomas  Basin  was  delineated  in  order  to  define  Lake  Mapanuepe. 

Definition  of  sub-basin  boundaries  to  provide  output  at  desired  sites  sometimes  resulted  in  the 
creation  of  some  very  small  basins.  These  areas  were  not  given  sub-basin  numbers  or 
explicitly  identified  in  the  HEC-1  models.  Instead,  for  modeling  purposes,  these  small  areas 
were  re-allocated  to  nearby  sub-basins  such  that  schematically  there  would  be  no  intervening 
areas  or  routing  reaches  (e.g.,  the  small  areas  just  upstream  of  sub-basin  S7  in  the  Sacobia- 
Bamban  Basin  (see  Figure  3.5.4)  were  re-allocated  to  upstream  sub-basins  S3,  S4,  and  S6 
such  that,  schematically,  simulation  sites  S3DS,  S4DS,  S6DS,  and  S7US  all  exist  near  site 
S7US  without  any  intervening  areas  or  routing  reaches).  Routing  reaches  through  these 
small  areas  were  ignored  because  flow  travel  times  are  significantly  shorter  than  the  one-hour 
model  time  step. 

Parameter  Estimation.  Physical  and  computed  parameters  at  simulation  output 
sites  were  obtained  as  described  beio.v. 

•  Stream  elevations  were  determined  from  1:50,000  topographic  maps  dated 
1986. 

•  Basin  areas  upstream  of  each  site  were  computed  as  the  sum  of  all 
contributing  sub-basin  areas. 

•  Basin  times  of  concentration  were  computed  using  the  formula  presented  in 
Section  3.2.2  based  on  the  slope  and  length  characteristics  of  the  entire  basin. 

•  Average  annual  basin  rainfall  amounts  were  estimated  based  on  isohyets  of 
mean  annual  rainfall  over  the  study  area  as  shown  by  Figure  2.3.4. 
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•  Average  annual  streamflows  at  each  site  were  estimated  on  the  basis  of 
average  annual  basin  rainfall,  basin  area,  and  a  formula  that  will  be  presented 
in  Section  3.3. 

•  In  some  basins  where  simulation  output  sites  do  not  share  a  common 
headwater  area,  more  than  one  storm  pattern  is  required  in  order  to  maximize 
flows  at  all  sites  within  a  given  basin.  Maximum  flows  at  each  site  were 
assumed  to  result  from  storms  centered  over  the  headwater  sub-basin  with  the 
highest  rainfall  as  determined  from  the  isohyetal  map. 

Physical  and  computed  parameters  for  sub-basins  were  obtained  as  described  below. 

•  Physical  parameters  of  area,  internal  flow  path,  channel  length,  and 
elevation  changes  were  measured  from  1 :50,000  sc<de  maps  dated  1986. 
Internal  flow  path  lengths  and  elevation  changes  are  provided  for  all  sub¬ 
basins,  as  required  'Vr  determining  time  of  concentration  for  local  sub-basin 
runoff.  Channel  lengths  through  sub-basins,  and  corresponding  elevation 
changes,  are  applicable  only  in  cases  where  there  is  •»  upstream  basin 
generating  incoming  channel  flow  to  be  routed  througi:  the  sub-basin. 

•  As  noted  in  Section  2.4.5,  the  ratio  of  sub-basin  event  rainfall  to  summit 
rainfall  was  determined  to  be  approximately  constant  for  all  return  periods  and 
for  all  durations  greater  than  six  hours.  For  convenience,  the  design  storm 
parameter  of  sui  oasin  rainfall  as  a  percentage  of  summit  rainfall  was 
determined  from  50-year  24  hour  rainfall  isopluvials  for  the  study  area  as 
shown  by  Figure  2.4.95.  This  parameter  is  used  to  correct  for  elevation- 
related  rainfall  variations  when  constructing  sub-basin  storm  hyetographs. 

•  The  iistance  of  each  sub-basin  from  the  center  of  the  design  storm(s)  is  used 
to  determine  the  depth-area  (depth-distance)  correction  when  constructing  sub¬ 
basin  storm  hyetographs.  These  distances  were  measured  from  the  middle  of 
the  sub-basin  with  the  storm  center  to  the  middle  of  each  other  sub-basin. 
Distances  to  other  storm  locations  are  given  only  if  the  sub-basin  contributed 
to  the  flows  at  the  sites  requiring  alternative  storm  locations. 

•  Runoff  parameters  of  time  of  concentration,  storage  coefficient,  infiltration 
loss,  and  base  flow  are  all  required  for  input  to  the  HEO  1  model  to  define 
runoff  from  each  sub-basin.  Times  of  concentration  were  computed  using  the 
formula  presented  in  Section  3.2.2  together  with  the  slope  and  length 
characteristics  of  each  sub-basin.  The  storage  coefficient  for  the  Clark  unit 
hydrograph  for  each  sub-basin  was  computed  as  25  times  the  time  of 
concentration.  Infiltration  rates  were  assumed  to  be  constant  at  3  mm/hour  for 
all  sub-basins  except  for  (1)  the  perched  S7,  P7,  and  G19  sub-basins  for  which 
an  artificially  high  1,000  mm/hour  loss  rate  was  applied  to  ensure  no  rainfall 


A-44 


excess  from  the  basin;  (2)  T7  (Lake  Mapanuepe)  for  which  no  infiltration  loss 
was  assumed;  and  (3)  0.8  km2  of  A2  and  2.6  km2  of  A4  which  were  assumed 
impervious  due  to  the  effects  of  urbanization.  The  high  loss  rate  was  applied 
to  the  perched  sub-basins  to  reflect  high  but  unquantified  losses  expected 
through  highly  permeable  perched  channel  reaches.  Base  flows  from  each  sub¬ 
basin  were  estimated  on  the  basis  of  average  annual  sub-basin  rainfall,  area, 
and  the  formula  presented  in  Section  3.2.2. 

Design  Storms.  The  2-,  10-,  50-,  100-,  and  500-year  storms,  developed  as 
described  in  Section  2.4.5,  were  used  for  the  HEC-1  modeling. 

Effects  of  Eruption.  Due  to  lack  of  justification  for  changes  to  any  other 
parameters  such  as  loss  rate  or  storage  coefficient,  eruption  impacts  reflected  in  the  HEC-1 
models  of  post-eruption  basins  were  limited  to  (1)  physical  changes  in  basin  areas,  (2)  the 
formation  of  Lake  Mapanuepe,  and  (3)  perched  channel  reaches  on  the  Bamban,  Potrero,  and 
Gumain  Rivers. 

3-2.5  Confidence  Limits  on  Computed  Frecuencv  Events.  The  5  and  95  percent 
confidence  limits  on  the  HEC-1  computed  instantaneous  peaks,  maximum  1-day  volumes, 
and  maximum  3-day  volumes  were  estimated  as  follows: 

1)  The  5  and  95  percent  confidence  limits  were  obtained  on  all  streamflow 
gage  data  sets  on  which  frequency  analyses  were  conducted  (see 
Section  2.4.4).  These  confidence  limits  were  obtained  in  accordance  with 
guidelines  described  in  Bulletin  17B  of  the  U.S.  Water  Resources  Council. 

2)  For  each  streamflow  gage  data  set,  the  change  in  percent  from  the  expected 
value  to  the  5  and  95  percent  confidence  limits  were  calculated  and  averaged. 
The  average  percentage  changes  obtained  from  the  streamflow  gage  data  sets 
were  applied  to  the  HEC-1  generated  instantaneous  peaks,  maximum  1-day 
volumes,  and  maximum  3-day  volumes  to  obtain  the  5  and  95  percent  limits. 

33  Flow  Duration  Curve  Modeling  Methodology 

3.3.1  Flow  Duration  Curve.  Flow  duration  curves  at  gaged  locations  were  presented 
in  Section  2.4.2  of  the  regional  analyses  portion  of  this  appendix.  Computed  flow  duration 
curves  at  hydrologic  simulation  output  sites  were  computed  on  the  basis  of  the  following  flow 
data  and  plotting  positions: 

1)  Average  annual  flow  above  each  simulation  output  site  was  estimated  from 
the  normal  annual  precipitation  (NAP)  map  of  the  study  area.  Rainfall  shown 
by  that  map  for  the  mountain  watershed  areas  had  been  estimated  using  a 
water  balance  approach  to  match  observed  average  flows  from  the  watersheds. 
Thus,  for  the  mountain  watersheds,  the  NAP  map  is  functionally  equivalent, 
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after  correction  for  evapotranspiration,  to  a  map  of  mean  annual  runoff  (and 
flows)  over  the  study  area. 


Conversion  from  average  annual  basin  rainfall  shown  by  the  NAP  map  to 
average  annual  flows  for  specific  basins  in  the  vicinity  of  Mount  Pinatubo  was 
accomplished  with  the  following  approximate  formula  derived  from  data  at 
gaged  basins: 

Q„t  =  A  (R/31,500  -  0.05) 

where:  Q„,  =  Average  annual  discharge,  m’/s 

A  -  Basin  area,  knf 

R  =  Average  annual  rainfall  over  basin,  mm/yr. 

In  deriving  the  formula  (and  also  the  NAP 
map),  average  annual  basin  rainfall  was 
assumed  to  be  equal  to  average  annual  runoff 
plus  1,550  mm  for  evapotranspiration  losses. 

Normalized  daily  flow  duration  curves  from  gaged  basins  indicate  that  the 
average  annual  flows  are  exceeded  from  15  to  40  percent  of  the  time, 
averaging  about  25  percent  (see  Figure  2.4.52).  For  the  computed  flow 
duration  curves,  the  average  annual  flow  was  plotted  to  be  exceeded  about 
25  percent  of  the  time. 

2)  The  2-year  24-hour  (average)  flow  was  computed  by  HEC-1  simulations. 
The  2-year  flow  is  expected  to  occur  once  every  two  years,  on  average,  based 
on  frequency  characteristics  of  an  annual  series  winch  considers  only  the  single 
highest  data  point  for  each  year  of  record.  However,  flow  duration  curves 
assume  a  partial  duration  series  which  consider  all  data  points  within  each 
year.  According  to  Kite12,  the  2-year  return  period  from  an  annual  series  is 
equivalent  to  a  1.44-year  return  period  in  a  partial  duration  series.  For  the 
computed  flow  duration  curves,  the  2-year  24-hour  (average)  flow  was  plotted 
to  be  exceeded  about  0.19  percent  of  the  time,  corresponding  to  one  day  in 
1.44  years. 

3)  The  10-year  24-hour  (average)  flow  was  computed  by  HEC-1  simulations. 
The  10-year  flow  is  expected  to  occur  once  every  10  years,  on  average.  No 
adjustment  was  made  for  the  minor  difference  in  partial  duration  vs.  annual 
series  annual  return  period,  9.5  years  vs.  10  years.  For  the  computed  flow 
duration  curves,  the  10-year  24-hour  (average)  flow  was  plotted  to  be 


Kite,  G.W.,  1977,  Frequency  and  Risk  Analysis  in  Hydrology,  Water  Resources  Publications. 
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exceeded  about  0.027  percent  of  the  time,  corresponding  to  one  day  in  10 
years. 

Flows  and  plotting  positions  were  interpolated  between  and  extrapolated  beyond  the  three 
points  described  above  considering  the  shape  of  flow  duration  curves  for  gaged  streams  in 
the  study  area,  as  shown  by  Figures  2.4.52  and  2.4.53. 

3.3.2  Confidence  Limits  on  Flow  Duration  Curves.  The  5  and  95  percent 
confidence  limits  on  the  flow  duration  curves  were  estimated  as  follows: 

1)  The  5  and  95  percent  confidence  limits  of  the  2-  and  10-year  24-hour 
(average)  flow  and  the  mean  annual  flow  were  obtained  on  all  streamflow  gage 
data  sets  on  which  frequency  analyses  were  conducted  'see  Section  2.4.4). 
These  confidence  limits  were  obtained  in  accordance  with  guidelines  described 
in  Bulletin  17B  of  the  U.S.  Water  Resources  Council. 

2)  For  each  streamflow  gage  data  set,  the  change  in  percent  from  the  expected 
value  to  the  5  and  95  percent  confidence  limits  were  calculated  and  averaged. 
The  average  percentage  changes  obtained  from  the  streamflow  gage  data  sets 
were  applied  to  the  HEC-1  generated  2-  and  10-year  24-hour  (average)  flow 
and  to  the  mean  annual  flow  that  was  obtained  as  described  in  Section  3.3.1. 
This  resulted  in  the  5  and  95  percent  confidence  limit  on  the  2-year  24-hour 
flow,  10-year  24-hour  flow,  and  average  annual  flow.  The  5  and  95  percent 
confidence  limit  flows  thus  obtained  were  used  to  develop  the  5  and  95  percent 
confidence  limit  flow  duration  curves  using  the  same  methodology  described  in 
Section  3.3.1. 

3.4  River  Hydraulic  Modeling  Methodology 

Flow  depths  and  velocities  were  required  for  the  hydraulic  design  of  mitigation  measures  in 
seven  of  the  eight  major  basins.  These  depths  and  velocities  were  obtained  through  use  of 
the  HEC-2,  Water  Surface  Profile  Model,  normal  depth  calculations,  and  critical  depth 
calculations.  Channel  and  cross-section  variations,  transport  of  sediment,  bed  and  bank 
roughness,  and  spill  resistance,  all  of  which  create  turbulence  and  energy  losses,  tend  to 
increase  with  increasing  discharge.  These  energy  losses  result  in  flow  conditions  on  steep 
natural  streams  that  may  approach  but  do  not  exceed  critical  flow  except  in  very  localized 
areas  of  the  channel.”  On  reaches  that  modeling  results  indicated  were  of  supercritical 
slope,  velocities  were  obtained  from  a  supercritical  flow  analysis  and  depths  from  a 
subcritical  flow  analysis.  The  velocities  and  depths  thus  obtained  were,  in  most  cases, 
approximately  equal  to  the  critical  depth  and  velocity  because  most  of  these  reaches  were 


“  Jarrett,  R.D.,  November  1984,  Hydraulics  of  High-Gradient  Streams,  Journal  of 
Hydraulic  Engineering,  Vol  110,  No.  11,  ASCE. 
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near  a  critical  slope.  On  reaches  that  modeling  results  indicated  were  of  subcritical  slope, 
velocities  and  depths  were  obtained  from  a  subcritical  flow  analysis.  A  Mannings’  *n”  value 
of  0.25  was  used  for  all  calculations. 

c.oss-sectional  data  for  both  HEC-2  modeling  and  normal/critical  depth  calculations  were 
obtained  from  a  number  of  sources.  In  a  few  instances,  detailed  surveyed  cross-sectional 
data  were  available.  In  many  instances,  cross-sections  were  obtained  from  a  digital  terrain 
model  (DTM).  At  some  locations  that  lacked  adequate  surveyed  or  DTM  data,  the  data 
required  for  input  to  the  HEC-2  model  were  obtained  from  aerial  photographs,  channel 
centerline  profiles,  and  professional  judgement  from  team  members  that  observed  the  site  of 
interest.  At  some  locations  no  depths  and  velocities  were  obtained  for  existing  conditions 
because  adequate  data  could  not  be  obtained. 

3.5  Hydrologic  Results 

3.5.1  Unit  Hvdrogiaphs.  Table  3.5.1  presents  the  peak  discharge  of  the  unit 
hydrograph  for  each  sub-basin.  Hie  unit  hydrograph  peaks  correspond  to  a  rainstorm  of  one 
hour  duration  with  a  total  rain  depth  of  1  mm.  Figure  3.5.1  shows  a  typical  unit  hydrograph 
computed  by  the  HEC-1  program.  Most  unit  hydrographs  are  very  similar  to  the  shape 
typified  in  Figure  3.5.1.  The  exceptions  are  presented  in  Figures  3.5.2  and  3.5.3.  These 
plots  for  sub-basins  A2  and  A4  show  composite  unit  hydrographs  representing  the 
combination  of  the  SCS  unit  hydrograph  (with  no  storage  coefficient)  applied  for  effectively 
impetvious  areas  and  the  Clark  unit  hydrograph  applied  for  non-impervious  areas.  These 
two  figures  dramatically  illustrate  both  the  peak  flow  increases  expected  from  urban 
development  in  the  Pinatubo  region,  and  the  effect  of  a  large  storage  coefficient  on  the  shape 
of  a  unit  hydrograph. 

3.5.2  Sacobia-Bamban  Basin.  The  Sacobia-Bamban  basin  is  146  knr  in  area, 
extending  northeasterly  from  the  base  of  Mount  Pinatubo  to  the  interior  lowlands  of  the 
island  of  Luzon.  Plate  1  shows  the  location  of  the  basin  relative  to  the  study  area  and 
hydrometeorological  data  stations.  Figure  3.5.4  shows  the  basin  at  a  larger  scale. 

The  basin  headwater  area  consists  of  steep  and  narrow  parallel  valleys  drained  by  the 
Sacobia,  Sapang-Cauayan,  Marimla,  and  Malago  Rivers.  Of  these,  only  the  Sacobia  and 
Malago  extend  to  near  the  base  of  Mount  Pinatubo;  the  other  rivers  originate  at  lower 
elevations  down  the  mountain’s  northeastern  slope.  The  Bamban  River  begins  at  the 
confluence  of  the  Sacobia  and  Marimla  Rivers  about  25  km  northeast  of  the  crater  rim,  just 
upstream  of  the  Highway  3  road  crossing  near  the  village  of  Bamban. 

Elevations  for  the  Sacobia  River  and  other  headwater  tributaries  range  from  about 
1,100  meters  in  the  headwaters  of  the  Sacobia  and  Malago  Rivers  to  55  meters  at  the 
confluence  defining  the  start  of  tne  Bamban  River  above  the  Highway  3  crossing.  The 
Bamban  River  component  of  the  basin  is  relatively  flat,  dropping  only  about  23  meters  over 
its  12  km  long  reach.  Most  ot  the  Bamban  River  is  contained  within  a  diked  channel  section 
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which  is  now  perched  above  the  surrounding  topography.  Perching  of  the  Bamban  River  is  a 
consequence  of  significant  aggradation  resulting  from  the  June  1991  eruption.  This  condition 
presumably  leads  to  a  net  water  loss  resulting  from  percolation  through  the  channel  bed  and 
under  the  levee  system.  For  modeling  purposes,  all  event  rainfall  for  sub-basin  S7  was 
assumed  to  be  lost  to  infiltration  (i.e.,  the  infiltration  rate  was  set  to  an  arbitrarily  high  value 
greater  than  the  maximum  rainfall  rate). 

The  1:50,000  scale  maps  dated  1986  indicate  no  significant  population  centers  or  urban 
development  within  the  basin.  Clark  AFB  is  located  immediately  to  the  south  of  the  basin, 
but  does  not  affect  basin  hydrology. 

The  Sacobia-Bamban  basin,  located  on  the  northeast  slopes  of  Mount  Pinatubo,  is  in 
Pinatubo’s  rain  shadow  during  the  Southwest  Monsoon  or  rainy  season.  Annual  rainfall 
amounts  over  the  basin  vary  from  a  maximum  of  about  4.000  mm/yr  in  the  upper  headwater 
areas  near  the  summit  of  Mount  Pinatubo  to  a  minimum  of  about  1,800  mm/yr  at  the 
downstream  end  of  the  basin  in  the  interior  lowlands  of  the  island  of  Luzon.  Similar 
variations  in  rainfall  over  the  basin  are  expected  during  single  storm  events.  Figure  2.3.4 
shows  isohyets  of  average  annual  precipitation  over  the  study  area. 

Sub-Basin  and  Ourout  Site  Parameters.  Physical  and  computed  parameters  of 
sub-basins  within  the  Sacobia-Bamban  River  Basin  are  summarized  in  Table  3.5.2.  Physical 
and  computed  parameters  at  simulation  output  sites  in  the  Sacobia-Bamban  River  Basin  are 
summarized  in  Table  3.5.3.  As  discussed  in  Section  2.4.5,  the  flow  at  a  given  site  was 
obtained  by  centering  the  storm  over  the  contributing  sub-basin  that  has  the  highest  rainfall 
(usually  the  highest  elevation  headwater  sub-basin).  Because  the  sites  listed  in  Table  3.5.3 
do  not  all  share  common  headwaters,  the  storms  had  to  be  centered  over  three  sub-basins  in 
order  to  obtain  flows  at  all  sites.  The  column  in  Table  3.5.3  labeled  "Critical  Storm 
Location’'  identifies  the  Sacobia-Bamban  sub-basin  over  which  the  storm  was  centered  to 
obtain  flows  at  the  indicated  site. 

Design  Flood  Hvdrographs.  Design  flood  hydrographs  computed  by  the  HEC- 
1  model  at  Sacobia-Bamban  basin  simulation  output  sites  for  each  of  the  2-,  10-,  50-,  100-, 
and  500-year  return  period  hypothetical  storms  are  shown  by  Figures  3.5.5  through  3.5.11. 
The  hydrographs  presented  for  each  site  correspond  to  model(s)  with  the  storm(s)  centered 
over  the  sub-basin  location  identified  by  Table  3.5.3  r'  maximum  flows  at  each  site. 

Design  Discharge'’  _  ^encv  Curves.  Peak  discharge  and  maximum 

24-hour  an '  3-day  flow  volurr  -i  irom  each  of  the  2-,  10-.  50-,  100-,  and  500-year 
design  flood  hydrographs  r-  .mmarized  by  Table  3.5.4. 
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Flow  Duration  Curves.  Daily  flow  duration  curve  data  for  each  of  the 
Sacobia-Bamban  basin  simulation  output  sites,  computed  following  the  method  described  in 
Section  3.3,  are  summarized  by  Table  3.5.5  The  data  shown  for  site  S7DS  do  not  account 
for  seepage  losses  during  low-flow  periods  from  the  perched  upstream  channel,  and  probabiy 
over-estimate  the  low-flow  characteristics. 

3.5.3  Ahacan  River  Basin.  The  Abacan  basin  is  51  km'  in  area,  originating  about 
4  km  east  of  the  crater  rim  of  Mount  Pinatubo  and  extending  easterly  to  the  interior  lowlands 
of  the  island  of  Luzon.  Plate  1  shows  the  location  of  the  basin  relative  to  the  study  area  and 
hydrometeorological  data  stations.  Figure  3.5.12  shows  the  basin  at  a  larger  scale. 

The  basin  headwater  area  consists  of  two  steep  and  narrow  parallel  valleys  drained  by  the 
Abacan  River  and  one  major  tributary,  Sapang-Bayo  Creek.  The  basin  headwaters  originate 
on  Mount  Pinatubo’s  eastern  slope  at  elevations  about  1,000  meters  below  the  crater  rim. 
Sapang-Bayo  Creek  joins  the  Abacan  River  about  4  km  upstream  of  the  Highway  3  crossing 
and  about  2  km  south  of  Clark  Air  Base.  The  lower  portion  of  the  basin  below  Highway  3 
is  mostly  confined  within  dikes. 

Elevations  for  Abacan  River/Sapang-Bayo  Creek  range  from  abort  500  meters  in  the  upper 
headwater  areas  to-130  meters  at  Sapang-Bayo/Abacan  confluence  to  10  meters  at  the  end  of 
the  dike-confined  channel  section.  Unlike  the  Bamban  River,  the  dike-confined  channel 
section  is  not  perched  above  the  surrounding  landscape. 

The  1:50,000  scale  maps  dated  1986  indicate  significant  urban  development  and  "densely 
built  up’  areas  in  the  Abacan  basin  in  the  vicinity  of  Clark  AFB  and  Angeles  City.  Portions 
of  Clark  AFB’s  runways  and  hangers  extend  into  the  basin  and  also  presumably  drain  to  the 
Abacan  River.  The  extent  of  urban  development  relative  to  the  basin  size  is  believed  to  be 
sufficient  to  have  a  noticeable  impact  on  basin  hydrology. 

The  Abacan  basin,  located  on  the  eastern  slopes  of  Mount  Pinatubo,  is  in  Pinatubo’s  rain 
shadow  during  the  Southwest  Monsoon  or  rainy  season.  Annual  rainfall  amounts  over  the 
basin  vary  from  a  maximum  of  about  3,000  mm/yr  in  the  upper  headwater  areas  on  the 
slopes  of  Mount  Pinatubo  to  a  minimum  of  about  1 ,800  mm/yr  at  the  downstream  end  of  the 
basin  in  the  interior  lowlands  of  the  island  of  Luzon.  Similar  variations  in  rainfall  over  the 
basin  are  expected  during  single  storm  events.  Figure  2.3.4  shows  isohyets  of  average 
annual  precipitation  over  the  study  area. 

Sub-Basin  and  Output  Site  Parameters.  Physical  and  computed  parameters  of 
sub-basins  within  the  Abacan  River  Basin  are  summarized  in  Table  3.5.6.  Physical  and 
computed  parameters  at  simulation  output  sites  in  the  Abacan  River  Basin  are  summarized  in 
Table  3.5.7.  As  discussed  in  Section  2.4.5,  the  flow  at  a  given  site  was  obtained  by 
centering  the  storm  over  the  contributing  sub-basin  that  has  the  highest  rainfall  (usually  the 
highest  elevation  headwater  sub-basin).  Because  the  sites  listed  in  Table  3.5.7  all  share 
common  headwaters,  the  storms  had  to  be  centered  over  only  one  subbasin  in  order  to 
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obtain  flows  at  all  sites.  The  column  in  Table  3.5.7  labeled  'Critical  Storm  Location” 
identifies  the  Abacan  sub-basin  (Al)  over  which  the  storm  was  centered. 


Design  Flood  Hvdrogranhs.  Design  flood  hydrographs  computed  by  the  HEC- 
1  model  at  Abacan  basin  simulation  output  sites  for  each  of  the  2-,  10-.  50-.  1 00-.  and  500- 
year  return  period  hypothetical  storms  are  shown  by  Figures  3.5.13  through  3.5.15. 

Design  Discharge/Volume  Freouencv  Curves.  Peak  discharge  and  maximum 
24-hour  and  3-day  flow  volume  data  from  each  of  the  2-,  10- ,  50-,  100-.  and  500-year 
design  flood  hydrographs  are  summarized  by  Table  3.5.8. 

Flow  Duration  Curves.  Daily  flow  duration  curve  data  for  each  of  the  Abacan 
basin  simulation  output  sites,  computed  following  the  method  described  in  Section  3.3.  are 
summarized  by  Table  3.5.9. 

3.5.4  O’Donnell  Basin.  The  study  area  considered  herein  for  the  O’Donnell  basin 
includes  two  major  rivers,  the  O’Donnell  and  the  Bulsa.  The  O’Donnell  River  drains  the 
northern  slopes  of  Mount  Pinatubo,  and  has  a  basin  area  upstream  of  the  confluence  with  the 
Bulsa  of  about  266  knf .  The  Bulsa  River  primarily  drains  the  eastern  slopes  of  the  Zambales 
mountains,  and  has  a  basin  area  upstream  of  the  confluence  with  the  O’Donnell  of  about 
510  km".  The  entire  basin  extends  about  2  km  below  the  O’Donnell-Bulsa  confluence  and 
has  a  total  area  of  about  817  km'.  It  is  the  largest  of  all  basins  being  assessed  under  the 
present  work. 

Plate  1  shows  the  location  of  the  basin  relative  to  the  study  area  and  hydrometeorological 
data  stations.  Figure  3.5.16  shows  the  basin  at  a  larger  scale. 

Basin  headwater  areas  on  Mount  Pinatubo  consist  of  steep  and  narrow  parallel  valleys 
drained  primarily  by  the  O’Donnell,  Apalong,  and  Bangat  Rivers.  Of  these  three  tributaries, 
only  the  O’Donnell  sub-basin  extends  folly  to  Pinatubo’s  crater  rim  where  the  post-eruption 
elevation  is  about  1,200  meters.  The  Apalong  and  Bangat  Rivers  originate  from  a  secondary 
peak  on  Mount  Pinatubo  which,  witn  a  pre-emption  summit  elevation  of  about  1,500  meters, 
may  now  be  the  highest  point  on  the  mountain. 

Basin  headwater  areas  for  the  Bulsa  River  on  the  east  slopes  of  the  Zambales  Mountains 
reach  a  maximum  elevation  of  about  1,600  meters.  These  headwater  areas  include  numerous 
steep  and  narrow  stream-cut  valleys  which  seem  generally  less  entrenched  into  the  mountain 
slopes  than  those  on  Pinatubo. 

The  stream  elevation  at  the  confluence  of  the  Bulsa  and  O’Donnell  Rivers.  located  near  the 
downstream  end  of  the  basin  study  area,  is  about  40  meters. 
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The  1:50,000  scale  maps  dated  1986  indicate  several  small  population  centers  (O’Donnell, 
Santa  Lucia,  Moriones),  but  none  of  sufficient  size  or  scale  to  appreciably  affect  the 
hydrology  at  sites  of  interest  within  the  basin. 

The  O’Donnell  basin,  which  generally  drains  in  a  northeasterly  direction,  is  in  the  rain 
shadow  of  the  Zambales  Mountains  and,  to  a  lesser  extent,  of  Mount  Pinatubo,  during  the 
Southwest  Monsoon  or  rainy  season.  Annual  rainfall  amounts  over  the  basin  vary  from  a 
maximum  of  about  6,000  mm/yr  in  the  upper  headwater  areas  of  the  Bulsa  River  in  the 
Zambales  Mountains,  and  5,000  mm/yr  in  the  upper  headwater  areas  of  O’Donnell  River 
tributaries  draining  Mount  Pinatubo,  to  a  minimum  of  about  1,800  mm/yr  at  the  downstream 
end  of  the  basin  in  the  interior  lowlands  of  the  island  of  Luzon.  Similar  variations  in  rainfall 
over  the  basin  are  expected  during  single  storm  events.  Figure  2.3.4  shows  isohyets  of 
average  annual  precipitation  over  the  study  area. 

Sub-Basin  and  Output  Site  Parameters.  Physical  and  computed  parameters  of 
sub-basins  within  the  O’Donnell  River  Basin  are  summarized  in  Table  3.5.10.  Physical  and 
computed  parameters  at  simulation  output  sites  in  the  O’Donnell  River  Basin  are  summarized 
in  Table  3.5.11.  As  discussed  in  Section  2.4.5,  the  flow  at  a  given  site  was  obtained  by 
centering  the  storm  over  the  contributing  sub-basin  that  has  the  highest  rainfall  (usually  the 
highest  elevation  headwater  sub-basin).  Because  the  sites  listed  in  Table  3.5.1 1  do  not  all 
share  common  headwaters,  the  storms  had  to  be  centered  over  three  sub-basins  in  order  to 
obtain  flows  at  all  sites.  The  column  in  Table  3.5.11  labeled  "Critical  Storm  Location" 
identifies  the  O’Donnell  sub-basin  over  which  the  storm  was  centered  to  obtain  flows  at  the 
indicated  site. 


Design  Flood  Hvdrographs.  Design  flood  hydrographs  computed  by  the  HEC- 
1  model  at  basin  simulation  output  sites  for  each  of  the  2-,  10-,  50- ,  100-,  and  500-year 
return  period  hypothetical  storms  are  shown  by  Figures  3.5.17  through  3.5.27. 

Design  Discharge/Volume  Frequency  Curves.  Peak  discharge  and  maximum 
24-hour  and  3-day  flow  volume  data  from  each  of  the  2-,  10-,  50-,  100-,  and  500-year 
design  flood  hydrographs  are  summarized  by  Table  3.5.12. 

Flow  Duration  Curves.  Daily  flow  duration  curve  data  for  each  of  the 
O’Donnell  basin  simulation  output  sites,  computed  following  the  method  described  in 
Section  3.3,  are  summarized  by  Table  3.5.13. 

3.5.5  Gumain/Porac  Basin.  The  Gumain/Porac  basin  is  302  knr  in  area,  extending 
in  a  generally  southeasterly  direction  from  Mount  Pinatubo  to  the  Pampanga  Delta.  Plate  1 
shows  the  location  of  the  basin  relative  to  the  study  area  and  hydrometeorological  data 
stations.  Figure  3.5.28  shows  the  basin  at  a  larger  scale. 

The  Gumain/Porac  basin  includes  two  major  rivers,  the  Gumain  and  the  Porac.  The 
headwaters  of  the  Gumain  River  consist  of  steep,  well-incised  tributaries  originating  near  the 


A-52 


crater  rim  of  Mount  Pinatubo  and  along  the  ridge  which  extends  south  from  Mount  Pinatubo, 
separating  the  Gumain/Porac  basin  from  the  westerly-flowing  Santo  Tomas  tributaries.  The 
Gumain  River  flows  approximately  32  km  southeast  from  the  crater  rim  of  Mount  Pinatubo 
to  its  confluence  with  the  Porac  River  at  the  head  of  the  Gumain  Floodway.  Elevations 
within  the  basin  range  from  about  1,600  meters  on  the  ridge  line  approximately  4  km  south 
of  Mount  Pinatubo  and  about  1,200  meters  at  the  crater  rim  to  about  10  meters  at  the  head  of 
the  Gumain  Floodway. 

The  headwaters  of  the  Porac  River  originate  on  the  southeast  slopes  of  Mount  Pinatubo. 
approximately  5  km  southeast  of  the  crater  rim.  The  river  flows  east  and  then  south  some 
39  km  to  its  confluence  with  the  Gumain  River  at  the  head  of  the  Gumain  Floodway. 
Elevations  within  the  Porac  basin  range  from  1,150  meters  at  the  high  point  to  10  meters  at 
the  head  of  the  Gumain  Flood  way. 

The  lower  reaches  of  the  Gumain  and  Porac  Rivers  contain  a  number  of  major  irrigation  and 
flood  control  projects,  including  the  Gumain  Floodway.  One  major  aspect  of  these  projects 
was  the  diversion  of  the  Porac  River  into  the  Gumain  Floodway  system;  the  Porac’s  natural 
course  appears  to  be  in  a  channel  which  flows  about  4  km  north  of  the  floodway.  However, 
it  has  not  been  possible  to  obtain  information  on  the  exact  configuration  and  operating 
policies  for  these  projects  or  their  current  (post-eruption)  condition.  For  purposes  of 
hydrologic  modeling  of  post-eruption  conditions,  it  was  assumed  that  flood  flows  from  both 
the  Gumain  and  Porac  Rivers  will  be  directed  to  and  confined  within  the  Gumain  Floodway. 

The  Gumain  Floodway  begins  at  the  confluence  of  the  Gumain  and  Porac  Rivers  and 
continues  downstream  approximately  8  km  to  its  outlet  in  the  Pampanga  Delta  at  an 
approximate  elevation  of  5  meters.  The  floodway,  represented  as  sub-basin  G19,  has 
aggraded  significantly  since  the  eruption  and  is  now  perched  above  the  surrounding 
landscape.  This  condition  presumably  leads  to  a  net  water  loss  resulting  from  percolation 
through  the  channel  bed  and  under  the  levee  system.  For  modeling  purposes,  all  event 
rainfall  for  sub-basin  G19  was  assumed  to  be  lost  to  infiltration  (i.e,  the  infiltration  rate  was 
set  to  an  arbitrarily  high  value  greater  than  the  maximum  rainfall  rate). 

The  drainage  area  delineated  during  this  study  for  the  Gumain/Porac  basin  at  the  downstream 
end  of  the  Gumain  Floodway  is  302  knf.  This  is  less  than  the  370  knv  basin  area  previously 
published  for  a  stream  gage  located  along  the  floodway.  The  previously  published  value  is 
believed  to  have  included  areas  which  can  no  longer  drain  to  the  lower  Gumain  because  of 
extreme  channel  aggradation. 

The  1:50,000  scale  maps  dated  1986  indicate  several  small  population  centers  (i.e.  Pabanlag, 
Del  Carmen,  and  Santa  Rita),  but  none  are  of  sufficient  size  to  appreciably  affect  the 
hydrology  at  sites  of  interest  within  the  basin. 

Annual  rainfall  amounts  over  the  basin  vary  from  a  maximum  of  about  5,000  mrr./yr  in  the 
Gumain  River  headwater  region  on  Mount  Pinatubo  to  2,000  mm/yr  at  the  downstream  end 
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of  the  basin  at  the  western  edge  of  the  Pampanga  Delta.  Similar  variations  in  rainfall  over 
the  basin  are  expected  during  single  storm  events.  Figure  2.3.4  shows  isohyets  of  average 
annual  precipitation  over  the  study  area. 

Sub-Basin  and  Output  Site  Parameters.  Physical  and  computed  parameters  of 
sub-basins  within  the  Gumain/Porac  River  Basin  are  summarized  in  Table  3.5.14.  Physical 
and  computed  parameters  at  simulation  output  sites  in  the  Gumain/Porac  River  Basin  are 
summarized  in  Table  3.5.15.  As  discussed  in  Section  2.4.5,  the  flow  at  a  given  site  was 
obtained  by  centering  the  storm  over  the  contributing  sub-basin  that  has  the  highest  rainfall 
(usually  the  highest  elevation  headwater  sub-basin).  Because  the  sites  listed  in  Table  3.5.15 
do  not  all  share  common  headwaters,  the  storms  had  to  be  centered  over  two  sub-basins  in 
order  to  obtain  flows  at  all  sites.  The  column  in  Table  3.5.15  labeled  "Critical  Storm 
Location’  identifies  the  Gumain/Porac  sub-basin  over  which  the  storm  was  centered  to  obtain 
flows  at  the  indicated  site. 

Design  Flood  Hvdrographs.  Design  flood  hydrographs  computed  by  the  HEC- 
1  model  at  basin  simulation  output  sites  for  each  of  the  2-,  10-,  50- ,  100-,  and  500-year 
return  period  hypothetical  storms  are  shown  in  Figures  3.5.29  through  3.5.35. 

Design  Discharge/Volume  Frequency  Curves.  Peak  discharge  and  maximum 
24-hour  and  3-day  flow  volume  data  from  each  of  the  2-,  10-,  50-,  100-,  and  500-year 
design  flood  hydrographs  are  summarized  in  Table  3.5.16. 

Flow  Duration  Curves.  Daily  flow  duration  curve  data  for  each  of  the 
Gumain/Porac  basin  simulation  output  sites,  computed  following  the  method  described  in 
Section  3.3,  ate  summarized  in  Table  3.5.17.  The  data  shown  for  site  G19US  do  not 
account  for  seepage  losses  from  the  perched  upstream  channel  during  low-flow  periods  and 
probably  over-estimate  the  low-flow  characteristics. 

3.5.6  Pasig-Potrero  Basin.  The  Pasig-Potrero  basin  is  77  km-  in  are:-.,  originating  at 
the  Mount  Pinalubo  crater  rim  and  extending  first  in  an  easterly  direction  and  then,  further 
downstream,  in  a  southeasterly  direction  to  the  Pampanga  Delta.  Plate  1  shows  the  location 
of  the  basin  relative  to  the  study  area  and  hydrometeorological  stations.  Figure  3.5.36  shows 
the  basin  at  a  larger  scale. 

The  basin  headwater  area  Is  drained  by  five  streams:  the  Bucbuc,  Yangca,  Timbu,  and 
Papatac  Rivers,  and  a  stream  that  prior  to  the  eruption  was  the  uppermost  headwater  stream 
of  the  Sacobia  River.  The  Papatac  River  is  formed  at  the  confluence  of  the  Bucbuc  and 
Yangca.  The  Pasig  River  is  formed  at  the  confluence  of  the  Papatac  and  Timbu  Rivers. 
Below  the  former  site  of  the  Mancatian  Bridge,  the  Pasig’s  name  changes  to  Potrero. 
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Elevations  in  the  basin  range  from  about  1,200  meters  near  the  crater  rim  to  near  zero  at  the 
confluence  of  the  Potrero  River  with  the  Guagua  River.  The  Potrero  River  component, 
which  comprises  almost  half  of  the  basin  length,  is  relatively  flat,  dropping  about  100  meters 
over  its  18  km  length. 

The  1:50,000  scale  maps  dated  1986  do  not  indicate  any  urban  development  that  would 
significantly  affect  basin  hydrology.  The  town  of  Bacolor  near  the  downstream  end  of  the 
basin  has  been  affected  by  shallow  flooding  from  the  Potrero  River,  but  is  not  considered  to 
lie  in  the  basin  because  the  Potrero  River  is  isolated  from  the  surrounding  topography  by 
levees  and  because  the  channel  from  about  4-1/2  km  below  the  former  site  of  the  Mancatian 
Bridge  to  the  confluence  with  the  Guagua  River  is  perched  above  the  surrounding 
topography. 

The  Pasig-Potrero  basin,  located  on  the  eastern  slopes  of  Mount  Pinatubo.  is  in  Pinatubo's 
rain  shadow  during  the  Southwest  Monsoon  or  rainy  season.  Annual  rainfall  amounts  vary 
from  a  maximum  of  about  5,000  mm/yr  in  the  upper  headwater  areas  on  the  slopes  of  Mount 
Pinatubo  to  a  minimum  of  about  1,800  mm/yr  at  the  downstream  end  of  the  basin  in  the 
Pampanga  Delta.  Similar  variations  in  rainfall  over  the  basin  arc  expected  during  single 
storm  events.  Figure  2.3.4  shows  isohyets  of  average  annual  precipitation  over  the  study 
area. 


Sob-Basin  and  Output  Site  Parameters.  Physical  and  computed  parameters  of 
sub-basins  within  the  Pasig-Parero  River  Basin  are  summarized  in  Table  3.5.18.  Physical 
and  computed  parameters  at  simulation  output  s;tes  in  the  Pasig-Potrero  River  Basin  are 
summarized  in  Table  3.5.19.  As  discussed  in  Section  2.4.5,  the  flow  at  a  given  site  was 
obtained  by  centering  the  storm  over  the  contributing  sub-basin  that  has  the  highest  rainfall 
(usually  the  highest  elevation  ix^  1  water  sub-basin).  Because  the  sites  listed  in  Table  3.5.19 
do  net  all  share  common  headwaters,  the  storms  had  to  be  centered  over  two  sub-basins  ir. 
order  to  obtain  flows  at  ai:  sites.  The  column  in  Table  3.5.19  labeled  “Critical  Storm 
Location'  identifies  -he  Pssig-P'Xrero  sub-basin  over  which  the  stonn  was  centered  to  obtain 
flows  at  the  indicated  site. 

Design  Flood  Hvdrograchs.  Design  flood  hydrographs  computed  by  the  HEC- 
1  medei  at  Pasig-Potrero  simulation  output  sites  for  each  of  the  2-.  10- ,  50- .  100-,  and  500- 
year  return  period  hypothetical  storms  arc  shown  by  Figures  3.5.37  through  3.5.45. 

Design  Discharve'Volun'.c  Frequency  Curves.  Peak  discharge  and  maximum 
24-hour  and  3-day  flow  volume  data  from  each  of  the  2-.  i0-.  50-.  100-,  and  500  year 
design  flood  hydrographs  are  summarized  in  Table  3.5.20. 

Flow  Duration  Corves.  Daily  flow  duration  curve  data  ior  each  of  the  Pasig- 
Potrero  basin  simulation  sitev,  computed  following  the  method  described  in  Section  3.3,  are 
summarized  in  Table  3.5.21.  The  data  shown  for  site  F7DS  do  no’  account  for  seepage 


losses  from  the  perched  upstream  channel  during  low-flow  periods  and  probably  over¬ 
estimate  the  low-flow  characteristics. 

3.5.7  Santo  Tomas  Basin.  The  Santo  Tomas  basin  is  approximately  262  knv  in  area, 
extending  in  a  southwesterly  direction  from  Mount  Pinatubo  to  the  South  China  Sea.  Plate  1 
shows  the  location  of  the  basin  relative  to  the  study  area  and  hydrometeorological  data 
stations.  Figure  3.5.46  shows  the  basin  at  a  larger  scale. 

The  Santo  Tomas  River  system  incorporates  two  major  tributaries,  the  Marella  River  and  the 
Mapanuepe  River,  which  join  to  form  the  Santo  Tomas.  The  headwaters  of  the  Marella 
River  originate  near  the  crater  rim  of  Mount  Pinatubo  at  an  elevation  of  about  1 ,500  meters 
and  along  the  ridge  extending  south  from  Mount  Pinatubo  which  separates  the  Santo  Tomas 
basin  from  the  easterly  flowing  Gumain  River  tributaries.  The  Marella  River  drains  the 
southwest  slopes  of  Mount  Pinatubo  and  combines  with  the  Mapanuepe  River  at  an  elevation 
of  about  90  meters.  The  reach  length  from  the  confluence  of  the  Marella  and  Mapanuepe 
Rivers  to  the  crater  rim  is  about  28  km. 

The  headwaters  of  the  Mapanuepe  River  originate  near  the  divide  between  the  Santo  Tomas 
and  Gumain  basins  at  an  deration  of  around  1,000  meters.  The  Mapanuepe  River  sub-basin 
includes  a  large  mine  site,  a  mine  tailings  dam,  and  Lake  Mapanuepe.  Approximately 
4.2  km-  of  the  mine  she  does  not  contribute  surface  runoff  to  the  watershed  and  hence  was 
no;  included  in  the  hydrologic  model.  The  impoundment  behind  the  tailings  dam  is  small  in 
comparison  to  Lake  Mapanuepe  and  no  flow  routing  was  done  through  this  facility. 

Lake  Mapanuepe  was  formed  following  the  June  1991  eruption  of  Mount  Pinatubo  as  a  result 
of  blockage  of  the  Mapanuepe  River  outlet  by  recurrent  lahars  and  severe  aggradation  on  the 
Marella.  Under  current  conditions,  the  Mapanuepe  River  joins  the  Marella  River 
approximately  1.5  km  downstream  from  the  outlet  of  Lake  Mapanuepe.  The  surface  area  of 
Lake  Mapanuepe  at  the  invert  elevation  of  its  current  outlet,  121.6  meters,  is  about  8.0  knf; 
topographic  contours  show  that  the  lake  surface  area  would  be  about  17  km'  at  a  water  level 
elevation  of  140  meters.  Lake  storage  volumes  between  these  two  elevations  were  computed 
using  the  conic  method  for  reservoir  volumes  as  described  in  the  HEC-1  manual.  The 
following  stage/discharge  characteristics  were  used  for  modeling  purposes: 


Lake  Elevation  (ml  Outlet  Discharge  fm’/sl 


<  121.6 
121.6  to  129.6 
<  129.6 


0 

23.57d£10d/(10-t-2d)]0“fI 

2.36a{a/[26+(2.08)(d-8)]}'’“7 
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where:  d  =  overflow  depth  =  lake  elevation  - 121.6 
a  =  flow  area  =  80  +  (7.78+0.28d)(d-8) 


The  Santo  Tomas  River  begins  at  the  Marella-Mapanuepe  confluence.  The  Santo  Tomas  is 
joined  by  the  Santa  Fe  River  approximately  10  km  downstream  from  gage  W094A.  The 
Santo  Tomas  then  flows  a  further  12  km  through  coastal  lowlands  to  Highway  7,  and  an 
additional  1  km  to  the  South  China  Sea. 

The  1:50,000  scale  maps  dated  1986  indicate  several  small  population  centers  (San  Rafael, 
Dalanaon  and  Aglao),  but  none  are  of  sufficient  size  or  scale  to  appreciably  affect  the 
hydrology  at  sites  of  interest  within  the  basin. 

Annual  rainfall  amounts  over  the  basin  vary  from  a  maximum  of  about  5,000  mm/yr  in  the 
Marella  River  headwater  region  on  Mount  Pinatubo  to  a  minimum  of  about  3,600  mm/yr  at 
the  downstream  end  of  the  basin  near  the  coastline.  Similar  variations  in  rainfall  over  the 
basin  are  expected  during  single  storm  events.  Figure  2.3.4  shows  isohyets  of  average 
annual  precipitation  over  the  study  area. 

Sub-Basin  and  Output  Site  Parameters.  Physical  and  computed  parameters  of 
sub-basins  within  the  Santo  Tomas  River  Basin  are  summarized  in  Table  3.5.22.  Physical 
and  computed  parameters  at  simulation  output  sites  in  the  Santo  Tomas  River  Basin  are 
summarized  in  Table  3.5.23.  As  discussed  in  Section  2.4.5,  the  flow  at  a  given  site  was 
obtained  by  centering  the  storm  over  the  contributing  sub-basin  that  has  the  highest  rainfall 
(usually  the  highest  elevation  headwater  sub-basin).  Because  the  sites  listed  in  Table  3.5.23 
do  not  all  share  common  headwaters,  the  storms  had  to  be  centered  over  two  sub-basins  in 
order  to  obtain  flows  at  all  sites.  The  column  in  Table  3.5.23  labeled  "Critical  Storm 
Location"  identifies  the  Santo  Tomas  sub-basin  over  which  the  storm  was  centered  to  obtain 
flows  at  the  indicated  site. 

Design  Flood  Hvdrographs.  Design  flood  hydrographs  computed  by  the  HEC- 
1  model  at  basin  simulation  output  sites  for  each  of  the  2-,  10-,  50- ,  100-,  and  500-year 
return  period  hypothetical  storms  are  shown  by  Figures  3.5.47  through  3.5.51. 

Design  Discharge/Volume  Frequency  Curves.  Peak  discharge  and  maximum 
24-hour  and  3-day  flow  volume  data  from  each  of  the  2-,  10-,  50-,  100-,  and  500-year 
design  flood  hydrographs  are  summarized  by  Table  3.5.24. 

Flow  Duration  Curves.  Daily  flow  duration  curve  data  for  each  of  the  Santo 
Tomas  basin  simulation  output  sites,  computed  following  the  method  described  in 
Section  3.3,  are  summarized  by  Table  3.5.25. 
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3.5.8  Bucao  Basin.  The  Bucao  basin  is  656  km5  in  area,  extending  in  a  generally 
northwesterly  direction  from  Mount  Pinatubo  and  southwesterly  from  the  Zambales 
Mountains  to  the  South  China  Sea.  Plate  1  shows  the  location  of  the  basin  relative  to  the 
study  area  and  hydrometeorological  data  stations.  Figure  3.5.52  shows  the  basin  at  a  larger 
scale. 

The  Bucao  basin  incorporates  the  Bucao  River  and  its  two  major  tributaries,  the  Balin- 
Buquero  River  and  the  Balintawak  River.  The  central  portion  of  the  basin  includes  a  large 
area  of  relatively  flat  and  low-lying  terrain  nestled  between  the  mountains  which  define  the 
basin  perimeter  Mount  Pinatubo,  the  Zambales  Mountains,  and  the  coastal  mountains  located 
between  Mount  Pinatubo  and  the  South  China  Sea. 

The  headwaters  of  the  Bucao  River  originate  on  the  northwest  slopes  of  Mount  Pinatubo  2  to 
5  km  north  of  the  crater  rim  at  an  elevation  of  about  900  meters.  The  river  flows  in  a 
generally  westerly  direction  through  rugged  terrain  for  approximately  28  km  to  its  confluence 
with  the  Balintawak  River  at  an  elevation  of  about  50  meters.  The  Bucao  then  enters  a  broad 
flat  valley  and  continues  to  flow  west  approximately  4  km  to  its  confluence  with  the  Balin- 
Buquero  and  a  further  12  km  to  the  Highway  7  crossing.  The  Bucao  enters  the  South  China 
Sea  approximately  2  km  below  Highway  7. 

The  headwaters  of  the  Balin-Buquero  River  originate  to  the  south  of  the  Bucao  River 
headwater  areas  and  extend  to  the  crater  rim  of  Mount  Pinatubo  at  an  elevation  of  about 
1,500  meters.  The  Balin-Buquero  and  its  principal  tributaries  (such  as  the  Maronut)  drain 
the  western  slopes  of  Mount  Pinatubo  and  the  northeastern  slopes  of  the  coastal  mountain 
range  which  lies  between  Mount  Pinatubo  and  the  South  China  Sea.  The  Balin-Buquero 
flows  in  a  generally  northwesterly  direction  for  approximately  20  km  from  the  crater  rim  of 
Mount  Pinatubo  to  its  confluence  with  the  Maronut  River  at  an  elevation  of  about  90  meters. 
Below  the  confluence  with  the  Maronut,  the  Balin-Buquero  enters  a  broad  flat  valley  and 
continues  to  flow  northwest  for  a  further  12  km  to  its  confluence  with  the  Bucao  at  an 
elevation  of  about  40  meters.  The  drainage  area  of  the  Balin-Buquero  above  its  confluence 
with  the  Bucao  is  approximately  217  knr. 

The  headwaters  of  the  Balintawak  River  originate  to  the  north  of  the  Bucao  River  headwater 
areas  and  drain  the  southern  slopes  of  the  Zambales  Mountains  at  elevations  of  up  to 
1,670  meters.  The  Balintawak  River  flows  in  a  generally  southwesterly  direction  through 
rugged  terrain  for  approximately  20  km  to  its  confluence  with  the  Bucao  River  at  an 
elevation  of  90  meters.  The  drainage  area  of  the  Balintawak  upstream  of  its  confluence  with 
the  Bucao  is  approximately  166  km5. 

The  headwater  areas  of  the  Bucao  and  Balin-Buquero  Rivers  were  severely  disturbed  by  the 
June  1991  eruption  of  Mount  Pinatubo,  with  massive  deposits  of  pyroclastic  material  filling 
in  entire  river  channels  and  destroying  much  of  the  pre-eruption  drainage  system.  Post¬ 
eruption  changes  to  drainage  boundaries  occurred  to  most  of  the  Bucao  and  Balin-Buquero 
headwater  drainages.  The  Maronut  River  was  the  most  affected  by  the  eruption,  with  a 
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reduction  in  catchment  area  from  31  km3  to  1 1  km:  as  a  result  of  gross  change*  in  catchment 
topography. 

The  1:50,000  scale  maps  dated  1986  indicate  several  small  population  centers  (San  Juan, 
Poonbato  and  Maguiguis),  but  none  are  of  sufficient  size  or  scale  to  appreciably  affect  the 
hydrology  at  sites  of  interest  within  the  basin. 

Annual  rainfall  amounts  over  the  basin  vary  from  a  maximum  of  about  6,000  mm/yr  in  the 
Zambales  Mountains  and  5,000  mm/yr  in  upper  headwaters  on  Mount  Pinatubo  to 
3,800  mm/yr  in  the  coastal  lowlands.  Similar  variations  in  rainfall  over  the  basin  are 
expected  during  single  storm  events.  Figure  2.3.4  shows  isohyets  of  average  annual 
precipitation  over  the  study  area. 

Sub-Basin  and  Output  Site  Parameters.  Physical  and  computed  parameters  of 
sub-basins  within  the  Bucao  River  Basin  are  summarized  in  Table  3.5.26.  Physical  and 
computed  parameters  at  simulation  output  sites  in  the  Bucao  River  Basin  are  summarized  in 
Table  3.5.27.  As  discussed  in  Section  2.4.5,  the  flow  at  a  given  site  was  obtained  by 
centering  the  storm  over  the  contributing  sub-basin  that  has  the  highest  rainfall  (usually  the 
highest  elevation  headwater  sub-basin).  Because  the  sites  listed  in  Table  3.5.27  do  not  all 
share  common  headwaters,  the  storms  had  to  be  centered  over  six  sub-basins  in  order  to 
obtain  flows  at  all  sites.  The  column  in  Table  3.5.27  labeled  "Critical  Storm  Location” 
identifies  the  Bucao  sub-basin  over  which  the  storm  was  centered  to  obtain  flows  at  the 
indicated  site. 


Design  Flood  Hydrogranhs.  Design  flood  hydrographs  computed  by  the  HEC- 
1  model  at  simulation  output  sites  for  each  of  the  2-,  10-,  50- ,  100-,  and  500-year  return 
period  hypothetical  storms  are  shown  by  Figures  3.5.53  through  3.5.63. 

Design  Discharge/Volume  Frequency  Curves.  Peak  discharge  and  maximum 
24-hour  and  3-day  flow  volume  data  from  each  of  the  2-,  10-,  50-,  100-,  and  500-year 
design  flood  hydrographs  are  summarized  by  Table  3.5.28. 

Flow  Duration  Curves.  Daily  flow  duration  curve  data  for  each  of  the  Bucao 
basin  simulation  output  sites,  computed  following  the  method  described  in  Section  3.3,  are 
summarized  by  Table  3.5.29. 

3.5.9  Maloma  Basin.  The  Maloma  basin  is  150  km'  in  area,  originating  about  7  km 
southwest  of  Mount  Pinatubo  and  extending  in  a  westerly  direction  to  the  South  China  Sea. 
Plate  1  shows  the  location  of  the  basin  relative  to  the  study  area  and  hydrometeorological 
data  stations.  Figure  3.5.63  shows  the  basin  at  a  larger  scale. 

The  Maloma  basin  includes  two  major  rivers,  the  Maloma  River  and  the 
Gorongoro/Kakilingar  River,  which  join  before  discharging  into  the  South  China  Sea.  The 
basin  primarily  drains  the  coastal  mountains  to  the  west  of  Mount  Pinatubo;  drainage  of 
Mount  Pinatubo  itself  is  limited  to  the  extreme  eastern  headwaters  of  the  Maloma  River 
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which  extend  to  the  lower  southwest  slopes  ot  Mount  Pinatubo  at  an  elevation  of  only  about 
600  meters. 


The  Maloma  River  flows  west  from  Mount  Pinatubo  in  a  narrow  canyon  through  the  coastal 
mountain  range  which  lies  between  Mount  Pinatubo  and  the  South  China  Sea.  It  is  then 
joined  by  the  Gorongoro/Kakilingar  River  about  6  km  upstream  of  the  Highway  7  bridge  at 
an  elevation  of  less  than  10  meters.  Elevations  within  the  Maloma  basin  range  from  sea 
level  to  about  1 ,000  meters,  with  the  highest  elevations  occurring  within  the  coastal 
mountains. 

The  Gorongoro/Kakilingar  River  originates  entirely  from  the  coastal  mountains  to  the  west  of 
Mount  Pinatubo,  and  flows  westward  in  a  deep  narrow  valley  through  the  coastal  mountains. 
Elevations  within  the  Gorongoro/Kakilingar  catchment  range  from  about  800  meters  in  the 
upper  headwater  areas  to  less  than  10  meters  at  the  confluence  with  the  Maloma. 

The  1:50,000  scale  maps  dated  1986  and  1991  indicate  several  small  population  centers 
(Payodpod,  Maquineng,  and  Maloma),  but  none  are  of  sufficient  size  to  appreciably  affect 
the  hydrology  at  sites  of  interest  within  the  basin. 

Annual  rainfall  amounts  over  the  basin  vary  from  a  maximum  of  about  5,000  mm/yr  in  the 
coastal  mountains  to  a  minimum  of  about  4,000  mm/yr  in  the  low  lying  area  between  Mount 
Pinatubo  and  the  coastal  mountains,  and  along  the  coast  near  the  South  China  Sea.  Similar 
variations  in  rainfall  over  the  basin  are  expected  during  single  storm  events.  Figure  2.3.4 
shows  isohyets  of  average  annual  precipitation  over  the  study  area. 

Sub-Basin  and  Output  Site  Parameters.  Physical  and  computed  parameters  of 
sub-basins  within  the  Maloma  River  Basin  are  summarized  in  Table  3.5.30.  Physical  and 
computed  parameters  at  simulation  output  sites  in  the  Maloma  River  Basin  are  summarized  in 
Table  3.5.31.  As  discussed  in  Section  2.4.5,  the  flow  at  a  given  site  was  obtained  by 
centering  the  storm  over  the  contributing  sub-basin  that  has  the  highest  rainfall  (usually  the 
highest  elevation  headwater  sub-basin).  Because  the  sites  listed  in  Table  3.5.31  do  not  all 
share  common  headwaters,  the  storms  had  to  be  centered  over  three  sub-basins  in  order  to 
obtain  flows  at  all  sites.  The  column  in  Table  3.5.31  labeled  "Critical  Storm  Location” 
identifies  the  Maloma  sub-basin  over  which  the  storm  was  centered  to  obtain  flows  at  the 
indicated  site. 


Design  Flood  Hvdrographs.  Design  flood  hydrographs  computed  by  the  HEC- 
1  model  at  basin  simulation  output  sites  for  each  of  the  2-,  10-,  50-,  100-,  and  500-year 
return  period  hypothetical  storms  are  shown  by  Figures  3.5.65  through  3.5.68. 

Design  Discharge/Volume  Frequency  Curves.  Peak  discharge  and  maximum 
24-hour  and  3-day  flow  volume  data  from  each  of  the  2-,  10-,  50-,  100-,  and  500-year 
design  flood  hydrographs  are  summarized  by  Table  3.5.32. 
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Flow  Duration  Curves.  Daily  flow  duration  curve  data  for  each  of  the 
Maloma  basin  simulation  output  sites,  computed  following  the  method  described  in 
Section  3.3,  are  summarized  by  Table  3.5.33. 

3.5.10  Confidence  Limits. 

Frequency  Events.  Table  3.5.34  presents  the  2-,  10-,  50-,  100-,  and  500-year 
peak  discharges  estimated  with  the  use  of  the  HEC-2  model  and  the  5  and  95  percent 
confidence  limits  about  these  peak  discharges  at  etch  hydrologic  site  considered  in  this  study . 
On  Figure  3.5.69  are  plotted  the  HEC-1  estimated  peak  discharges  for  the  2-  through  500- 
year  events  and  the  5  and  95  percent  confidence  limits  about  these  peaks  at  a  typical 
hydrologic  site  (Site  09US  in  the  O’Donnell  River  basin). 

Table  3.5.35  presents  the  2-,  10-,  50-,  100-,  and  500-year  maximum  24-hour  volumes 
estimated  with  the  use  of  the  HEC-2  model  and  the  5  and  95  percent  confidence  limits  about 
these  volumes  at  each  hydrologic  site  considered  in  this  study.  On  Figure  3.5.70  are  plotted 
the  HEC-1  estimated  peak  discharges  for  the  2-  through  500-year  events  and  the  5  and 
95  percent  confidence  limits  about  these  peaks  at  a  typical  hydrologic  site  (Site  09US  in  the 
O’Donnell  River  basin). 

Table  3.5.36  presents  the  2-,  10-,  50-,  100-,  and  500-year  maximum  3-day  volumes 
estimated  with  the  use  of  the  HEC-2  model  and  the  5  and  95  percent  confidence  limits  about 
these  volumes  at  each  hydrologic  site  considered  in  this  study.  On  Figu.„  3.5.71  are  plotted 
the  HEC-1  estimated  peak  discharges  for  the  2-  through  500-year  events  and  the  5  and 
95  percent  confidence  limits  about  these  peaks  at  a  typical  hydrologic  site  (Site  09US  in  the 
O’Donnell  River  basin). 

Flow  Duration  Curves.  Table  3.5.37  presents  the  data  for  the  5  percent 
confidence  limit  about  the  flow  duration  curve  at  each  hydrologic  site  considered  in  this 
study.  Table  3.5.38  presents  the  data  for  the  95  percent  confidence  limit  about  these  same 
flow  duration  curves.  On  Figure  3.5.72  are  plotted  the  computed  flow  duration  curve  and 
the  5  and  95  percent  confidence  limit  about  the  computed  curve  at  a  typical  hydrologic  site 
(Site  S7DS  in  the  Sacobia-Bamban  River  basin). 
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3.6  River  Hydraulic  Modeling  Results 

Table  3.6.1  presents  Clearwater  flow  depths.  Tabic  3.6.2  presents  bulked  flow  (i.e., 
sediment  +  water  flow)  depths.1'  In  Table  3.6.3  are  Clearwater  flow  velocities 
Table  3.6.4  presents  bulked  flow  velocities.  Note'  contained  in  Table  3.6.5  provide 
information  on  each  reach  that  is  useful  for  the  correct  interpretation  of  the  data  in 
Tables  3.6.1  through  3.6.4. 

Figures  3.6.1  through  3.6.11  present  the  results  of  the  hydraulic  modeling  of  clearwatcr 
flows  at  bridges  located  in  the  reaches  indicated  on  Tables  3.6.1  through  3.6.5. 

Figures  3.6.12  through  3.6.22  present  the  results  at  these  same  bridges  but  for  flows  that 
have  been  increased  for  suspended  sediment. 

3.7  HEC-1  Input 

Enclosure  1  provides  a  listing  of  the  HEC-1  input  used  for  the  modeling  of  the  100-year 
event  in  the  Sacobia-Bamban  River  basin.  Storms  were  centered  over  three  different  sub¬ 
basins  of  the  Sacobia-Bamban  Basin  in  order  to  obtain  flow  estimates  at  all  required  sites  (see 
Section  3.5.2,  Sub-Basin  and  Output  Site  Parame'ers).  Therefore,  Enclosure  1  includes  three 
separate  HEC-1  input  files,  one  for  each  of  the  three  assumed  storm  centers.  HEC-1  input 
for  the  2-,  10-,  50-,  and  500-year  events  on  the  Sacobia-Bamban  Basin  are  identical  to 
Enclosure  1  except  the  values  on  the  PI  cards,  which  represent  incremental  storm 
precipitation,  are  different. 

The  HEC-1  input  for  other  basins  was  similar  to  the  input  provided  in  Enclosure  1. 


"  Technical  Appendix  B,  Sedimentation,  contains  information  on  flow  bulking  for  suspended 
sediment. 
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ENCLOSURE,  PLATE,  AND  EXHIBIT  A  FOR 
TECHNICAL  APPENDIX  A 
HYDROLOGY  AND  HYDRAULICS 


ENCLOSURE  1 

HEC-1  INPUT  FOR  THE  100-YEAR  EVENT  ON  THE  PASIG-POTRERO 

BASIN: 


ID  PASIG-POTRERO  RIVER:  100-YEAR  EVENT 

ID  POST-ERUPTION  CONDITIONS 

ID  US  ARMY  CORPS  OF  ENGINEERS 
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•DIAGRAM 
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IT  60  01JAN00  0100  241 
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100-Year  Storm,  5-Day  Event: 

Ratio  of  Storm  Volume  from  Isopluvial  Maps  (Sub-Basin  to  Summit)  :  .55 
Distance  from  Sub-Basin  to  Assumed  Stoim  Center  :  .0  km 
Distance  Assumed  for  Depth-Area  Correction:  0.92  km 
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End  1-1 


100-Year  Storm,  5-Day  Event: 

Ratio  of  Storm  Volume  from  Isopluvial  Maps  {Sub-Basin  to  Summit)  :  .93 
Distance  from  Sub-Basin  to  Assumed  Storm  Center  :  .0  km 
Distance  Assumed  for  Depth-Area  Correction:  1.7  km 
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100-Year  Storm,  5-Day  Event: 

Ratio  of  Storm  Volume  from  Isopluvial  Maps  {Sub-Basin  to  Summit)  :  .70 
Distance  from  Sub-Basin  to  Assumed  Storm  Center  :  4.2  km 
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End  1-2 
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100-Year  Storm,  5-Day  Event: 

Ratio  of  Storm  Volume  from  Isopluvial  Maps  (Sub-Basin  to  Summit)  :  .65 
Distance  from  Sub-Basin  to  Assumed  Storm  Center :  6.70  km 
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tOJ-Year  Storm,  5-Day  Event: 

Ratio  cf  Storm  Volume  from  Isopluvial  Maps  (Sub-Basin  U  Summit)  :  .55 
Distance  from  Sub-Basin  io  Assumed  Storm  Center  :  S.7  km 


Eaci  1-3 
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100-Year  Storm,  5-Day  Event: 

Ratio  of  Storm  Volume  from  Isopiuvia!  Maps  (Sub-Basin  to  Summit)  :  .50 
Distance  from  Sub-Basin  to  Assumed  Storm  Center :  9.7  km 
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End  1-4 


100-Year  Storm,  5-Day  Event: 

Ratio  of  Storm  Volume  from  Isopluvial  Maps  (Sub-Basin  to  Summit)  :  .42 
Distance  from  Sub-Basin  to  Assumed  Storm  Center  :  15.0  km 
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100-Year  Storm,  5-Day  Event: 

Ratio  of  Storm  Volume  from  isopluvial  Maps  (Sub-Basin  to  Summit)  :  .40 
Distance  from  Sub-Basin  to  Assumed  Storm  Center :  20.8  km 


PGP6 


PI 

2.97 

2.97 

2.97 

2.97 

2.97 

2.97 

2.97 

2.97 

2.97 

Pi 

2.97 

2.97 

2.97 

2.97 

2.97 

2.97 

2.97 

2.97 

2.97 

PI 

2.97 

2.97 

2.97 

2.S7 

2.97 

2.97 

2.97 

2.97 

2.97 

PI 

2.97 

2.97 

2.97 

2.97 

2.97 

2.97 

2.97 

2.97 

2.97 

PI 

2.97 

2.97 

2.97 

2.97 

2.97 

2.97 

2.97 

4.79 

4.79 

PI 

4.79 

4.79 

4.79 

4.79 

4.79 

4.79 

4.79 

4.79 

4.79 

Enel  1-5 
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100-Year  Storm,  5-Day  Event: 

Ratio  of  Storm  Volume  from  isopiuvia)  Maps  (Sub-Basin  to  Summit)  :  .39 
Distance  from  Sub-Basin  to  Assumed  Storm  Center  •  ,°.S  km 
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*  KM  For  the  following  hydrographs:  Stations  that  begin  with  a  "P" 

*  KM  (e.g.,  P1DS,  P5US.  etc.)  are  the  total  hydrograph  at  the 

*  KM  indicated  location  but  may  (on  the  upper  sub-basins)  also  be 

*  KM  a  sub-basin  runoff  hydiograph.  Stations  that  begin  with  an  *RF" 

*  KM  are  sub-basin  runoff  hydrographs  for  the  indicated  area  and  do  not 

*  KM  represent  a  hydrograpn  at  any  location.  Stations  that  begin  with 


End  1-6 


*  KM  an  ’RT'  are  the  resulting  hydrograph  after  being  routed  from  the 

*  KM  indicated  location  to  the  next  downstream  location  and  do  not 

*  KM  represent  a  hydrograph  at  any  location. 

*  KM 

*  KM  Two  storms  were  applied  over  P3  (Timbu).  One  to  obtain  the 

*  KM  hydrograph  at  P3DS,  the  other  for  use  in  obtaining  hydrographs 

*  KM  at  other  locations  since  a  storm  centered  over  PO  gives  greater 

*  KM  flow  at  these  other  locations  than  a  storm  centered  over  P3 


KKP3DS 

KM  Hydrograph  for  P3DS 
BA  6.0 
BF  0.3 
PR  PT 
LU  0,3 

ZW  A=PlNATUBO_DESiG\l  B=PASIG  P3DS  C=FLOW  F=100-YR  COMPUTED 
UC  1.0.25.5 


KK  PODS 

KM  Hydrograph  calculation  for  PODS. 

BA  21.3 
BF  3.0 
PR  P0 
LU  0.3 

ZW  A=PINATUBO_DESlGN  B=PASIG  PODS  C=FLOW  F=100-YR  COMPUTED 
UC  0.7,17.5 


KK  RTP0DS 

KM  Muskingum  routing  of  PODS  to  P1DS. 
RM  1..7..2 


KKRFP1 

KM  Sub-basin  runoff  calculation  for  P1DS. 

BA  9.3 

BF  0.7 

PR  PI 

LU  0,3 

UC  .8,20.8 


KKP1DS 


Enel  1-7 


KM  Hydrograph  calculation  for  P1DS. 

ZW  A=PINATUBO_DESIGN  B=PAS!G  P1DS  C=FLOW  F=1  DO-YR  COMPUTED 
HC2 


KKP2DS 

KM  Hydrograph  calculation  for  P2DS. 

BA  4.4 
BF  0.3 
PRP2 
LU  0,3 

ZW  A=PINATUBO_DESIGN  B=PASIG  P2DS  C=FLOW  F=100-YR  COMPUTED 
UC  .4,11.0 


KKP4US 

KM  Hydrograph  calculation  for  P4US. 

ZW  A=PINATUB0_DES1GN  B=PASIG  P4US  C=FLOW  F=100-YR  COMPUTED 
HC2 


KK  RTP4US 

KM  Muskingum  routing  from  P4US  to  P4DS 
RM  1..4..2 


KKRFP4 

KM  Sub-basin  runoff  calculation  for  P4. 

BA  3.1 

BF  0.1 

PRP4 

LU  0,3 

UC  .9,23.5 


KKP4DS 

KM  Hydrograph  calculation  for  P4DS. 

ZW  A=PINATUBO_DESIGN  B=PASIG  P4DS  C=FLOW  F=100-YR  COMPUTED 
HC2 


KKP3DS 

KM  Hyd.  calc,  for  P3DS  for  use  at  pts.  other  than  P3DS 
BA  6.0 
BF  0.3 
PRP3 


End  1-8 


LU  0,3 
UC  1.0,25.5 


KK  P5US 

KM  Hydrograph  calculation  for  P5US. 

ZW  A=PINATUBO_DESIGN  B=PASIG  P5US  C=FLOW  F=100-YR  COMPUTED 
HC2 


KK  RTP5US 

KM  Muskingum  routing  from  P5US  to  P5DS 
RM  1,1. 1..2 


KK  RFP5 

KM  Sub-basin  runoff  calculation  for  P5. 

BA  12.6 

BF  0.0 

PRP5 

LU  0,3 

UC  1.9.48.3 


KKP5DS 

KM  Hydrcgraph  calculation  for  P5DS. 

ZW  A=PINATUBO_DESIGN  B=PASIG  P5DS  C=FLOW  F=100-YR  COMPUTED 
HC2 


KK  RTP5DS 

KM  Muskingum  routing  from  P5DS  to  P6DS. 
RM  1  ,.5„2 


KK  RFP6 

KM  Sub-basin  runoff  calculation  for  P6. 

BA  17.7 

BF  0.0 

PR  P6 

LU  0,3 

UC  2.2,55. 


KK  P6DS 

KM  Hydrograph  calculation  for  P6DS 
HC2 


End  1-9 


KK  RTP6DS 

KM  Muskingum  routing  from  P6DS  to  P7DS. 
RM  2,1.56, .2 


KKRFP7 

KM  Sub-basin  runoff  calculation  for  P7. 

BA  2.4 

BF  0.0 

PRP7 

LU  0,1000 

UC  4.82,121 


KKP7DS 

KM  Hydrograph  calculation  for  P7DS 

ZW  A=?INATUB0_DESI6N  B=PASIG  P7DS  C=FLOW  F=100-YR  COMPUTED 
HC2 


72. 


Encl  Mo 
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TECHNICAL  APPENDIX  A 


EXHIBIT  A 

DATA  PRIOR  TO  THE  OCTOBER  1993  CHANGE 
IN  PASIG/SACOBIA  BASIN  CONFIGURATION 


A  headwater  basin  of  the  Sacobia-Bamban  River  {formerly  identified  as  sub-basin  SI) 
was  captured  by  the  Pasig-Potrero  River  in  October  of  1993  and  is  now  assigned  a  sub¬ 
basin  identifier  of  PO.  In  the  Sacobia-Bamban  Basin,  this  change  in  basin  configuration 
reduces  the  runoff  through  sub-basins  S2,  S3,  and  S4  thereby  reducing  the  runoff  at 
simulation  output  sites  S2DS,  S3DS,  S3DS  +  S4DS,  S7US,  and  S7DS.  In  the  Pasig- 
Potrero  Basin,  this  change  in  basin  configuration  increases  the  runoff  through  sub¬ 
basins  PI,  P4,  P5,  P6,  and  P7  thereby  increasing  the  runoff  at  simulation  output  sites 
P1DS,  P4US,  P4DS,  P5US,  P5DS,  and  P7DS.  Runoff  was  not  affected  through  sub¬ 
basins  or  at  simulation  output  sites  not  listed  above. 

Figures  A-1  through  A-4  are  flood  hydrographs  at  two  simulation  output  sites  on  the 
Sacobia-Bamban  River  and  two  simulation  output  sites  on  the  Pasig-Potrero  River  that 
were  affected  by  the  change  in  basin  configuration.  These  hydrographs  were 
computed  for  conditions  that  existed  prior  to  the  October  1 993  capture  of  a  Sacobia- 
Bamban  headwater  basin  by  the  Pasig-Potrero  River.  Relative  changes  in  the  hydrology 
of  these  sub-basins  and  other  affected  sub-oasins  can  be  estimated  by  comparing 
Figures  A-1  through  A-4  to  the  appropriate  post-October  1993  figures  described  in  the 
text  of  Technical  Appendix  A  (Figures  3.5.6,  3.5.1 1,  3.5.43,  and  3.5.45). 

Table  A-1  shows  flow  duration  curves  at  the  same  simulation  output  sites  indicated  cn 
Figures  A-1  through  A-4  computed  for  conditions  that  existed  prior  to  the  October 
1993  capture  of  a  Sacobia-Bamban  headwater  basin  by  the  Pasig-Potrero  River. 

Relative  changes  in  the  hydrology  of  these  sub-Dasins  and  other  affected  sub-basins 
can  be  estimated  by  comparing  the  data  in  Table  A-1  to  the  appropriate  post-October 
1993  data  in  tables  described  in  the  text  of  Technical  Appendix  A  (Tables  3.5.5  and 
3.5.21). 


1 


Prior  to  October  1993  Loss  of 
Headwater  Basin  to  Pasig  -  Potrero  River 
Computed  Flood  Hydrographs 
Sncobia-Bamban  Basin  at  Site  S3DS 


Tine  in  Dnys 


Prior  to  October  19.93  Loss  of 
Headwater  Basin  to  Pasig  -  Potrero  River 
Computed  1‘lood  Hydrographs 


October  1993  Capture  of  Sacobia  Headwater  Basin 
Computed  Flood  Hydrographs 
Pnsig-Potrero  Basin  at  Site  P5US 


puocos  jsd  SJ313UJ  oiqno  ‘^Sicuosiq 


Time  in  Days 


Prior  to  October  1993  Capture  of  Sacobia  Headwater  Basin 


time  in  Days 


Table  A-1 

Prior  to  October  1993  Change  in  Basin  Configuration 

Sample  Computed  Flow  Duration  Curves 
Sacobia-Bamban  and  Pasig-Potrero  Basins 
Data  in  cubic  meters  per  second 


%  OFTrME 

EXCEEDS) 

S3DS 

S7QS 

P5US 

P7DS 

100 

0 

0 

0 

0 

50 

1.5 

3.4 

0.5 

0.7 

25 

* 

• 

1.2 

1.7 

20 

4.8 

10.7 

1.6 

2.2 

10 

7.9 

17.8 

2.8 

3.7 

5 

11.8 

25.9 

5 

5.5 

2 

20.1 

43.4 

8.8 

9.3 

1 

29.8 

63.1 

12.6 

13.7 

0.5 

42.4 

88.3 

17.3 

19.4 

0.2 

68.1 

142 

26.5 

32.1 

0.1 

91.7 

195.4 

36.5 

45 

0.05 

115.4 

249.3 

46.6 

58.2 

0.02 

140.6 

305.9 

57.2 

71.8 

*  Not  calculated- 
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TECHNICAL  APPENDIX  A 
HYDROLOGY  AND  HYDRAULICS 


Table  1.3.1 

Conversion,  SI  to  English  Units 


TABLE  2.3.1 

Daily  Rainfall  and  Evaporation  Data 


i  or  purposes  01  tu.s  :,tuuy,  USA01  and  USA02  are  station  identifiers  given  to  the  daily  raingauges  operated  by  the  U.S.  Navy  and  U.S.  Air  Force. 


TABLE  2.3.2 
PHTVOL.CS  Rainfall  Stations 


HUH 

DATE' 

INSTALLED 

|  RG-i  (201) 

Ml  Cu3drado 

South 

1070 

July  1991 

j  RG-2  (202) 

BUGZ 

Southwest 

600 

July  1991 

!  RG-3  (203) 

PI2 

Nonheast 

640 

July  1991 

|  RG-4  (204) 

Ml  Culianan 

Northwest 

550 

Aug.  1991 

!  RG-S  (205) 

Gumain 

Southeast 

820 

Aug.  1991 

|  RG-6  (206) 

Sacobia 

East/Northeast 

510 

Sept  1991 

TABLE  23.3 
Hourly  Rainfall  Data 


m 

rL^tTfriife- 

^  JrObta3red 

D324 

Iba,  Zaiclats 

15*20’  N 

119*58’  E 

1972417-14  to  19724)3-10 
1976415-17  to  19764)5-31 
19ED4J74B  to  19804)7-11 

A016 

Hadenda  Lnisita, 
Tailac 

ms'Tj 

120*36' E 

1 970-03-23  to  15704)9-06 
19724)7-14  to  1972-03-10 
1974-03-1210.19744)8-19  j 

JPORAC’ 

Santa  Cmz.  Potac 
Ratnpanga 

15*05*  K 

12D*33'H 

1970458-28  to  1970-094)3  1 
1972-07-14  to  19704)84)5  { 
19744)8-12  to  19744)8-19 
19764)5-19  to  .1976415-29 
1977-ll-:i  10  1577-13-16 

Official  Station  ID-  unknown. 


TABLE  2.3.4 

Cross  Correlations  of  Daily  Rainfall  Data 


Station  Key 

ID _ File  Name  Station  Name 

1  D3245192.WAT  Iba,  Zambaies 

2  A0166892.WAT  Hacienda  Luisita,  Tarlac 

3  A0177492.WAT  CLSU.  Munoz,  Nueva  Ecija 

4  A0208892.WAT  Magalang,  Pampanga 

5  R3127792.WAT  Bai  Magalang,  Pampanga 

6  R313709I.WAT  Julian  Subd,  San  Fernando,  Pampanga 

7  R3147092.WAT  Masantol,  Pampanga 

8  R3157690.WAT  Camiling,  Tarlac 

9  R3167290.WAT  Mayantoc,  Tarlac 

10  R3187591.WAT  Palawig,  Zambaies 

11  R3197590.WAT  San  Felipe,  Zambaies 

12  R3207592.WAT  Sta  Rita  Hem  Sch,  Cabangan,  Zambaies 

13  R3227592.WAT  Znas,  San  Marcelino,  Zambaies 

14  USNRPMB.WAT  Cubi  Point  Naval  Air  Station 

15  CLARK.WAT  Clark  Atr  Force  Base 


TABLE  2.3.5 

Estimated  Annual  Pan  Evaporation  (mm) 
CLSU,  Munoz 


Year'" 

V"  .'/fan 

Evaporation 

•/ 

1974 

1855 

1975 

2084 

1976 

1972 

1977 

2109 

1978 

1731 

1979 

1842 

1980 

1968 

1981 

1796 

1982 

1921 

1983 

2153 

1984 

1895 

1985 

2010 

1986 

1852 

1987 

2006 

1988 

2008 

1989 

1943 

1990 

1914 

1991  I  2049 


Table  2.3.6 
Annual  Rainfall,  mm 


RAIN 

Iba 

Uitita 

CLSU 

Julian  Kasantol 

Sta  Rita 

Znas 

Cubl 

Clark 

CAGE 

RD324 

RA016 

RA017 

R0313 

R0314 

R0320 

R0322 

RUSA1 

R’JSA2 

YEAR 

1950 

2207 

1951 

3039 

1782 

1952 

4073 

1756 

1953 

2916 

2167 

1954 

2716 

1315 

195S 

1927 

1439 

1956 

3S31 

1761 

1957 

1441 

1958 

1732 

1959 

2500 

1097 

1960 

5088 

3829 

2331 

1961 

5481 

3794 

2079 

19(2 

3535 

1970 

1963  - 

3784 

3477 

2205 

1964 

*8594 

3075 

1966 

1965 

3739 

3539 

2054 

1966 

3747 

3302 

2773 

196-» 

5072 

4396 

2001 

1968 

1972 

1754 

2146 

1694 

1969 

3585 

1465 

2688 

1727 

1970 

4272 

1796 

2346 

1971 

2745 

les? 

1669 

1625 

2964 

2191 

1972 

4659 

3526 

3678 

3128 

4308 

4120 

1973 

3324 

1300 

1459 

13(8 

2572 

1604 

1974 

4124 

2384 

2526 

26*9 

2172 

4138 

2619 

1975 

2528 

2045 

1438 

1638 

2868 

1516 

1976 

4374 

2475 

2650 

2654 

2232 

4888 

4154 

4226 

2705 

1977 

3901 

1713 

1505 

1410 

1472 

4112 

4250 

3768 

1765 

1978 

5227 

2011 

1998 

2119 

2212 

5099 

5402 

2347 

1979 

3551 

1518 

1522 

1483 

1558 

4131 

3293 

4058 

1817 

1980 

3960 

1742 

1675 

1885 

2072 

3946 

2564 

258$ 

1742 

1981 

1692 

1788 

1393 

1580 

3233 

1775 

1982 

1562 

1789 

3446 

3857 

1389 

1983 

2120 

1305 

1347 

985 

1685 

2467 

2260 

2566 

1034 

1984 

4137 

2107 

2126 

2428 

4277 

3025 

4758 

1920 

1985 

4119 

2211 

2264 

5154 

3561 

4944 

2391 

1986 

4024 

1983 

2292 

2312 

4930 

3983 

4612 

2313 

1987 

2562 

1171 

1337 

875 

1864 

2849 

2590 

1446 

1988 

3874 

2042 

1060 

2222 

4025 

2984 

3303 

1807 

1989 

1524 

2056 

1381 

1989 

4133 

3424 

3556 

1971 

1990 

3509 

2294 

1884 

2575 

4833 

3688 

4245 

2298 

1991 

4021 

1760 

1334 

2038 

4355 

3928 

«  OBS 

34 

18 

17 

20 

19 

13 

14 

31 

41 

MEAN 

3832 

1840 

1932 

1802 

1982 

4162 

3547 

3S75 

1966 

STD 

1211 

536 

378 

(63 

426 

748 

703 

801 

526 

TABLE  1.3.7 

Annual  Runoff  from  Basin,  mm 


'Hi  2j;'5  JZS3  JW  2772  2(18  4514 

688  1270  674  928  824  824  811  1569  *01 


TABLE  2.4.1 
Strenmflow  Dam 


:d  (Into  exists  after  1972 

years  missing  wlililn  period  of  record. 

n  ID's  'nnD'  tire  unofficial  ID'S 


TABLE  2.4.2 

Peak  Discharges  In  in’/s,  Including  Records  of  Doublful  Quality 


MO.  otv. 


TABLE  2.4.3 

Normalized  Peak  Flows  In  m’/s  /  km1,  Including  Records  of  Doubtful  Quality 


1.1)  0.44  0.1)  ).»(  l.D  1.11  0.11  0.14  1,11  1.11  1.0)  1.(1  1.1) 

1.5?  0.45  0,57  1.49  2.94  1.14  0.04  O.Ot  1.14  1,39  3.91  0.93  4.14 


TABLE  2.4.4 

Peak  Discharges  in  m’/s,  Excluding  Records  of  Doubtful  Quality 


SfO.  DCV, 


TABLE  2.A.S 

Normalized  Peak  Flows  In  m*/s  /  km2,  Excluding  Records  of  Doubtful  Quality 

XtVM  SULIA  O'OONtfCt  O'OOmttl  urn?  camiuhq  e amium  pasio-p  pasio-p  pdpac  popac  cuhaim  tu  CUMAJM  CAUIAHAV  COLO  IACSIT  SUCAO  Sto  TOMAS  KALCMA 


M  OIS.  IS  11  7  6  9  6  1  4  15  25  24  33  IB  24  13  16  21  7 

MEAN  3.33  0.51  0.53  4.50  3.71  2.22  0  25  «  “3  1.46  **9  1.36  1.65  6.26  1.28  1.57  1.64  1.98  2.73 

STO.  OIV.  1.45  0.45  0.57  0.75  2.33  1.34  0.04  .)  1.14  7  0.68  0.95  4.00  0.31  0.74  0.96  0.89  2.85 


Stream? low  Data 


7ht  parlod  of  avallabla  dally  racord  la  gartarally  ahortar  than  tha  parlod  of  racord  for  ptak  Inatantanaoua  flow*. 
In  abianco  of  ataga  racordar,  watar  laval  datamlnad  by  ataff  gauga  raad  t  or  J  tf«M  par  day. 


Table  2.4.7 

Daily  Data  Rainfall.  Evaporation,  and  Streamflow  Stations 


Type 


io  # 


Station  Nue 


Record  Considered 


RAIN 


R0324  Iba,  Zanhales 

RAOii  Hacienda  Lulslta,  Tar lac 

**017  C1S0,  Munoz,  heava  Ecija 

R0313  Julian  s&bd.,  San  Fernando.  Parpanga 

*0314  Kasantol,  Paspanga 

*0320  Sta  Rita  Elea  Scb,  Caban gan,  Zasbales 

*6322  2nas,  San  Harcelino,  2axbales 

KJSA1  Cubi  Pciat  Naval  Air  Station 

R0SA2  Clark  Air  Force  Base 


1951  -  1992 
1968  -  1992 

1974  -  1992 
1970  -  1991 
1970  -  1992 

1975  -  1992 
1975  -  1992 
1958  -  1990 
1950  -  199’ 


SVAP 


EAQ17  CLSO,  Munoz,  Neava  Eel Ja 


1974  -  1992 


FLOW 


W010A 
W011A 
WO  118 
W012A 
W023A 
K023S 
WQ81A 
W082A 
K084A 
W086A 
N088A 
W092A 
W093A 
K.94A 
K099B 


Bulsa  River,  Villa  Agllpay;  (405  k»*) 
O'Donnell  River,  Palublob;  (240  kn*) 
O'Donnell  River,  Patllng;  (112  Job*) 

Bangat  River,  sta  Lucia;  (90  ks*i 
Canlllng  River,  Kanbalen;  (142  ks?) 
Canlling  River,  Poblaclon.*  (280  ks*) 
Paslg-Potrero  Rlv--,  Cabetlcan;  (242  k»*) 
Pasig  Potrero  River,  Hda  Dolores;  (28  kn») 
Porac  River,  Del  Caraen;  (111  k=*) 

Ccnaln  River,  Pabanlag;  (128  kn*> 

Colo  River,  San  Benito;  (76  kn*) 

Bagslt  River,  Daapal;  (68  k*») 

Bucao  River,  San  Joan;  (615  k**l 
Santo  7out  River,  Dalanavan;  (177  k*2) 
Maloaa  River,  Kaloaa;  (151  k*») 


1960  -  1972 

1965  -  1972 
1958  -  1967 

1966  -  1972 

1964  -  1972 
1957  -  1966 

1965  -  1969 

1966  -  1972 
1957  -  1972 
1957  -  1571 
1957  -  1971 
1960  -  1972 
1957  -  1973 
1957  -  1967 
1984  -  1991 


Table  2.4.9 

Frequency  Analysis  of  1-Day  Maximum  Annual  Rainfall 


Rain  Cages 
(Rainfall,  as) 

RD324  Iba 

RAO 16  Kda  Lulslta 

RAO 17  COSO 

R0313  Julian  subd. 

R0314  Hasantol 

R0320  Sta  Rita  Sch 

R0322  Znas 

RUSA1  Cubi  Point  KAS 
RXJSA2  Clark  AFB 


1-day  data 

RD324  E* 

RAO 16  Kda  Lulslta 
RAO 17  ClSU 
R0313  Julian  Subd. 
R0314  Kasantol 
R0320  Sta  Rita  Sch 
R0322  Znas 
RUSA1  CublPolnt  KAS 
RCSA2  Clark  ATS 


•turn  Period;  (Chance  of  Exceedance  In  any  one  year) 


Years  of 

2-years 

10-years 

50-years 

100-years 

500-years 

Record 

(50%) 

(10%) 

(2%) 

(1%) 

(0.2%) 

38 

247 

399 

545 

610 

772 

20 

123 

210 

292 

329 

420 

18 

127 

204 

286 

324 

425 

22 

144 

273 

417 

468 

678 

22 

151 

240 

331 

374 

482 

17 

231 

419 

608 

696 

915 

18 

200 

346 

504 

580 

783 

33 

246 

375 

489 

537 

653 

41 

140 

2S5 

376 

433 

580 

(above)  multiplied  by  1.13  to  obtain  24-hour  duration  aaounts 

38 

279 

451 

626 

689 

872 

20 

139 

237 

330 

372 

475 

18 

144 

231 

323 

366 

480 

22 

163 

308 

471 

551 

766 

22 

171 

271 

374 

423 

545 

17 

261 

473 

687 

786 

1034 

18 

226 

391 

570 

655 

885 

33 

278 

424 

553 

607 

738 

41 

158 

288 

425 

489 

655 

Table  2.4.10 

Frequency  Analysis  of  2-Day  Maximum  Annual  Rainfall 


Return  Period;  (Chance  of  Exceedance  In  any  one  year) 


Rain  Cages 

Years  of 

2-years 

10-years 

50-years 

100-years 

500-years 

(Rainfall,  ra) 

Record 

(50%) 

(10%) 

(2%) 

(1%) 

(0.2%) 

R0324  Iba 

38 

375 

600 

820 

920 

1173 

RA016  Pda  Lcisita 

20 

170 

297 

417 

470 

599 

RAC  17  CX.SU 

18 

168 

298 

447 

521 

724 

R0313  Julian  Subd. 

22 

2C3 

402 

625 

735 

1026 

R0314  Kasar.tol 

22 

219 

359 

494 

554 

70S 

R0320  Sta  Rita  Sch 

17 

369 

617 

851 

955 

1209 

R0322  2nas 

18 

323 

S24 

723 

ei< 

1044 

RUSAl  Cubi  Point  NAS 

33 

359 

543 

709 

78  2 

957 

RUSA2  Clark  ATB 

<1 

M6 

364 

543 

628 

850 

2-day  data 

(above)  asltlplied  by  1.04 

to  obtain  48’ 

■hour  duration  aaounts 

R3324  Iba 

38 

390 

624 

853 

957 

1220 

RA016  Kda  Luis It a 

20 

177 

309 

434 

489 

623 

RAO 17  CLSU 

18 

175 

310 

465 

542 

753 

R0313  Julian  Subd. 

22 

211 

418 

650 

764 

1067 

R0314  Kasantol 

22 

228 

373 

514 

576 

733 

R0320  Sta  Rita  sch 

17 

384 

642 

885 

993 

1257 

R0322  2nas 

18 

336 

545 

752 

847 

1086 

RUSAl  Cubl  Point  NAS 

33 

373 

565 

737 

813 

995 

RUSA2  Clark  ATB 

41 

204 

379 

565 

653 

884 

Table  2.4.11 

Frequency  Analysis  of  5-Day  Maximum  Annual  Rainfall 


Return  Period;  (Chance  of  Exceedance  in  any  one  year) 


Rain  Cages 

Years  of 

2-years 

10-year: 

s  50-years 

100-years 

500-years 

(Rainfall,  ca) 

Record 

(50%) 

(10%) 

(2%) 

(1%) 

(0.2%) 

RD324  Iba 

38 

590 

892 

1164 

1282 

1567 

RA016  Hda  Luisita 

20 

239 

428 

624 

715 

950 

RAO 17  CLSU 

ie 

241 

435 

660 

773 

1084 

R0313  Julian  Subd. 

22 

2e7 

567 

892 

1055 

1498 

R0314  Masaatol 

22 

317 

560 

808 

924 

1219 

R0320  Sta  Rita  Sch 

17 

597 

939 

12S6 

1396 

1740 

R0322  Znas 

18 

502 

811 

1126 

1274 

1652 

RUSAl  Ccbl  Point  HAS 

33 

570 

872 

1137 

1250 

1520 

RUSA2  Clark  AFB 

41 

271 

502 

765 

895 

1251 

b-t lay  data 

(above)  aultlplled  by  1.02 

to  obtain 

120-hour  duration  asounts 

(0324  Iba 

38 

t02 

910 

1187 

1308 

1598 

RA016  Kda  Luisita 

20 

244 

437 

636 

729 

969 

RAO 17  CLSU 

18 

246 

444 

673 

788 

1106 

R0313  Julian  Subd. 

22 

293 

578 

910 

1076 

1528 

R0314  Kasantol 

22 

323 

571 

824 

942 

1243 

R0320  Sta  Rita  Sch 

17 

609 

958 

1281 

1424 

1775 

R0322  Znas 

18 

512 

827 

1149 

1299 

1685 

RUSAl  Cubl  Point  NAS 

33 

581 

889 

1160 

1275 

1550 

KUSA2  Clark  AFB 

41 

276 

512 

780 

913 

1276 

Table  2.4.12 

Frequency  Analysis  of  iO-Day  Maximum  Annual  Rainfall 


Return  Period;  (Chance  of  Exceedance  in  any  one  year) 


Rain  Cages 

Years  of 

2-years 

10-years 

50-years 

100-years 

500-years 

(Rainfall.  sr>) 

Record 

(50%) 

(10%) 

(2%) 

(1%) 

(0.2%) 

RD324  Iba 

38 

834 

1291 

1710 

1895 

2344 

RAO 16  Kda  Luisita 

20 

342 

585 

818 

922 

1180 

RAO  17  CLS'J 

28 

338 

547 

766 

871 

1145 

R0313  Julian  Subd. 

22 

392 

801 

1270 

1755 

2127 

R0314  Kasantol 

22 

450 

760 

1076 

1224 

1603 

R0320  Sta  Rita  Sch 

17 

909 

1422 

1914 

2137 

2691 

R0322  Znas 

18 

729 

1108 

1463 

1622 

2012 

rusai  cabi  Point  HAS 

33 

860 

1289 

1646 

1794 

2136 

R0SA2  Clark  AFB 

<1 

365 

662 

983 

1137 

1546 

10-day  data  (above)  rsaltlplied  by  1.01  to  obtain  2 < 0-hoar  duration  asocnts 


RD324  Iba 

38 

842 

1304 

1727 

1914 

2367 

RAQ16  Hda  Luisita 

20 

345 

591 

826 

931 

1192 

RA017  CLSU 

12 

341 

552 

774 

880 

1156 

R0213  Julian  Subd. 

22 

396 

809 

1283 

1773 

2148 

80314  Kasantol 

22 

455 

768 

1087 

1236 

1619 

R0320  St  a ‘-Rita  Set 

17 

318 

1436 

1933 

2158 

2718 

R0322  Zras 

18 

736 

1119 

1473 

1638 

2032 

RUSAI  Cubi  Point  HAS 

33 

869 

1302 

1662 

1812 

2157 

RUSA2  Clark  AFB 

41 

369 

669 

933 

1148 

1561 

Table  2.4.13 

Frequency  Analysis  of  15-Day  Maximum  Annual  Rainfall 


Return  Period;  (Chance  of  Exceedance  In  «ny  one  year) 


Rain  cages 

Years  of 

2-years 

10-years 

50-years 

100-years 

500-ycar: 

(Rainfall,  =) 

Record 

(50%) 

<1C%) 

(2%) 

(1%) 

(0.2%) 

RD324  Zba 

38 

1020 

1561 

2014 

2204 

2645 

RAO 16  Kda  Luislta 

20 

<00 

67< 

353 

1083 

1417 

RAO 11  CtSU 

18 

<0< 

620 

839 

9<0 

1201 

R0313  Julian  Subd. 

22 

442 

897 

1<14 

1670 

2357 

Rtn<  Kasantol 

22 

528 

883 

1230 

1388 

1784 

RO320  Sta  Rita  Sch 

17 

1100 

1737 

2332 

2597 

3246 

R0322  Znas 

28 

9C5 

1373 

1750 

1903 

7251 

RUSAl  Cubl  Point  KAS 

33 

1030 

1528 

19<1 

2113 

2507 

ROSA?  Clark  AF8 

<1 

<27 

779 

1188 

1395 

1967 

Table  2.4.14 

Streamfiow  Data  Frequency  Analyses  Summary 
Expected  Probability  Discharges  in  m3/s 

(ranges  correspond  to  alternative  data  sets) 


Certain  River  (KC86A) 


Table  2.4.15 

Composite  Analysis  of  Rainfall  Frequency  Data 


Return  Period.-  (Chance  o t  Exceedance  It.  any  one  year) 


AVERAGE  RAINFALL 

2-years 

10-years 

50-years 

100-years 

SCO-years 

(S3)  FOR 

DURATION 

(50%) 

(10%) 

(2%) 

(1%) 

(0.2%) 

COASTAL  CASES 

24— HR 

261 

435 

606 

684 

882 

COASTAL  CAGES 

2— DAY 

371 

594 

807 

902 

1140 

COASTAL  CAGES 

S-DAY 

S76 

896 

1194 

1327 

1652 

COASTAL  CAGES 

ID-DAY 

841 

1290 

1700 

1881 

2319 

COASTAL  CAGES 

1 S-DAY 

1014 

1550 

2C09 

2204 

2662 

INTERIOR  CAGES 

24-KR 

15S 

267 

385 

440 

584 

INTERIOR  CAGES 

2-DAY 

199 

358 

525 

60S 

812 

INTERIOR  CAGES 

S-DAY 

276 

508 

765 

890 

1224 

INTERIOR  CAGES 

ID-DAY 

381 

678 

992 

1194 

1535 

INTERIOR  CAGES 

1 S-DAY 

440 

771 

1125 

1295 

1745 

Ratio  of 

Interior  to  Coastal 

Rainfall 

for  same  Frequency-Dus 

atlon  Events 

INTERIOR/COASTAL 

24— HR 

0.S9 

0.61 

0.63 

0.64 

0.66 

INTERIOR/COASTAL 

2-DAY 

0.53 

0.60 

0.65 

0.67 

0.71 

INTERIOR/ COASTAL 

S-DAY 

0.47 

0.S6 

0.64 

0.67 

0.74 

INTERIOR/COASTAL 

10-DAY 

0.45 

0.52 

0.56 

0.63 

0.66 

INTEF  :  JR/ COASTAL 

1 S-DAY 

0.43 

0.49 

0.55 

0.S8 

0.65 

INTERIOR/  COASTAL 

Average 

0.49 

0.S6 

0.61 

0.64 

0.68 

NOTES: 


1)  Ratio  of  Interior/Coastal  Mean  Annual  Rain:  1904  es  /  3779  cm  -  C-50 

2)  Peak  hlgh-elevation  rainfall  on  Mount  Plnatubo  estimated  to  be  approximately  1.4  tines  the 
average  coastal  rainfall  for  all  frequency-duration  asour.es.  Higher  rainfalls,  tc  approximately 
1.7  tines  the  average  coastal  rainfall,  are  estimated  to  occur  at  higher  elevatlors  coastal 
soon  tains  located  approximately  40  ka  north  of  Mount  Plnatubo. 


Cages;1  RD324  Iba 

R0320  St a.  Rita  Elea.  School 

R0322  Zcas 

RUSAl  Ccbl  Point  NAS 


Interior 

Cages:  RAO 16  Hda  Lulslta 

RA017  CLSU 

R0313  Julian  Subdivision 
R0314  Kasantol 
RUSA2  Clark  AF8 


2 


For 


Table  2.4.16 

Short -Duration  Rainfall  Frequency  Data 
Representative  Interior,  Coastal,  and  High-Elevation  Stations 

Data  in  mm 


Return  Period; 

(Chance  ci 

Exceedance  in 

any  one  year) 

2 -years 

10-years 

50-years 

100-years 

SOO-y- 

Station 

DURATION 

(50%) 

(10%) 

(2%) 

C%) 

(0.2 

LOKLAXO  INTERIOR 
Kda.  Lulsita 

1  HOUR 

52 

85 

114 

12b 

n/a 

Hda-  Lulsita 

2  HOURS 

60 

102 

139 

154 

n/a 

Hda.  Lulsita 

3  HOURS 

69 

1)5 

155 

132 

n/a 

Hda.  Lulsita 

6  HOURS 

84 

130 

131 

188 

r/a 

Kda.  Lulsita 

12  HOURS 

101 

139 

248 

233 

e/a 

Hda.  Lulsita 

24  HOURS 

123 

246 

350 

394 

n/a 

HDDKK SO  COAST 

Iba 

1  HOUR 

61 

83 

109 

119 

n/a 

Iba 

2  HOURS 

83 

126 

164 

180 

n/a 

Iba 

3  HOURS 

98 

152 

2C0 

22C 

n/a 

Iba 

6  HOURS 

138 

224 

300 

332 

n/a 

Iba 

12  HOURS 

183 

332 

533 

601 

n/a 

Iba 

24  HOURS 

238 

482 

693 

388 

n/a 

RICH  ELEVATIONS 
Baguio  City 

1  HOUR 

56 

80 

102 

111 

n/a 

Baguio  City 

2  HOURS 

80 

145 

203 

223 

n/a 

Baguio  City 

3  HOURS 

100 

193 

281 

313 

n/a 

Baguio  City 

6  HOURS 

153 

340 

500 

568 

n/a 

Baguio  City 

12  HOURS 

231 

528 

389 

900 

n/a 

Baguio  City 

24  HOURS 

319 

634 

985 

1113 

n/a 

X3»z»  extracted  froe  Hydrology  and  Flood  Forecast  Center,  fASASA,  1981:  "Rainfall  Ir.tenslty- 
Duxatlon-Freqaency  Data  of  tbe  Philippines"  Volin*  1,  First  Edition. 


Baguio  City  is  located  approximately  140  la  NNE  of  Kt .  Pinatubo,  at  elevation  1310  n  (4500 
feet). 


Table  2.4.17 

Analysis  of  Short-Duration  Rainfall  Frequency  Data 
From  Representative  Interior,  Coastal,  and  High-Elev.  Stations 
For  Same-Frequency-Duration  Events 


Return  Period;  (Chance  of  Exceedence  in  any  one  year) 


Ratio  of  Station* 

DURATION 

2-years 

(50%) 

10-year* 

(10%) 

50-years 

(2%) 

ICO-years 

(1%) 

AVERAC 

INTERIOR  TO  COAST 

Luisite/lba 

1  HOUR 

0.86 

0.98 

1.04 

1.06 

0.98 

Luisita/Xba 

2  HOURS 

0.T2 

0.81 

0.85 

0.86 

0.83 

Luisita/Xba 

3  HOURS 

0.7C 

0.75 

0.77 

0-78 

0.7$ 

Luisita/Xba 

6  HOURS 

0.61 

0.58 

0.57 

0.S7 

0.58 

Luisita/Xba 

12  HOURS 

0.54 

0.48 

0.46 

0.46 

0.48 

Luisita/Xba 

24  HOURS 

0.53 

0.51 

0.50 

0.50 

0.51 

INTERIOR  TO  HIGH  ELEVATIONS 
Luisita/fiaguio  1  HOUR 

0.94 

1.06 

1.12 

1.13 

1.06 

Luis lta/Begs io 

2  HOURS 

C-75 

0.70 

0.68 

0.68 

0.70 

Lui* ita/Baguio 

3  HOURS 

C.68 

C.58 

0.55 

0.54 

0.58 

Luis ita/Baguio 

6  HOURS 

0.53 

0.38 

0.34 

0.33 

0.39 

Luixita/Baguio 

12  HOURS 

O.M 

0.34 

0.31 

0.31 

0.35 

Luixita/Baguio 

24  HOURS 

0.40 

0.36 

0.36 

0.35 

0.36 

COAST  TO  RICH  ELEVATIONS 

Iba/Baguio  1  HOCR 

1.10 

1.C8 

1.07 

1.07 

1.07 

Iba/Baguio 

2  HOURS 

1.04 

0.87 

0.81 

0.79 

0.87 

Iba/Baguio 

3  HOURS 

0.38 

0.77 

0.71 

C.69 

0.78 

Iba/Baguio 

6  HOURS 

C.88 

0.66 

0.60 

0.58 

0.68 

Iba/Baguio 

12  HOURS 

0.81 

0.7C 

0.68 

0.67 

0.71 

Iba/Baguio 

24  HOURS 

0.75 

0.72 

0.71 

0.7i 

0.71 

Table  2.4.18 


Depth-Duration  Data  for  Design  Storms 
Hypothetical  Station  Located  near  Summit  of  Mount  Pinatubo 

Data  in  mm 


Return  Period;  (Chance  of  Exceedance  in  any  one  year) 


2-years 

10-yrs 

50-yrs 

100-yrs 

500-yrs 

DURATION 

(50%) 

(10%) 

(2%) 

(«) 

C.2X) 

1-Hour 

57 

81 

103 

113 

146 

2*  Hour* 

82 

140 

193 

21S 

280 

3-Hour* 

102 

186 

261 

293 

381 

6- Hours 

159 

311 

446 

505 

655 

12-Hours 

240 

478 

694 

789 

1022 

24-Hours 

352 

587 

850 

936 

1196 

2-Days 

501 

802 

1101 

1230 

1544 

5-Days 

777 

1209 

1624 

1804 

2236 

NOTES 

1-  The  depth-durat ion- frequency  data  summarized  above  are  for  a  hypothetical  station 
located  at  approximately  1000  meters  elevation  near  the  sunit  of  Notait  Pinatubo. 
This  station  has  a  mean  annual  rainfall  of  about  5000  mo. 

2.  The  1-hour  duration  depth -frequency  data  are  effectively  independent  of  station 
location.  The  1-hour  data  presented  above  are  directly  applicable  to  any  site 
(coastal,  mountain,  or  interior)  throughout  the  study  area. 

3.  Depth-frequency  data  for  durations  of  6  hours  and  longer  may  be  transposed  to  other 
sites  throughout  the  study  area  after  multiplication  by  a  site-specific  factor. 

The  factor  is  equal  to  the  ratio  of  the  desired  site's  mean  annual  rainfall  in 
millimeters  to  3000  wm,  the  mean  annual  rainfall  for  the  hypothetical  station. 

4.  Same- frequency  durations  are  not  fully  f^irdded  in  a  single  storm;  i.e.,  a  100- year 
24-hour  storm  is  not  expected  to  include  a  100-year  1-hour  duration.  The  available 
data  shov-d  a  maxioja  20-year  1-hour  rain  within  a  200-year  24-hour  storm,  and  a 
maxima  2-year  1-hour  rain  within  a  25-year  24-hour  storm.  This  relationship  of  1- 
hour  to  24-hour  reum  periods  can  be  approximated  by  asswing  a  constant  rainfall 
intensity  over  the  peak  six-hour  duration. 


Table  3.2.1 

Sub-Basin  Parameters 
Ouumiti  Basin,  Pre-Eruption 
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Plow  Routing  Parameters 

Travel  Tlmo  (iirs)  N/A  0.9  N/A  0,8  0.5  N/A  N/A 

Number  of  Sub-Rcaclics  N/A  1  N/A  I  1  N/A  N/A 


Table  3.2.2 

Sub-Basin  Parameters 
Bucao  Basin,  Pre-Emption 


Flow  Routing  Parameters 

Travel  Time  (lirs)  N/A  N/A  0.3  1.4  N/A  N/A 

Number  of  Sub-Reaches  N/A  N/A  0  1  N/A  N/A 


Table  3.2.3 

Sub-Basin  Parameters 
Santo  Tomas  Basin,  Pre-Emption 
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Plow  Routing  Parameters 

Travel  Time  (hrs)  ■  N/A  0.6  0.7  N/A 

Number  of  Sub-Reaches  N/A  1  I  N/A 


Table  3.2.4 


HEC-l  Calibration  and  Sensitivity  Analysis 
Simulations  with  Clark  Unit  Hydrograph 


Target 

180-230  ! 

490-670 

11-13 

29-36  1 

(n/a)-28 

8l-(n/a) 

Gumain  (W086A) 

A 

220 

940 

12 

54 

18 

78 

B 

220 

790 

16 

53 

28 

96 

C 

220 

680  ] 

17 

49 

36 

101 

Target 

890-960 

37504350 

5660 

108-235 

135-142 

230484 

Bucao  (W093A) 

A 

731 

3370 

52 

232 

93 

406 

B 

786 

2772 

64 

214 

140 

471 

C 

811 

2358 

67 

191 

16$ 

467 

Target 

310-360 

800-1560 

22-26 

45-100 

54-57 

145-280 

Santo  Tomas  (W094A) 

A 

320 

1580 

21 

102 

31 

174 

B 

330 

1190 

26 

84 

49 

364 

C 

320 

970 

27 

73 

60 

154 

•  Target  values  (range)  from  HEC-FFA  analyses  of  recorded  data;  ranges  correspond  to  alternative  data  sets. 


Constant 

Loss  Rate  Unit  Hydrograph 
mm/hr  Storage  Coefficient 


(T<  =  Basin  Time  of  Concentration) 


Table  3.5.1 

Sub-Basin  Unit  Hydrograph  Peaks 


!-«.  «D  rt  IT)  O 

cm  -r-  in  r->  c\j 

in  cj  -r-  co 

o'  o'  o'  o  d 


cm  co  tn  o  CD 

O  CD  CO  N  03 

*-  ©  o  o  o 

o  o  o  d  o 


TfO'CCDCMCDCOCONTf 

CDt-cONIOOCDO’-W 

o^cor-^nococoto 

dddddddddo 


O  CVJ  CO  ^  to  CO 
Q.  0.  0.  0.  D.  Q.  Q. 

N  _  , 

TTOOOOOOO  ! 

UOGffi®ffi®fflWCOMlDCD 

~i:^±:i=i:i=i5C/)C/)C0a>CO 
O  O  O  O  O  O  O  O  -  -  -  *  - 

c  Dl.  Q  Q  n.  n  r>u  to  B  c  g  ffl 

QajOTmwwwcgooooo 

—  —  ~  Q  Q  Q  0  Q 

CO  CO  CO  CO  <0 


o^iflmwoocD’-T-wcoi-cosiflco^toin 

CDC0U)0)10WCD0)C0Tra)C0OU)OC0®C0NN 

qq^-^^^^-r-q-r-i-cMCM-r-co^^^qq 

dddddddddo  dddddddddo 


o  •*-  CM 

§555 


cohn-to 

o  o  o 


CD  ffl  ffl  ffl 

c  c  c  c 
c  c  c  c 
o  o  o  o 
p  o  p  p 


o 

00  o>  TT 

O  O  O 


a>  ©  © 
c  c  c 
c  c  c 
o  o  o 

a  o  o 


OaT-CO©t“OCDC9©COOWO^©NS©»*WNr- 

OdJCDCMCOCOCOCD^-O^r^rCD^T^'CM-r-OlOlOTrCOCM 

©v-qcoq^rco©cotocM'»*;CM<r-u>Tsrc,>f-c\jtococMcu 

do"  dddddddddo*  o' dddddddddo 


o  o  o 
a  eg  © 
o  o  o 

3  3  3 

cd  m  co 


0.221  O'Donnell,  011 _ 0.383  Santo  Tomas,  T8 _ 0.181 

0.246  O'Donnell,  012 _ 0.224  Santo  Tomas,  T9 _ 0,160 

0.170  O'Donnell,  013 _ 0,447  Santo  Tomas,  T10  0.029 


Table  3.5.2 
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Table  3.5.4 

Computed  Maximum  Annua!  Peak  Discharge  and  Voiumc  Frequency  Data 
Saeotwi-Samban  Basin 


Site  I  D  Maximum  Annua!  Peak  Discharge  (cms) 


2-year 

10-year 

SO-year 

100-year 

500- year 

S2DS 

25 

54 

82 

95 

125 

S3DS 

32 

71 

109 

125 

166 

S4DS 

13 

29 

45 

52 

70 

S3DSrS4DS 

44 

98 

152 

175 

232 

S60S 

62 

138 

211 

243 

321 

S7US 

102 

233 

358 

413 

547 

S7CS 

102 

230 

354 

409 

541 

Site  ID 

Maximum  Annua!  24-hr  Volume  (cubic  decameters) 

2-year 

10-year 

50-year 

100-year 

500-year 

S2DS 

2,100 

4,200 

6,500 

7,400 

9,800 

S3DS 

2.700 

5,eoo 

8,300 

£.800 

13,000 

S4DS 

1,100 

2.300 

3,600 

4,200 

5,5C0 

S3DStS4DS 

3,700 

7,700 

11,900 

13.800 

18,300 

S6DS 

S.OCO 

10,500 

10,000 

18,400 

24,300 

S7US 

8.300 

17,900 

27,600 

31.800 

42.000 

3703 

8,300 

17,900 

27,500 

31,700 

41.SOO 

Site  !  D 

Maximum  Annual  3-Day  Volume  (cubic  decameters) 

2-year 

1  C-year 

50-year 

100-year 

500-year 

S2DS 

4,500 

9,400 

14,500 

16,800 

22,230 

S303 

5,900 

12,200 

19,100 

22,100 

29,400 

S4DS 

2,500 

5,200 

3,400 

9,800 

13,100 

S3DS+S40S 

8,200 

17.100 

27,000 

31,400 

41,900 

SSOS 

10,200 

21,400 

33,400 

38,600 

51,200 

S7US 

17,600 

•33.000 

59,60o 

69,100 

52,000 

S70S 

17,600 

37,900 

59,600 

69,100 

92,000 
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Sub-Bnsln  Parameters 
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Plow  Routing  Parmatcrs 

Travel  Tlmo  of  Flood  Wavo  (lira)  N/A 

Number  ofSub-Reaclics  (1-lir  time  step)  N/A 


Table  3.5.7 

Simulation  Output  Sites 
Abnean  Busin 


Table  3.5.8 

Computed  Maximum  Annual  Peak  Discharge  and  Volume  Frequency  Data 
Aba  can  Basin 


Maximum  Annual  Peak  Discharge  (cms) 


Site  ID 

2-year 

10-year 

50-year 

100-year 

AIDS 

10 

23 

35 

40 

A3DS 

31 

69 

106 

122 

A5DS 

41 

90 

138 

159 

Maximum  Annual  24-hr  Volume  (dam3) 


Site  ID 

2-year 

10-year 

50-year 

100-year 

AIDS 

800 

1,600 

2,400 

2,800 

A3DS 

2,500 

5,100 

7,800 

9,000 

A5DS- 

3,100 

6,500 

10,000 

11,500 

Maximum  Annual  3-Day  Volume  (dam3) 


Site  ID 

2-year 

10-year 

50-year 

100-year 

AIDS 

1300 

2,800 

4300 

5,000 

A3DS 

4,900 

10,300 

16300 

18,800 

A5DS 

6,300 

13,100 

21,000 

23300 

500-year 

52 

161 

211 


500-year 

3,600 

11,900 

15,300 


500-year 

6,700 

25,000 

21,900 


Table  3.5.9 
Computed  Flow  Duration  Curves 
Abacan  Basin 

Data  in  cubic  meters  per  second 


%  OF  TIME  AIDS  A3DS  A5DS 

EXCEEDED 

100  0.G  0.0  0.0 

50  0.2  0.6  0.8 

20  0.5  2.1  2.6 

10  0.9  3.4  4.3 

5  1.7  5.1  6.3 

2  3.0  8.6  10.7 

1  4.3  12.7  15.8 

0-5  6.0  18.0  22.4 

0-2  9.3  29.3  36.6 

0.1  12.7  40.3  50.8 

0.05  16.1  51.3  65.2 

0.02  19.7  62.9  80.2 
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Clntk  Sloinue  CoefndoiK  (Ills)  37.5  45,0  20,0  45.0  50.0  17  5  62.5  55.0  65  0 

Unirorm  Iiiflltrnliott  Loss  (mm/lir)  333333333 

Basedow  (ems)  2,7  3.7  1.8  2,0  2.5  1.3  2.7  1,7  0.6 
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Flow  Fouling  Pninmclcw 

Tfdvel Time  of  Flood  Wove  (Iks)  I.*  N/A  1.1  N/A 

Number  of  Stib'Rcndics  ( I  -hr  lime  slep)  2  N/A  1  N/A 


Table  3.5.12 

Computed  Maximum  Annual  Peak  Discharge  and  Volume  Frequency  Data 
0*  Donnell  Basin 


Site  ID 

OIDS 

05US 

OSDS 

07DS 

09U3 

OlOtJS 

O10DS 

015U5 

016DS 

017DS 

018DS 


Site  ID 

OIDS 

OSOS 

OSDS 

OTDS 

OSDS 

OlODS 

010DS 

015DS 

016DS 

OlTDS 

Q18DS 


Site© 

OIDS 

05US 

OSDS 

07DS 

09US 

OlODS 

OlODS 

015US 

016DS 

017DS 

018DS 


Maximum  Annual  Peal:  Discharge  (cms) 


2-year 

10-year 

50-year 

100-year 

500-year 

31 

65 

97 

110 

144 

121 

258 

389 

446 

585 

150 

321 

490 

563 

741 

47 

102 

154 

177 

233 

63 

137 

210 

241 

318 

218 

468 

714 

822 

1,084 

223 

477 

730 

841 

UlO 

264 

546 

815 

932 

1316 

392 

799 

1301 

1,377 

1,802 

409 

830 

1351 

1,435 

1,881 

616  1,273  1,933 

Maximum  Annual  24-hr  Volume  (dam3) 

2321 

2,922 

2-year 

10-year 

50-year 

100-year 

500-year 

2,600 

5,000 

7,400 

8300 

11.000 

9300 

19,800 

29,700 

34,000 

44300 

12300 

25,000 

38,000 

43.700 

57,400 

3,81. 

7,600 

11300 

13300 

17,400 

5,100 

10300 

16,000 

18,400 

24300 

17,800 

36300 

55,700 

64,000 

84300 

18,200 

37,600 

57,400 

66,100 

87300 

21300 

42,600 

63,100 

72,000 

93,700 

32,400 

64,200 

96,100 

110.000 

143,600 

33,700 

67300 

101,100 

115,800 

151,500 

50,600  102,600  155,300 

Maximum  Annual  3-Day  Volume  (dam?) 

178,400 

234300 

2-year 

10-year 

50-year 

100-year 

500-year 

5,600 

11.000 

16,400 

18,700 

24300 

20,900 

42,400 

C4300 

73,800 

56,800 

26,300 

54,100 

£3,600 

96,400 

127300 

7,800 

16,000 

24300 

28300 

37300 

10,900 

22,700 

35300 

40,600 

53,700 

38,300 

79,400 

123,100 

142300 

188300 

39,600 

82300 

128300 

148,400 

196,900 

48,200 

95300 

141,900 

162,000 

211,100 

74,400 

149,000 

224,900 

258,000 

337,600 

78^00 

157,600 

239,100 

274,800  ' 

360,400 

115300 

236,100 

362,100 

417300 

549,800 

Table  3.5.13 
Computed  Flow  Duration  Curves 
O’Donnell  Basin 
Data  in  cubic  meters  per  second 


%  OF  TIME 
EXCEEDED 

OIDS 

05US 

OSDS 

07DS 

09US 

OlOUS 

100 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

50 

0.8 

3.0 

4.0 

1.2 

1.8 

6.0 

20 

2.4 

9.4 

12.1 

3.8 

5.3 

18.2 

10 

3.9 

15.6 

20.1 

6.3 

8.8 

30.3 

5 

5.7 

22.4 

28.6 

8.9 

12.4 

42.8 

2 

9.5 

36.9 

46.8 

14.6 

19.9 

69.4 

1 

13.7 

53.0 

66.5 

20.7 

28.0 

98.0 

0.5 

19.0 

73.2 

91.2 

28.3 

37.9 

133.4 

0.2 

28.9 

111.3 

138.5 

42.6 

56.8 

201.5 

0.1 

38.5 

149.5 

187.6 

57.5 

77.1 

273.4 

0.05 

49.3 

193.0 

243.2 

74.3 

100.1 

354.9 

0.02 

62.5 

245.3 

309.8 

94.5 

127.5 

452.2 

%  OF  TIME 
EXCEEDED 

O10DS 

015US 

016DS 

017DS 

018DS 

100 

0.0 

0.0 

0.0 

.0 

0.0 

50 

6.4 

7.4 

12.6 

13.6 

20.0 

20 

19.2 

22.5 

37.1 

39.8 

58.7 

10 

32.1 

37.5 

61.9 

66.7 

98.2 

5 

45.0 

52.9 

85.5 

91.5 

135.2 

2 

72.5 

85.7 

135.5 

144.1 

213.6 

1 

101.7 

120.8 

186.9 

197.5 

294.0 

0.5 

137.5 

164.3 

248.5 

260.9 

389.8 

0.2 

206.5 

244.0 

363.0 

379.0 

571.1 

0.1 

280.8 

324.5 

486.2 

508.5 

771.1 

0.05 

364.8 

416.3 

626.3 

655.7 

998.0 

0.02 

465.1 

527.5 

795.2 

832.9 

1270.0 
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Sub-Basin  Parameters 
Gumaln/Pornc  Basin 
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Flow  Routing  Parameters 

Travel  Time  of  Flood  Wave  (lus)  N/A  0.9  N/A  0.8  0.5  N/A  N/A 

Number  of  Sub-Readies  (I -hr  lime  slop)  N/A  I  N/A  I  I  N/A  N/A 


Gtimain/Porac  Basin 


O  ci  co  O 

s  2  S  "  " 


d  T  'f'  <  < 
2  °  “  2  2 


N  <n  «  wi 


~  “=  2  Z 


or-*  oo  o 

£  ^  do 


r~  v*  n  r~ 
o'  t-'  o’ 


conn 
—  o'  —• 


oo  o  n  -* 
o  d  — 


oo  O  n  *o 
n  v*  — * 


n  ^  o  «*>  i 

«  «5  2  ^  ' 


r"*  O  O  <N  O 

o’  —  g  d  r: 


—  *o  m 

-  K  « 


2  *  -= 

2  2^  ®  O 
S  C  U-  c  ' 
a  ;  -  o  1 

n  jt  w  — - 
**•  W  TT 

I  s  §!. 

«  <  J  w  < 


2 

u 

|  <J  - 
-  2  c  r«  — 

|loo 

;  2  n  c 

O  CO  w  w 

u  —  a  <3 

,  =  3  ffl  OQ 

~  E  -o  x> 

!  i2  3  5  = 

.5  Z  *2  ™ 

;  <5  ■<  c  ra 

;  £  o  -  « 

°§EE 

c  c  o  o 

;  o  2  w  w 


O  c  S 

c^9 

=  e  c 

£  o  .2 

2  I  3  2  -. 

u  o  C)  —  ^ 

I  §  g  !  5 

5  S  e  S 

J=  O  W  c  O 

£  H  ■£  §  ? 

fc  =  .2  =  S 

§  t-  u  a  a 


Flow  Routing  Parameters 

Travel  Time  of  Flood  Wave  (lirs)  1.0  0.6  1.1  N/A  1.3  N/A 

NumbcrofSub-Rcaclics(l-hrtlmestcp)  1  1  1  N/A  1  N/A 
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Sub-Basin  Parameters 
Gumain/Porac  Basin 
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Flow  Routing  Parameters 

Travel  Time  of  Flood  Wave  (brs)  N/A 

Number  of  Sub-Readies  (l-lir  lime  slop)  N/A 


Table  3.5.16 

Computed  Maximum  Annual  Peak  Discharge  and  Volume  Frequency  Data 
Gumam/Porac  Basrn 


Maximum  Annual  peak  Discharge  (ems) 


Sue  ID 

2-vear 

10-year 

50-year 

100-veai 

500-vear 

G9DS 

220 

470 

700 

801 

1,048 

G1CDS 

228 

485 

723 

828 

1.0S3 

G12DS 

43 

95 

145 

166 

219 

G17DS 

128 

283 

430 

494 

651 

GISDS 

146 

320 

488 

561 

740 

G19US 

367 

">94 

1.196 

1.371 

1,801 

G19DS 

366 

790 

1,191 

1.366 

1,793 

Maximum  Annual  24-hr  Volume  (dam5) 


Site  ID 

2-year 

10-year 

50-year 

100-year 

500-year 

G9DS 

16,200 

31,800 

47,200 

o 

o 

•n 

70.300 

G10DS 

16,900 

33,300 

49,500 

56,600 

73,800 

G12DS 

3,400 

7,000 

10,600 

12^00 

16,000 

G17DS 

10,000 

20,700 

31,400 

36,000 

47,400 

G18DS 

11,500 

23,400 

36,300 

41,700 

54,900 

G19US 

28,000 

56,400 

84,800 

97,200 

127.300 

G19DS 

28,000 

56,300 

84,700 

97.000 

127,100 

Maximum  Annual  3-Day  Volume  (dam3) 


Site  ID 

2-year 

10-year 

50-year 

100-year 

500-year 

G9DS 

28,600 

57,300 

85,700 

98.000 

127,700 

G10DS 

30,500 

61,100 

91,800 

105.100 

137,200 

G12DS 

6,600 

13,800 

21,400 

24.600 

32,600 

G17DS 

19,400 

40,300 

62.500 

72.200 

95,600 

G18DS 

24,000 

48,000 

74,900 

S6.500 

114,600 

G19US 

52,800 

107.700 

164.500 

1S9.300 

248,900 

G19DS 

52,700 

107,000 

164,500 

1S9.300 

24S.900 

Table  3.5.17 

Computed  Flow  Duration  Curves 
Gumain/Porac  Basin 
Data  in  cubic  meters  per  second 


G17DS 

0.0 

2.7 

8.7 

14.3 

21.0 

35.3 

51.7 

72.7 
117.7 
162.6 
207.9 
255.3 


%  OF  TIME 

G18DS 

G19US 

G19DS 

EXCEEDED 

100 

0.0 

0.0 

0.0 

50 

3.4 

7.5 

7.5 

20 

10.7 

24.0 

24.0 

10 

17.7 

39.6 

39.6 

5 

25.6 

58.4 

58.4 

2 

42.4 

98.4 

98.4 

1 

61.1 

144.2 

144.2 

0.5 

84.7 

203.2 

203.2 

0.2 

129.9 

327.1 

326.9 

0.1 

175.7 

447.1 

446.5 

0.05 

227.6 

568.1 

567.2 

0.02 

289.7 

695.6 

694.3 

%  OF  TIME 

G9DS 

G10DS 

G12DS 

EXCEEDED 

100 

0.0 

0.0 

0.0 

50 

3.6 

4.1 

0.9 

20 

11.6 

13.0 

2.9 

10 

21.0 

23.6 

4.8 

5 

36.6 

40.7 

7.1 

2 

64.6 

71.0 

12.0 

1 

92.0 

99.8 

17.5 

0.5 

125.3 

133.5 

24.7 

0.2 

187.9 

196.2 

40.0 

0.1 

254.3 

266.0 

55.1 

0.05 

321.3 

336.3 

70.3 

0.02 

392.2 

410.7 

86.3 

Q. 


P70S  Yes  PoHero  River  above  confluence  with  Guagua  River  0  76.7  5.5  3,150 


Table  3.5.20 


Computed  Maximum  Annual  Peak  Discharge  and  Volume  Frequency  Data 
Pasig-Potrero  Basin 


Site  ID 

Maximum 

2-year 

10-year 

POOS 

55 

112 

PlDS 

69 

143 

P2DS 

8.5 

19 

P3DS 

5.9 

13 

P4US 

77 

161 

P4DS 

79 

166 

reus 

84 

178 

P5DS 

87 

186 

P7DS 

90 

195 

Annual  Peak  Discharge  (cms) 


50-year 

100-year 

500-yr 

166 

189 

246 

211 

241 

314 

28 

32 

42 

21 

24 

32 

238 

272 

355 

247 

282 

369 

265 

304 

397 

279 

319 

419 

294 

33"’ 

444 

Site  ID 


PODS 

P1DS 

P2DS 

P3DS 

P4US 

P4PS 

pais 

P5DS 

P7DS 


Maximum  Annual  24-hr  Volume  •  OIC  decameters) 


2-year 

4.100 

5.100 
500 
400 

5.700 
5.800 
6.200 
6.400 

6.700 


10-year 

7.700 

9.800 

1.100 

900 

10,900 

11.300 

12.100 

12.800 

13,600 


50-year 

11.200 

14.460 

1.700 

1.500 

16,000 

16.600 

18,000 

19.100 

20.400 


100-year 

12.700 

16.400 
1.900 

1.700 
18.300 
19.000 
20.500 
21.900 

23.400 


500-yr 

16.500 
21.300 

2.500 

2.200 

23.700 

24.700 

26.700 
28.600 
30.800 


Site  ID 


PODS 

P1DS 

P2DS 

P3DS 

P4US 

P4DS 

P5US 

P5DS 

pros 


Maximum  Annual  3-Day  Volume  (cubic  decameters) 


2-year 


10-year 


50-year  100-year  5C0-yr 


7.400 

13.900 

9.200 

17.700 

800 

1.600 

800 

1.700 

10,000 

19.300 

10.300 

19,900 

10.900 

21.400 

11.500 

22.800 

12.100 

24.500 

20.200 

23.000 

26.000 

29,600 

2.500 

2.900 

2.700 

3.100 

28.500 

32.500 

29.600 

33,800 

32.000 

36.600 

34.400 

39.500 

37.400 

43.000 

29.600 

38.400 

3.900 

4.200 

42.200 
44.000 
47.800 
51.900 
57.000 


Table  3.5.22  (sheet  1  o!  2) 

Sub-Unsin  Parameters 
Snnto  Tonmt  Dnsin 


['tow  Uonliiti*  Parameters 

TinvcITimc  ofPlood  Wave  (brs)  N/A  0.0  0.7  N/A  N/A  N/A 

Nimtbcr  of  Sub-Haicbcs/ Mir  lime  step)  N/A  l  I  N/A  N/A  N/A 


Table  3.5.24 

Computed  Maximum  Annual  Peak  Discharge  and  Volume  Frequency  Data 
Santo  Tomas  Basin 


Maximum  Annual  Peak  Discharge  (cms) 


Site  ID  - 

2-year 

10-year 

50-year 

100-year 

T2DS 

156 

320 

472 

538 

T5DS 

169 

347 

512 

584 

T7DS 

84 

172 

259 

296 

T9US 

224 

457 

684 

783 

T10DS 

285 

578 

867 

993 

Maximum  Annual  24-hr  Volume  (dam3) 


Site  ID 

2-year 

10-year 

50-year 

100-year 

TZDS 

12,300 

23,200 

34,000 

38,600 

T3DS 

13.400 

25,400 

37,200 

42,400 

T7DS 

7,000 

14,500 

21,800 

25.000 

TODS 

18,600 

36,800 

54,700 

62^00 

HODS 

23,600 

47,000 

70,100 

80.200 

Maximum  Annual  3-Day  Volume  (dam3) 


Site  ID 

2-year 

10-year 

50-year 

100-year 

12DS 

24,100 

45,700 

66,900 

76,000 

T3DS 

26,400 

50,500 

74,000 

84,200 

T7DS 

18,100 

37,600 

56,600 

64,900 

T9US 

42,100 

84,200 

125,300 

143,200 

TIODS 

54,300 

109,300 

163,500 

187,000 

500-year 

699 

760 

384 

1.023 

U97 


500-year 

50.000 

54,900 

32p00 

81,300 

104,400 


500-year 

98.200 

108.800 

84,600 

186,100 

243,500 


Table  3.5.25 
Computed  Flow  Duiation  Cur-'es 
Santo  Tomas  Basin 
Data  in  cubic  meters  per  second 


;  OF  TIME 

T2DS 

T3DS 

T7DS 

T9US 

EXCEEDED 

100 

0.0 

0.0 

0.0 

0.0 

50 

2.9 

3.3 

3.1 

6.3 

20 

9.3 

10.7 

9.0 

19.1 

10 

16.9 

17.6 

15.1 

31.8 

5 

29.2 

26.4 

20.4 

44.9 

2 

51.1 

45.2 

31.6 

72.7 

1 

71.9 

67.2 

42.7 

102  -  6 

0.5 

96.6 

95.8 

55 . 4 

139.4 

0.2 

140.9 

153.7 

79.8 

208.1 

0.1 

188.2 

205.6 

103.4 

278.7 

0.05 

235.7 

257.8 

140.8 

358.9 

0.02 

286.6 

313.5 

179.5 

455.7 

TIODS 

O.C 

3.6 

25.6 

42.7 
S3. 6 
55.4 

133.2 

179.2 
265  • 
355. 
458. 
5S2 . 
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Table  3.5,27 
Simulation  Output  Sites 


BIOUS  No  Bucao  River  :it  sticamflow  cnco  VV093A  15  6111  4.8  4,200  50  9 


Tabla  3.5.28 

Computed  Maximum  Annual  Peak  Discharge  and  Volume  Frequency  Data 
Bacao  Basin 


Maximum  Annual  Peak  Discharge  tons) 


SiicID 

2-year 

10-ytar 

50-year 

100->tar 

500-ytar 

BIDS 

86 

175 

260 

297 

387 

B2DS 

63 

132 

195 

223 

290 

B3DS 

15 

32 

48 

55 

72 

B5US 

200 

421 

626 

716 

933 

B5DS 

283 

592 

881 

1,007 

1,314 

B6DS 

45 

93 

139 

160 

208 

B7DS 

26 

54 

80 

91 

115 

B13DS 

376 

7S7 

1,173 

1,342 

1.752 

B14US 

643 

1,345 

2,007 

2,295 

2.SS8 

B16US 

780 

1,623 

•>  <20 

2.76S 

3,613 

B16DS 

834 

1,737 

2,591 

2,963 

3,869 

Maximum  Annual  24-hj  Volume  (dam5) 


Sire  ID 

2-year 

10-year 

SO-ytar 

100-year 

500-year 

BIDS 

7,100 

13,900 

20,600 

23,400 

30.400 

B2DS 

5300 

10,100 

14.900 

17.000 

22.000 

B3DS 

1300 

2300 

3,600 

4,100 

5,400 

B5US 

16.400 

32300 

47,500 

54300 

70,400 

B5DS 

23300 

45,800 

67,800 

77,300 

100,500 

B6DS 

3,700 

7,300 

10.900 

12300 

16300 

B7DS 

2,100 

4,100 

6,000 

6.900 

8,900 

B13DS 

31,000 

61,200 

90,700 

103,500 

134,700 

B14US 

53,000 

104,700 

155300 

177300 

230,600 

B16US 

64,100 

126300 

187,100 

213,600 

277,800 

B16DS 

68.600 

135300 

200,400 

228,700 

297,600 

Maximum  Annual  3-Day  Volume  (dam3) 


Site  ID 

2-year 

10-ytar 

50-year 

100-year 

500-year 

BIDS 

16,100 

32,000 

47,500 

54,300 

70,500 

B2DS 

11.100 

21,800 

32300 

36,800 

47,700 

B3DS 

2,600 

5,100 

7,600 

8,700 

11,300 

B5US 

34,800 

69300 

103,000 

117,600 

152,800 

B5DS 

49,800 

99,400 

147,800 

168,900 

219.600 

B6DS 

8300 

16,500 

24,700 

28300 

36,800 

B7DS 

4,300 

8,600 

12,700 

14,500 

18,800 

B13DS 

66,900 

134,100 

199,800 

228,400 

297300 

EMUS 

114,000 

229300 

341,700 

390,600 

508300 

BlfUS 

137,500 

275300 

409,800 

468,400 

609,400 

B16DS 

147,100 

294,700 

439,000 

501,800 

652,900 

Table  3.5.29 

Computed  Flow  Duration  Curves 
Bueno  Basin 

Data  in  cubic  meters  per  second 


*  OF  TIME 

BIDS 

B2DS 

B3DS 

B5US 

BSDS 

B6DS 

EXCEEDED 

100 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

5C 

2.3 

1.5 

0.4 

4.9 

7.1 

1.2 

20 

7.1 

4.7 

1.1 

15.3 

22.1 

3.6 

10 

11.8 

7.7 

1.9 

25.4 

36.6 

5.9 

i) 

16.7 

11.2 

2.7 

36.6 

52.6 

8.5 

2 

27.2 

18.8 

4.5 

60.5 

86.8 

13.9 

1 

38.7 

27.3 

6.  £ 

87.2 

124.6 

19.9 

0.5 

53.0 

73-2 

9.2 

120.9 

172.3 

27.4 

0.2 

79.4 

58.2 

14.0 

183.4 

261.0 

41.5 

0.1 

10S.8 

77.2 

18.7 

244.4 

348.1 

55.6 

0.05 

135.9 

98.3 

23.9 

314.0 

447.4 

71.6 

0.C2 

172.3 

125.3 

30.4 

398.1 

567.4 

91.0 

%  OF  TIME 

B7DS 

B13DS 

314US 

B16US 

B16DS 

EXCEEDED 

100 

0.0 

0.0 

0.0 

0.0 

0.0 

50 

0.6 

9.6 

17.0 

20.4 

21.5 

20 

1.8 

29.8 

52.0 

62.6 

66.2 

10 

3.0 

49.4 

86.4 

103.9 

109.9 

5 

4.4 

70.8 

123.2 

148.3 

157.2 

2 

7.4 

116.4 

201.5 

243.0 

258.2 

1 

10.8 

166.6 

287  .2 

346.7 

369  ^ 

0.5 

15.2 

229.8 

394.2 

476.5 

508.7 

0.2 

24.3 

347.5 

593.6 

717.9 

768.0 

0.1 

32.8 

464.3 

793.7 

958.7 

1026.0 

0.05 

41.3 

597.3 

1C31.4 

1232.8 

1319.6 

0.02 

50.4 

757.9 

1296.2 

1564.0 

1674.3 

Table  3.5.30 


Table  3.5.31 

Simulation  Output  Sites 
Mntomn  Bnsin 


Table  3.5.32 

Computed  Maximum  Annual  Peak  Discharge  and  Volume  Frequency  Data 
Maloma  Basin 


Maximum  Annual  Peak  Discharge  (cms) 


Site  ID 

2-year 

10-year 

50-year 

100-year 

M1DS 

12 

26 

38 

44 

M4DS 

119 

244 

360 

411 

M5DS 

69 

143 

212 

242 

M6DS 

212 

437 

646 

738 

Maximum  Annual  24-hr  Volume  (da m3) 

Site  ID 

2-year 

10-year 

50-year 

100-year 

M1DS 

900 

1.800 

2,600 

3,000 

M4DS 

9,700 

18,700 

27,400 

31500 

M5DS 

5,500 

10,700 

15,700 

17,900 

M6DS 

17.100 

33,000 

48,400 

555100 

Maximum  Annual  3- 

Day  Volume  (dam3) 

Site  ID 

2-year 

10-year 

50-year 

100-year 

M1DS 

1,700 

3,300 

4,800 

5,500 

M4DS 

21,500 

41,700 

61,200 

69,700 

M5DS 

11,100 

21,700 

32,100 

36,600 

M6DS 

36,000 

70,000 

103500 

117,600 

500-year 

57 

534 

315 

960 


500-year 

3,900 

40,400 

23^100 

71,600 


500-year 

75100 

905100 

47,500 

152,300 


Table  3.5.33 

Computed  Flow  Duration  Curves 
Maloma  Basin 

Data  in  cubic  meiers  per  secerns 


;  OF  TIKE 

M1DS 

M4DS 

M5DS 

M6DS 

EXCEEDED 

100 

0.0 

0.0 

0.0 

0.0 

50 

0.2 

2.7 

1.4 

5.3 

20 

0.6 

3.5 

4.6 

16.3 

10 

1.2 

14.1 

7.6 

27.1 

5 

2.0 

20.6 

11.3 

38.9 

2 

3.6 

34.6 

19.1 

64.0 

1 

5.2 

50.5 

28.2 

91.7 

0.5 

7.1 

71.0 

39.9 

126.7 

0.2 

10.7 

112.4 

63.6 

190.4 

0.1 

14.4 

150.8 

85.5 

252.1 

0.05 

18.2 

189.4 

107.5 

322.6 

0.02 

22.2 

230.6 

131.0 

'08.2 

Table  3.5.34  (sheet  1  of  2) 

Estimated  Peak  Instantaneous  Discharges  and  Confidence  Limits  at  Hydrologic  Output  Sites 


Table  3.5.34  (sheet  2  of  2) 

Estimated  Peak  Instantaneous  Discharges  and  Confidence  Limits  at  Hydrologic  Output  Sites 


T10DS 


Table  3.5.35  (sheet  1  of  2) 

Estimated  Maximum  Annual  24-Hour  Volume  and  Confidence  Limits  at  Hydrologic  Output  Sites 


Table  3.5.35  (sheet  2  of  2) 

Estimated  Maximum  Annual  24-Hour  Volume  and  Confidence  Limits  at  Hydrologic  Output  Sites 


Table  3.5.36  (sheet  1  of  2) 

Estimated  Maximum  3-Day  Volume  and  Confidence  Limits  at  Hydrologic  Output  Sites 


Table  3.5.36  (sheet  2  of  2) 

Estimated  Maximum  3-Day  Volume  and  Confidence  Limits  at  Hydrologic  Output  Sites 


Table  3.5.37  (sheet  1  of  2) 

Flow  Duration  Curve,  5  Percent  Confidence  Limit  Data  at  Hydrologic  Sites  (cms) 


M6DS 


Table  3.5.37  (sheet  2  ot  2) 

Flow  Duration  Curve,  5  Percent  Confidence  Limit  Data  at  Hydrologic  Sites  (cms) 


T10DS 


TABLE  3.5.38  (sheet  1  of  2) 

Flow  Duration  Curve,  95  Percent  Confidence  Limit  Data  at  Hydrologic  Sites  (cms) 

Percent  of  Time  Exceeded 


MiaEEaimnEi 


Table  3.6.1 

Estimated  Depths  o!  Clearwater  flows 


Table  3.6.2 

Estimated  Depths  for  Sediment  +  Water  Flows 


Abacan 

Downslr.  cf  Expressway  Bx. 
BeTw  Expr.  and  FnerxJsJxp  B 
Upsbeam  of  Friendsfrp  Bx. 

Bamban 

SacotxaandMarmta 
Sacotxa  »ow  orty 


downskaam  site 


03-0.6  (0.5)  O.S-O.9  (0.7)  0  6-1.1  (03) 

03-0.9  (0-S)  05-12  (OS)  0.7-1. 4  (0.9) 

0.4-0.6  (03)  06-0  9  (0.7)  03-1.0  (03) 


0  7-13(03)  0.8-13(11) 

0  8-1.4  (1.0)  0.9-1. 6  (13) 

0  6-1.1  (1.0)  0.9-13  (1.1) 


Bucao 

downstaam  of  bodge 
upskeara  cf  bodge 

Gtxnein 

Ooenstaaro  of  me:  km  16 
Upskeam  o I  mer  km  16 


0*3  1-3  2.1  2.6  2.9  3  4 

0-3  0.7-33  (2.1)  13-43  (33)  23-S3  (43)  33-57  (4.7)  43-63  (5.6) 


143-16.0  0.3  0.5  0.6  0.6  0.7 

16.0-16.0  03-23  (1.0)  04-2.6  (13)  03-2.7(1.4)  0.6-28(13)  0.7-3.1(16) 


P*»g 

Downskeam  o(  lower  bridge 
Beta,  lower  and  Bypass  8r. 
Upsteam  of  Bypass  Bridge 


Ptrac 


Dcwmswarr.  cl  bndge 
Upskeam  c!  bodge 


1-2.3  13-13  (1.4)  13*2.4  (23)  23-23  (23)  2.5-31  (2.7)  23*33  (3.1) 

23-4.1  0.7-1 .4  (1.1)  1.0-2 1  (1.7)  13-2.6  (2.1)  1 4-23  (23)  1.7-33  (2.6) 

4. 1-4.6  03*1.0  (03)  13-1.6  (1.4)  1.6-23  (1.7)  1.7-21  (13)  23-2.4  (22) 

I 

19*21  This  react)  was  not  modeled  tor  bufcad  6ow  as  P*  aamated  concentjaoo  of  tie  tuiw  J 
tow  was  such  fas  N  rcsdmo  mure  woUd  be  a  norvNewtcran  tud. 

III! 


5-7  J  1.1-20  (1.6)  1. 4-3.0  (23)  13-3.6  (23)  j  1.9-33  P3)  21-4.4  (3.4) 


16-20  03-13  (13)  03-26  (13)  0.4-3.1  (1.4)  0.5-32  (13)  0.6-3.6  (1.7) 

20*22.7  j  03-24  (1.0)  0.7-33  (13)  1.0-3. 9  (1.8)  1  0-4  l  (20)  1.1-43  (23) 


0*1.6  0.7-12  (33)  13*1.6  (13)  1.1-13  (13)  1.2-20  (13)  13-2 2  (1.6) 

1 .6-3.6  03-1.6  (1.1)  13-22  (1.7)  1 7-24  (23)  1.8-23  (2.1)  2023  (24) 


kMnafTartacj  03-13  (03)  03-13  (13)  13*22  (13)  1.1-24  (13)  13-23  (13) 

town  o?  Tartac  03-20  (1.4)  13-3.0  (21)  20-3.7  (2.7)  22-43  (29)  304.7  (33) 

town  Cf  Tartar  0.7-27  (1.4)  13-4  0  (23)  13-53  (33}  23-5  9  (33)  3.5-73  (S3) 


Table  3.6.3 

Estimated  Velocities  tor  Clear-Water  Rows 


ALacan 

Downs*,  of  Express*??  0r. 
Betw.  Expc.  and  Fnendshp  6r. 
Upstream  of  Fnendshp  Br. 


Sacotxa  and  Mararb  % 
Sacotxa  Sow  orfy 


H-15  (17) 
14-75  (2.0) 
1-6-27  (2.0) 


17-2.3  (2.1) 
0.3-1  1  (07) 
2.0-3  3  (56) 


152  6  '2  4) 
1.9-61  (28) 
2.4-3  4  (2.9) 


50-58  (2.5) 
1.9-6  0  (2.5) 
55-36  (31) 


51-25  (2  7) 
2. 1-6.1  (35) 
27-3.9  (3.4) 


Bangs! 

downstream  see 
RKSdesde 


dwmsteam  of  bodge 
upsteam  of  bridge 

Guman 

Downstaam  of  mer  km  16 
Ups*e»n  of  nver  km  16 


R»a 

Downsaeam  ol  lower  fridge 
Betw.  tower  and  Bypass  Br 
Upsaaam  of  Bypass  Bridge 


Poeac 

Near  Fkridatfanca 
Pone 

Downsaaam  of  bndge 
Upsaeam  of  bndge 

Sank:  Tomas 

Downstream  of  Maaolool  8*. 
Upsteam  of  Macctccl  Bndge 


i>rfac 

Daensaaaa  cf  Aqucx  Br  jiown  of  Tar 
Bel*.  Aquino  fni  Agana  Br.  flown  ef  Tar 
Upstream  of  Ajana  Bridge  {town  cf  Tar 


2.6 

0.7-2 5  (15) 


16.0-18.0  I  15-51  (1.6) 


1-2  5  09*1.1  (0.9) 

2.3*4. 1  0.7-15  (1.0) 

4.1-45  0  9-19  (15) 

19-21  1.9-35  (25) 


5-7  1.1-50  (1.6) 


18-2C  1.5-45  (2.0) 

20-22.7  1.7-25  (2  4) 


3.5  4.0  4.2 

1  0-35  (1.4)  05-4.0  (1.6)  0.9-4. 1  (15) 


0.4  0-5  0.6 

1 -6-2.6  (25)  15-3.0  (25)  50-3.1  (56) 


09-2  0  (15) 
5.3-3  7  (2.1) 


08-55  (1.3) 
15-50  (50) 
0.5-25  (15) 


1.2-5  4  (15) 
1.0-1.6  (1.4) 
1.1-53  *1.7) 


14-5.4  (25) 
50-45  (3.1) 


05-54  (17) 
1.3-25  (17) 


1  I -15  (15) 
15-45  (58) 
0  8-11  flB) 


1 .4-1.6  (15) 
15-53  (1.6) 
15-53  (15) 


15-55  (2  5) 
56-5.0  (35) 


06-57  (!5) 
15-2  4  (2.0) 


15-4  1  (2.1) 
15-45  (3.3) 
10-3.7  (15) 


15*1.7  (1.6) 
15-54  (1.7) 
14-24  (50) 


56-55  (3.7) 
15-55  (25) 


0.7-58  (2  0) 
1.6-56  (5*) 


15-15  (25) 
2.0-SI  (34) 
10-56  (1.9; 


4.7 

05-4.4  (1.8) 


0.7 

52-3  4  (58) 


1.6-15  (1  7) 

15- 25  (50) 

16- 25  (52) 


18-5  7  (58) 
26-5.7  (4.0) 


05-13  (19) 
I  £-2  6  (53) 


1  4-4  6  (55) 
25-5  v  (56) 
fl 


Table  3.6.4 

Estimated  Velocities  for  Sediment  -t  Water  Rows 


Table  3.6.5  (sheet  1  of  2) 

Hydraulic  Modeling  Notes 


'AO  sediment  concentrations  are  by  total  volume. 

"  DTM  (Digital  Terrain  Model). 

—  DPWH  (Department  ot  Public  Works  and  Highways). 


Table  3.6.5  (sheet  2  of  2) 

Hydraulic  Modeling  Notes 


"A1  sediment  concantrabo  ns  are  by  total  volume. 

—  DTM  (Digital  7  enain  Model). 

*•*  DPWH  (Department  of  Public  Works  and  Highways). 
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Monthly  Flow  Analysis 
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Daily  Discharge  Duration  Curve 
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Daily  Discharge  Duration  Curve 
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Daily  Discharge  Duration  Curve 
Station  W084A;  Porac  River,  Del  Carmen 
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Daily  Discharge  Duration  Curve 
Station  W088A;  Colo  River,  San  Benito 
Basin  Area  »  ?6  sq.  1cm 


Percent  of  Time  Exceeded 


Daily  Discharge  Durttion  Curve 
Station  W092A;  Bagsit  River,  Dampai 
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Annual  Basin  lUmofr 
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Basin  Area  177  si|.  km 


Figure  2.4.67 
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Figure  3.5.25 


Computed  Flood  Hydrographs 
C'Donnell  Basin  at  Site  018DS 


Figure  3.5.27 


Figure  1 


Computed  Flood  Hydrographs 
Gumain/Porac  Basin  at  Site  G9DS 


puooas  jsc  sjsqsm  oiqno  'sSjcqosiQ 


Figure  3.5.29 


Computed  Flood  Hydrograplis 
Gumain/Porac  Basin  at  Site  G10DS 


puooas  jsd  otqno  'acXBqosiQ 


nMfTrrtfnmnnnr 


Computed  Flood  Hydrographs 
Gumain/Porac  Basin  at  Site  G12DS 


Time  in  Days 


Flood  Hydrographs 
!  Basin  al  Site  G17DS 


Computed  1 
Gumain/Porac 


puooas  jsd  sjaqam  otqno  'aSjaqosiQ 


Figure  3.5.33 
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Figure  3.5.42 
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Figure  3.5.52 
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